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Wc report the first measurements of inclusive W and Z boson cross sections times tlie correspond- 
ing leptonic branching ratios for pp collisions at =1.96 TeV based on the decays of the W and Z 
bosons into electrons and muons. The data were recorded with the CDF II detector at the Fermilab 
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Tevatron and correspond to an integrated luminosity of 72.0 ± 4.3 pb~^. We test e-/i lepton univer- 
sality in W decays by measuring the ratio of the W —> to W ^ ev cross sections and determine a 
value of 0.991 ± 0.004(sfaf.) ± 0.011{syst.) for the ratio of W — /—!/ couplings {qij,/ Qe)- Since there is 
no sign of non-universality, we combine our cross section measurements in the different lepton decay 
modes and obtain aw x Br{pp ^ W £u) = 2.749 ± 0.010(stai.) ± 0.053{syst.) ± 0.165{lum.) nb 
and cr^-/z X Br(pp^ 7*/Z ^ = 254.9 ± 3.3{stat.) ± 4:.6{syst.) ± 15.2{lum.) pb for dilep- 
ton pairs in the mass range between 66 GeV/c^ and 116 GeV/c^. We compute the ratio R of the 
W —* Iv to Z ~* U cross sections taking all correlations among channels into account and obtain 
R = 10.84 ± 0.15(stat.) ± 0.14(sj/sf.) including a correction for the virtual photon exchange 
component in our measured 'y* /Z U cross section. Based on the measured value of R, we 
extract values for the W leptonic branching ratio, Br{W — > tv) = 0.1082 ± 0.0022; the total 
width of the W boson, V{W) = 2092 ± 42 MeV; and the ratio of W and Z boson total widths, 
V{W) /V{Z) — 0.838 ± 0.017. In addition, we use our extracted value of r(W) whose value depends 
on various electroweak parameters and certain CKM matrix elements to constrain the Vcs CKM 
matrix element, \Vcs\ ^ 0.976± 0.030. 

PACS numbers: 13. 38. Be, 13.38.Dg, 14.70.Fm, 14.70. Hp, 13.85.Qk, 12.38.Qk, 12.15.Ji 



I. INTRODUCTION 

Measurements of the production cross sections for 
both W and Z bosons in high energy pp collisions are 
important tests of the Standard Model (SM) of parti- 
cle physics. At hadron coUiders the W and Z bosons 
can most easily be detected through their leptonic 
decay modes. This paper presents measurements of 
aw ■ BriW tv), az ■ Br{Z U), and their ratio, 

az ■ Br{Z U) ^ ' 

for I = e and [i based on 72.0 pb^^ of pp collision 
data collected in 2002-2003 by the upgraded CoUider 
Detector at Fermilab (CDF) at a center-of-mass en- 
ergy of 1.96 TeV. These measurements are also de- 
scribed in [l[ . These measurements provide a test of 
SM predictions for the W and Z boson production 
cross sections, aw and cr^, as well as a precise indi- 
rect measurement of the total decay width of the W 
boson, r(VF), within the framework of the SM. This 
analysis is sensitive to deviations in r(VF) from the 
SM predictions at the level of about 2%. We also use 
our results to extract the leptonic branching fraction, 
Br(W iiy), and the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element, Vcs- Finally, we test the lep- 
ton universality hypothesis for the couplings of the 
W boson to e and /i leptons. 

A. W/Z production and decay 

The W and Z bosons, together with the massless 
photon (7), compose the bosonic fields of the uni- 
fied electroweak theory proposed by Weinberg 0, 
Salam [1], and Glashow The W and Z bosons 



q(q') l(l') q(q') KH 





i(v) 



FIG. 1: Diagrams for production and leptonic decay of 
a vector boson V — W, Z at leading (upper left) and 
next-to- leading order (others). 



were discovered in 1983 using the UAl and UA2 de- 
tectors d, H, 0, Q which were designed and built for 
this very purpose. The transverse momentum (pr) 
distribution of the reconstructed leptons in W ^ 
events was used to determine the W mass, while the 
Z mass was determined by directly reconstructing the 
invariant mass of dilepton pairs in Z i£ events. 

Present experimental measurements of electroweak 
parameters including vector boson masses and decay 
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widths are precise enough to provide tests of Quan- 
tum Chromodynamics (QCD) and the electroweak 
part of the Standard Model beyond leading order. 
These precise measurements not only test the elec- 
troweak theory but also provide possible windows to 
sectors of the theory at mass scales higher than those 
directly observable at current accelerator energies. 
These sectors enter into the electroweak observables 
through radiative corrections. While the parameters 
of the Z boson have been well studied , the proper- 
ties of the charged current carriers, the W bosons, are 
known with less precision. In hadron-antihadron col- 
lisions the W and Z are predominantly produced via 
the processes illustrated in Fig. [TJ The production of 
VV ~^ 7*/-^ where a quark in one hadron annihilates 
with an antiquark in the other hadron to produce 
the resulting vector boson is often referred to as the 
Drell-Yan production process. 

Calculations of the total production cross sec- 
tions for W and Z bosons incorporate parton cross 
sections, parton distribution functions, higher-order 
QCD effects, and factors for the couplings of the dif- 
ferent quarks and antiquarks to the W and Z bosons. 
Beyond the leading order Born processes, a vector bo- 
son V can also be produced by q{q)g interactions, so 
the Parton Distribution Functions (PDFs) of the pro- 
ton and antiproton play an important role at higher 
orders. Theoretical calculations of the W and Z pro- 
duction cross sections have been carried out in next- 
to-leading order (NLO) ELIiIL and next-to-next- 
to-leading order (NNLO) [l3,lli, [H [H, The 
NLO and NNLO computations used in this article are 
in the modified minimal-subtraction (MS) [l^, [l^ 
renormalization prescription framework. The full or- 
der ttg calculation has been made and includes fi- 
nal states containing the vector boson V and up to 
two additional partons. The two- loop splitting func- 
tion is used and the running of as includes thresh- 
olds for heavy flavors. The NLO cross section is ~ 
23^0 larger than the Born-level cross section, and the 
NNLO cross section is an additional ~ 3% higher. 
The main contribution to the calculated cross section 
is from qq interactions. The contribution of q{q)g in- 
teractions to the calculated cross section is negative 
at the Tevatron collision energy. 

The decay modes of the W boson are W ^ Iv (1 = 
e, /X, and r) and q(f , where the main modes ud, its, 
cs and cd have branching ratios proportional to their 
corresponding CKM matrix elements. The measured 
value for the branching fraction of the three combined 
leptonic modes is 32.0 ± 0.4% ^205, where the remain- 
ing fraction is assigned to the hadronic decay modes. 
The partial width into fermion pairs is calculated at 



lowest order to be [l^l 

ro(W^ ^ ff') = iFfff iVcGFM3,/(6\/2^), (2) 

where Vg' is the corresponding CKM matrix element 
for quark pairs or one for leptons. Mw is the W 
boson mass and Gp is the Fermi coupling constant. 
Nq is the corresponding color factor which is three 
for quarks and one for leptons. 

The expression for the partial decay widths into 
quark pairs also has an additional QCD correction 
due to vertex graphs involving giuon exchange and 
electroweak corrections due to next-to-leading order 
graphs which alter the effective coupling at the W- 
fermion vertex for all fermions. Within the context 
of the Standard Model, there are also vertex and 
bremsstrahlung corrections [2l[ that depend on the 
top quark and Higgs boson masses. The corrections 
can be summarized in the equation 

V{W ^ ff')sM = ro(VF ^ ff') 

x[l + ^v + '5w(o) (3) 

where 5w(o) is the correction to the width from loops 
at the M^-fermion vertex involving the Z boson or 
a SM Higgs boson, arises from the boson self- 
energies, and bp, is a correction required when the 
couplings are parametrized using the W mass and the 
value of Gp from muon decay measurements p^ . [23| . 
Since all of these corrections are small (~ 0.35%), the 
measurement of r(M^) is not very sensitive to these 
higher order effects. Higher order QCD corrections 
originating from quark mass effects are also small. 

B. Measurement of r(W) from the W and Z 
cross sections 

The width of the W boson can be extracted from 
the measurement of the ratio i?, which is defined in 
Eq. [1] This method was first proposed by Cabibbo 
in 1983 as a method to determine the number of 
light neutrino species [13] and has been adopted as 
a method to indirectly measure the branching ratio 
for the W ^ Iv decay mode. The ratio R can be 
expressed as 

aw V{W iv) r(Z) 

az r{z -> ei) r{w) ' ^ ' 

On the right hand side of Eq. 21 the ratio of the 
W and Z production cross sections can be calcu- 
lated from the boson couplings and knowledge of the 
proton structure. The Z boson total width, r(Z), 
and leptonic partial width, T(Z £i), have been 
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measured very precisely by the LEP experiments @ . 
With the measured value of R the branching ratio 
Br{W iv) ^ T{W -> iv)lT{W) can be extracted 
directly from Eq. [H The total width of the W bo- 
son, r(VF), can also be determined indirectly using 
the SM prediction for the partial width, V{W Iv). 
As shown in Eq. [51 r(VF) depends on electroweak pa- 
rameters and certain CKM matrix elements. We also 
use our measurement of the total W width to con- 
strain the associated sum over CKM matrix elements 
in the formula for r(M^) and derive an indirect value 
for Vca which is the least experimentally constrained 
element in the sum. Finally, the ratios of the muon 
and electron W cross section measurements are 

used to determine the ratios of the coupling constants 
of the W boson to the different lepton species, pro- 
viding a test of the lepton universality hypothesis. 
For reference. Table U provides a summary of previ- 
ous experimental results for ■ Br{W — > iv) and 
<Jz ■ Br{Z U) along with the measured values for 
R and the extracted values of r(W^). The most re- 
cent direct measurement of r(VF) obtained by LEP 
is 2.150 ± 0.091 GeV @. 

C. Overview of this measurement 

The signature of high transverse momentum lep- 
tons from W and Z decay is very distinctive in the 
environment of hadron collisions. As such, the de- 
cay of W and Z bosons into leptons provides a clean 
experimental measurement of their production rate. 
Experimentally, the cross sections times branching 
ratios are calculated from 



aw ■ Br{W -> £v) 



az ■ Br{Z ^ U) = 



Afobs 



ATbck 



Aw ■ ■ J ^dt 

^obs _ ^bck 

Az-ez- J Cdt' 



where and are the numbers of W 



(5) 



(6) 



bck 



and Z ^ II candidates observed in the data; iV^ 
and are the numbers of expected background 

events in the W and Z boson candidate samples; Aw 
and Az are the acceptances of the W and Z decays, 
defined as the fraction of these decays satisfying the 
geometric constraints of our detector and the kine- 
matic constraints of our selection criteria; tw and tz 
are the combined efficiencies for identifying W and 
Z decays falling within our acceptances; and / Ldt is 
the integrated luminosity of our data samples. 

In measuring the ratio of the cross sections some of 
the inputs and their experimental uncertainties can- 
cel. The strategy of this measurement is to select 



W and Z boson decays with one or both leptons (e 
or [l) falling within the central region of the CDF 
detector. This region is well instrumented and un- 
derstood and has good detection efficiencies for both 
lepton species. Using common lepton selection crite- 
ria (contributing to the factors tw and tz) for the W 
and Z channels has the great advantage of decreasing 
the systematic uncertainty in the measurement of R. 
The resulting smaller systematic uncertainty offsets 
the expected increase in statistical uncertainty orig- 
inating from the requirement of a common central 
lepton. For each lepton species, the selection criteria 
are optimized to obtain the least overall experimental 
uncertainty. 

The measurement of the ratio R is sensitive to new 
physics processes which change the W oi Z produc- 
tion cross sections or the W ^ iv branching ratio. 
The W ^ branching ratio could be directly af- 
fected by new decay modes of the W boson, such as 
supersymmetric decays that do not similarly couple 
to the Z boson. A new resonance at a higher mass 
scale that decays to or Z bosons may change the 
production cross sections. One example of a particle 
with a larger mass is the top quark at mt = 174.3 ± 
5.1 GeV/c^, which decays to aW boson and a bot- 
tom quark [20|. In pp collisions at ^/s = 1.8 TeV the 
production cross section for tt pairs is 6.5^i;4 pb [33l |. 
about 3000 times smaller than direct W boson pro- 
duction [2^ . The decays of ti pairs which result in 
the production of two W bosons should change the 
measured value of R by about 7 x 10""*, which is well 
below our sensitivity. The total width of the W bo- 
son can also get contributions from processes beyond 
the SM. For example, in supersymmetry, the decay 
W+ X^X*^ may be possible if the charginos and 
neutralinos are light \i% and so a precise measure- 
ment of r(W^) can constrain the properties of these 
particles. 



D. Outline of the paper 

This paper is organized as follows: in Sec. |TT] the 
CDF detector is described, with particular attention 
given to the subdetectors essential in the identifi- 
cation of charged leptons and the inference of neu- 
trinos. Section IIIII describes the data samples used 
in this analysis, and the selection of the W and Z 
candidate events is described in Sec. |TVl Section IVl 
describes the calculation of the geometric and kine- 
matic acceptances of our candidate samples, and the 
methods used to determine the efficiencies for iden- 
tifying events within our acceptances are presented 
in Sec. IVII The estimation of the contributions to 
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TABLE L Previous measurements of the W and Z production cross sections times branching ratios along with the 
measured values of R and the extracted values of r(V(^). 



Experiment 




Mode 


aw ■ Br{W eu) 


az ■ Br{Z -> U) 


R 


r(w^) 




(TeV) 




(nb) 


(Pb) 




(GeV) 


CDF(Run I) [25, 26, 27, 28, 29] 


1.80 


e 


2.49 ± 0.12 


231 ± 12 


10.90 ± 0.43 


2.064 ± 0.084 


D0(Run lA) [30] 


1.80 


e 


2.36 ± 0.15 


218 ± 16 






D0(Run lA) [30] 


1.80 




2.09 ± 0.25 


178 ± 31 






D0(Run lA) [30, 31] 


1.80 


e + fi 






10.90 ± 0.49 


2.044 ± 0.093 


D0(Run IB) [32] 


1.80 


e 


2.31 ± 0.11 


221 ± 11 


10.43 ± 0.27 


2.17 ± 0.07 



our candidate samples from background processes are 
discussed in Sec. IVIIl and finally the calculation of 
the cross sections along with the resulting value of 
R and other extracted quantities are summarized in 

Sec. lyml 



II. THE EXPERIMENTAL APPARATUS 

The data used for the measurements reported in 
this note were collected with the upgraded Collider 
Detector (CDF) [35] at the Fermilab Tevatron pp col- 
lider. Detector upgrades were made to accommodate 
the higher luminosities and new beam conditions re- 
sulting from concurrent upgrades to the Tevatron ac- 
celerator complex. In addition to the increases in 
luminosity, the pp center of mass energy was also in- 
creased from = 1.80 TeV to = 1.96 TeV. The 
relatively small change in beam energies leads to a 
substantial increase in the production cross sections 
for high-mass objects such as W/Z bosons (~ 9%) 
and top quark pairs (~ 30 %). We highlight the up- 
grades to the Run I detectors and electronics in the 
following sections. 



A. The CDF II Detector 

CDF is a general-purpose detector [H, [H, [s^l de- 
signed to detect particles produced in pp collisions. 
As illustrated in Fig. [21 the detector has a cylin- 
drical layout centered on the accelerator beamline. 
Tracking detectors are installed in the region directly 
around the interaction point to reconstruct charged- 
particle trajectories inside a 1.4 T uniform magnetic 
field (along the proton beam direction). The field 
is produced by a 5 m long superconducting solenoid 
located at the outer radius of the tracking region 
(1.5 m). Calorimeter modules are arranged in a 
projective tower geometry around the outside of the 
solenoid to provide energy measurements for both 




FIG. 2: Elevation view of half of the CDF Run II detector. 



charged and neutral particles. The outermost part 
of the detector consists of a series of drift chambers 
used to detect muons which are minimum-ionizing 
particles that typically pass through the calorimeter. 

The z-axis of the CDF coordinate system is defined 
to be along the direction of the incoming protons. A 
particle trajectory is then described by 9, the polar 
angle relative to the incoming proton beam; (j), the 
azimuthal angle about this beam axis; and zq, the 
intersection point of the particle trajectory with the 
beam axis. The pseudorapidity of a particle trajec- 
tory is defined as 77 = — ln(tan(6'/2)). The transverse 
momentum, px, is the component of the momentum 
projected on a plane perpendicular to the beam axis. 
Similarly, the transverse energy, Et, of a shower or 
an individual calorimeter tower is given by E ■ sin 9. 
The total transverse energy in an event is given by a 
sum over all calorimeter towers E^hi where E^ 
is the transverse energy measured in the zth tower 
and fii is the projection of the vector pointing from 
the event vertex to the ith calorimeter tower onto the 
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plane perpendicular to the beam axis (unit normal- 
ized). The vector sum of transverse energies mea- 
sured in the calorimeter is corrected to account for 
muons which deposit only a fraction of their energy 
in the calorimeter. The missing transverse energy in 
an event is the equal magnitude vector opposite to 
this vector sum of transverse energies. Fixed points 
on the detector are described using polar coordinates 
{r,(j),z) where r is the radial distance from the beam 
axis, (f> is the azimuthal direction about the beam 
axis, and z is the distance from the detector center 
in the direction along the beam axis. In some cases 
we also use a detector pseudorapidity variable, 77dot, 
to refer to fixed locations within the detector. This 
variable is based on the standard definition of pseu- 
dorapidity given above where the angle is redefined 
in the context of a fixed location as 9 = arctan(r/z). 
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FIG. 3: COT single-hit residual distribution obtained 
from Z fifi events. 



B. Tracking System 

All of the detectors in the inner tracking region 
have been replaced for Run II. The new silicon track- 
ing system consists of three concentric detectors lo- 
cated just outside the beam interaction region. In 
combination, these detectors provide high resolution 
tracking coverage out to | rj^ct I < 2. For the mea- 
surements presented here, silicon tracking informa- 
tion is incorporated solely to aid in the rejection of 
cosmic ray events from our muon samples. The rel- 
evant hit information comes from the Silicon Vertex 
Detector (SVX-II) [s^ which contains five layers of 
double-sided micro-strip detectors at radii of 2.4 to 
10.7 cm from the center of the detector. The SVX- 
II detector consists of three barrels divided into 12 
wedges in (/>. In total, the three barrels cover roughly 
45 cm along the z-axis on each side of the detector 
interaction point. 

The new open-cell drift chamber referred to as the 
Central Outer Tracker (COT) [H, sits directly 
outside of the silicon tracking detectors in the radial 
direction. The measured momenta and directions of 
the high px lepton candidates in our event samples 
are obtained from track reconstruction based solely 
on COT hit information. The chamber consists of 
eight superlayers of 310 cm length cells at radii be- 
tween 40 and 132 cm from the beam axis. Each su- 
perlayer contains 12 layers of sense wires strung be- 
tween alternating layers of potential wires. The wires 
in four of the superlayers (axial layers) are strung to 
be parallel to the beam axis, providing particle track 
reconstruction in the transverse plane. In the other 
four superlayers (stereo layers), the wires are strung 
at ± 2 degree angles with respect to the beam axis to 



allow also for particle tracking in the z-direction. The 
two superlayer types are alternated in the chamber 
within the eight radial layers starting with the inner- 
most stereo layer. The COT chamber has over 30,000 
readout channels, roughly five times the number in 
the Run I tracking chamber [4lj. Particles travers- 
ing the central region of the detector (| rjdct I <1) are 
expected to be measured by all eight superlayers. 

The COT is filled with a gas mixture of 5(Po argon 
and 50% ethane. This mixture was chosen to ensure a 
fast drift velocity (^ 50 //m/ns) compatible with the 
short interval between beam bunch crossings and the 
expected rise in instantaneous luminosity. The max- 
imum drift distance in the chamber is 0.88 cm corre- 
sponding to a drift time on the order of 200 ns. The 
single-hit resolution in the chamber has been studied 
using the high pT muon tracks in Z ^ nji candi- 
date events. The measured offset between the indi- 
vidual hits associated with these muons and the re- 
constructed path of the muon track is shown in Fig. [3] 
Based on this distribution, we measure a COT single- 
hit resolution of 180 /im. 

The solenoid produces a 1.4 T magnetic field inside 
the tracking volume that is uniform to 0.1 % in the 
region \z\ < 150 cm and \r\ < 150 cm. The transverse 
momentum of a reconstructed track, px (in GeV/c), 
is determined from px — 0.3qBrc where B (in T) is 
the magnetic field strength, the total particle charge 
is qe (e is the magnitude of the electron charge), and 
Tc (in m) is the measured radius of curvature of the 
track. The resolution of the COT track momentum 
measurement decreases for high px tracks which bend 
less in the magnetic field. The curvature resolution 
has been studied by comparing the inward and out- 
going track legs of reconstructed cosmic ray events. 
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FIG. 4: Distribution of the difference in curvature for 
the two tracks associated with a cosmic ray event as re- 
constructed by the COT, using cosmic ray data (top). 
Distribution of the E/p variable defined in Sec. IIVI for 
W ^ eu events (bottom). The mean and a are obtained 
from the Gaussian fit in the range 0.8 < E/p < 1.08. 



The difference in the measured curvature for the two 
track legs in these events is shown on the top of 
Fig. m We determine a COT curvature resolution of 
3.6 X 10 ~^cm~^, estimated from the a of this distri- 
bution divided by ^/2. This corresponds to a momen- 
tum resolution of (Jp^/PT — ^■'^ ^ 10 ^"^[GcV/c]^^. 
The COT track momentum resolution is also studied 
using the E/p distribution (see Sec. IIVP of electron 
candidates in W ev events. This distribution is 
shown on the bottom of Fig. 21 Since the COT track 
momentum resolution measurement is less precise at 
high pt than the corresponding calorimeter energy 
measurement, the Gaussian width of this distribution 
for 0.8 < E/p < 1.08 provides an additional measure 
of the curvature resolution. The resulting value is 
in good agreement with that obtained from studying 
cosmic ray events. 



C. Calorimeters 

Calorimeter modules used to measure the energy 
of both charged and neutral particles produced in 
pp collisions are arranged around the outer edges 
of the central tracking volume. These modules are 
sampling scintillator calorimeters with a tower based 
projective geometry. The inner electromagnetic sec- 
tions of each tower consist of lead sheets interspersed 
with scintillator, and the outer hadronic sections are 
composed of scintillator sandwiched between sheets 
of steel. The CDF calorimeter consists of two sec- 
tions: a central barrel calorimeter (| ?7dct | <1) and 
forward end plug calorimeters (1.1< | 77dot | <3.64). 
The scintillator planes in the central barrel lie paral- 
lel to the beam line, while those in the forward end 
plugs are arranged in the transverse direction. The 
central barrel consists of projective readout towers, 
each subtending 0.1 in rydot and 15° in (p. Each end 
plug also has projective readout towers, the sizes of 
which vary as a function of Tydet (0.1 in r^dot and 7.5° 
in at I Tydot | = 1-1 to 0.5 in 77dot and 15° in (p at 
I 77dct I = 3.64). 

The central barrel section of the CDF calorime- 
ter is unchanged from Run I. It consists of an in- 
ner electromagnetic (CEM) calorimeter and an outer 
hadronic (CHA) calorimeter (4^. The end- wall 
hadronic (WHA) calorimeter completes the cover- 
age of the central barrel calorimeter in the region 
0.6 < I 77dct I < 1-0 and provides additional forward 
coverage out to | rydot | = 1-3 As part of the CDF 
Run II upgrade, the original gas calorimetry of the 
end plug region (| r^dot | > 1-1) was replaced with scin- 
tillator plate calorimetry using scintillator tiles read 
out by wavelength shifting fibers embedded in the 
scintillator ^44 . i45] . The new design has an improved 
sampling fraction and reduces forward gaps that ex- 
isted in the old gas calorimeter system. The new plug 
electromagnetic (PEM) calorimeter provides cover- 
age in the 1.1 < | rydot | < 3.6 region and the new plug 
hadronic (PHA) calorimeter provides coverage in the 
1-3 < I 77dot I < 3.6 region Both the PEM and 

PHA incorporate the same polystyrene based scintil- 
lator and similar photomultiplier tubes used in the 
CEM. 

Calorimeter energy resolutions are measured us- 
ing test-beam data. The measured energy resolu- 
tions for electrons in the electromagnetic calorimeters 
are 14%/VE^ (CEM) and 16 ® 1% (PEM) 

j3S] where the units of Et and E are GeV. We 
also measure the single-particle (pion) energy reso- 
lution in the hadronic calorimeters to be 75 %/\fE 
(CHA), 80 (WHA), and 80 ® 5 % 

(PHA) [sHi ■ The energy resolution in the electromag- 
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netic calorimeters is also determined using Z ^ ee 
candidate events. The calorimeter energy scale is 
set so that the mean of the Gaussian fit to the di- 
electron invariant mass peak is 91.1 GeV/c^. This 
procedure results in a GEM energy resolution of 
13.5 %I\/Et © 1.5 %, in good agreement with the 
test-beam result [l^l- Jet energy resolution in the 
hadronic calorimeter sections 48 ] is determined using 
photon-jet balancing. In events in which a photon re- 
coils against a jet and no other activity is observed, 
the transverse energies associated with the two ob- 
jects must be equal and opposite. The photon en- 
ergy measured in the electromagnetic section of the 
calorimeter provides a reference point against which 
the energy deposition associated with the jet can be 
compared. The resolution of the large component of 
jet energy deposition in the hadronic calorimeters can 
be determined based on this comparison. The vast 
majority of hadronic particle showers are completely 
contained within the calorimeter. The combined lon- 
gitudinal depth of the central calorimeter module in 
interaction lengths is roughly 5.5 A and the equivalent 
depth in the plug modules is roughly 8.0 A. However, 
some small fraction of hadronic particle showers does 
leak out from the back end of the calorimeter, com- 
plicating muon identification. 

Proportional chambers (GES) are embedded in the 
electromagnetic section of the central barrel at a radi- 
ation length depth of roughly 6 Xq corresponding to 
the region of maximum shower intensity for electrons. 
These chambers are used to measure the profile of a 
shower and extract the location of the incident par- 
ticle within a given tower. The increased shower po- 
sition resolution provides additional selection criteria 
for electron candidates based on track-shower match- 
ing. The chambers, two per calorimeter wedge, uti- 
lize wires in the r-(^ view and cathode strips in the z 
view to determine the three-dimensional position of 
each shower. The resolution of the GES position mea- 
surement in r-(^ is roughly 0.2 cm. Each calorimeter 
module also has a second set of chambers (GPR) sit- 
uated on the front of the corresponding electromag- 
netic section which presamples each shower to pro- 
vide additional information useful in electron identi- 
fication and pion-photon separation. 

The first layer of the plug electromagnetic 
calorimeter is used as a preshower detector (PPR). Its 
scintillator is polyvinyltoluene-based, and it is twice 
as thick as the other sampling layers in the PEA/I. It 
has the same transverse segmentation as the PEM, 
but each scintillator tile in the PPR is read out in- 
dividually. The PEM also has a shower maximum 
detector (PES) embedded in it at a depth of ~ 6 Xq 
[4^ . The PES consists of two layers of 5 mm wide 



polyvinyltoluene-based scintillator strips, with each 
layer having a 45° crossing angle relative to the other. 
The PES provides coverage in the 1.1 < | rydot | < 3.5 
region. 



D. Muon detectors 

In order for a muon to pass through the calorimeter 
and into the most central portion of the GDF muon 
detector (| r/dct | < 0-6), it must have a minimum 
on the order of 1.4 GeV/c. In order to reach the outer 
portion of the central detector or the more forward 
detectors (0.6 < | 77dct | < 1-0), the muon is required 
to pass through an additional layer of steel absorber. 
Muons with a momentum above 3.0 GeV/c are essen- 
tially 100 % efficient for traversing the steel absorber 
over the entire solid angle of the combined muon de- 
tector coverage. The amount of energy deposited in 
the calorimeter by high muons produced in W and 
Z boson decays is observed to be Landau distributed 
with means of 0.3 GeV for deposits in the electromag- 
netic section and 2.0 GeV for those in the hadronic 
section. Reconstructed particle tracks in the GOT 
matched to minimum ionizing-like energy deposits in 
the calorimeter are treated as "stubless" muon candi- 
dates even in cases where the tracks are not matched 
with any hits in the muon detectors. The muon offline 
reconstruction forms stubs based on hit information 
in the muon detector and matches found stubs with 
the reconstructed GOT tracks to determine our high- 
est quality muon candidates. This final set of muon 
candidates includes only a small percentage of non- 
muon fakes originating from other hadronic particles 
that are not fully contained within the calorimeter 
(hadronic punchthrough) . Despite the fact that a 
non-negligible number of hadrons (on the order of 
1 in 220) pass through the entire calorimeter, the 
majority of those that enter the muon detector are 
absorbed in the filtering steel and do not produce 
associated hits in the outer sections of the detec- 
tors. Gonversely, "stubless" muon candidates include 
a substantially larger fraction of non-muon fakes, and 
the presence of additional physics objects (such as 
a second higher quality muon) associated with these 
candidates is typically required to increase the purity 
of the sample. 

The CDF muon detector is made up of four in- 
dependent detector systems outside the calorimeter 
modules. The Gentral Aluon Detector (GMU) [50|] is 
mounted directly around the outer edge of the central 
calorimeter module. The GMU is an original Run I 
detector component containing 2,304 single- wire drift 
chambers arranged in four concentric radial layers. 
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The drift chamber wires are strung parallel to the 
direction of the incoming beams, and wire pairs on 
layers 1 and 3 and layers 2 and 4 project radially back 
to the nominal beam interaction point, allowing for 
a coarse px measurement based on the difference in 
signal arrival times on the two wires within a pair. 
The CMU system provides symmetrical coverage in 
(j> in the central part of the detector (| ?7dot | < 0.6). 
The drift chambers have been upgraded to operate 
in proportional mode (in Run I these chambers were 
run in streamer mode). Operating in this mode re- 
duces the high voltage settings for the chambers and 
helps to prevent voltage sagging which is an issue due 
to the higher hit rates at Run II luminosities. Pre- 
cision position measurements in the (j) direction are 
made by converting signal arrival times into drift dis- 
tances in the plane orthogonal to the wire direction. 
The wires of cells in neighboring stacks are connected 
via resistive wires at the non-readout end of cells to 
also provide a coarse measurement of each hit posi- 
tion along the direction of the wire {z coordinate). 
The measurement is made by comparing time-over- 
threshold for the signals observed at the readout end 
of the two neighboring stacks. The maximum drift 
time within a CMU cell is 800 ns which is longer than 
the 396 ns spacing between bunch crossings in the ac- 
celerator. The ambiguity as to which beam crossing 
a particular CMU hit originates from is resolved in 
both the trigger and the offline reconstruction using 
timing information associated with a matched COT 
track and/or matching energy in the calorimeter. 

The Central Muon Upgrade Detector (CMP) and 
Central Muon Extension Detector (CMX) were also 
part of the CDF Run I configuration [5l|. The in- 
dividual wire drift chambers of these detectors are 
identical except for their lengths along the direction 
of the wire which is larger for CMP chambers. These 
drift cells are roughly a factor of two wider than 
those in the CMU detector resulting in a longer maxi- 
mum drift time of 1.8 /its. Matching scintillator detec- 
tors (CSPjCSX) installed on the outer edges of these 
systems can in principle provide timing information 
to resolve the three beam-crossing ambiguity arising 
from the long drift time. In practice, however, occu- 
pancies in these chambers are small enough at cur- 
rent luminosities to uniquely determine the appropri- 
ate beam-crossing from COT track matching. CSX 
timing information is used in the trigger to eliminate 
out-of-time hits from the beam halo associated with 
particle losses in the accelerator tunnel, but infor- 
mation from the scintillator systems is not currently 
utilized in muon reconstruction in this analysis. The 
CMP/CMX drift chambers are also run in propor- 
tional mode. The CMP chambers are arranged in a 



box-like structure around the outside of the CMU de- 
tector and an additional 3A of steel absorber which 
is sandwiched between the two detectors. The ad- 
ditional steel greatly reduces hadronic punchthrough 
into the CMP chambers and allows for cleaner muon 
identification. A total of 1,076 drift cells arranged in 
four staggered layers form the four-sided CMP struc- 
ture which provides additional coverage for the cen- 
tral part of the detector (| rj^et \ < 0.6) with variable 
coverage in </>. Drift cells in the CMX detector are ar- 
ranged in conical arrays of eight staggered layers to 
extend muon coverage up to | ij^et \ < 1-0. The par- 
tial overlap between drift tubes in the CMX conical 
arrangement allows for a rough hit position measure- 
ment in the z coordinate utilizing the different stereo 
angles of each cell with respect to the beam axis. 
The Run I configuration consisted of 1,536 drift cells 
arranged in four 120° sections providing coverage be- 
tween -45° to 75° and 105° to 225° in <j) on both ends 
of the detector. An additional 60° of CMX coverage 
on the bottom of the detector at both ends has been 
added for Run II, but these new components were still 
being commissioned in early running and are not uti- 
lized in the measurements reported here. The Bar- 
rel Muon Upgrade Detector (BMU) is another new 
addition for Run II which provides additional muon 
coverage in the regions 1.0 < | Ty^ot | < 1-5. This new 
detector system was also being commissioned in the 
initial part of Run II and is not used in these mea- 
surements. 



E. Cherenkov Luminosity Counters 

The small-angle Cherenkov Luminosity Counters 
(CLC) detector is used to measure the instantaneous 
and integrated luminosity of our data samples. This 
detector system is an additional Run II upgrade [s^ l 
that allows for high-precision luminosity measure- 
ments up to the highest expected instantaneous lu- 
minosities. 

The CLC consists of two modules installed around 
the beampipe at each end of the detector, which pro- 
vide coverage in the regions 3.6 < | rj^ct \ < 4.6. Each 
module consists of 48 long, conical gas Cherenkov 
counters pointing to the collision region. The coun- 
ters are arranged in three concentric layers of 16 coun- 
ters each, around the beam-pipe. The counters in the 
two outer layers are about 1.8 m and those in the in- 
ner layer are 1.1 m long. Each counter is made of 
highly reflective aluminized Mylar with a light col- 
lector that gathers the Cherenkov light into fast, ra- 
diation hard photomultiplier tubes with good ultravi- 
olet quantum efficiency. The modules are filled with 
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isobutane gas at about 22 psi which is an excellent 
radiator while having good ultraviolet transparency. 

The Cherenkov light cone half-angle, 9c, is 3.1° 
corresponding to a momentum threshold for light 
emission of 9.3 MeV/c for electrons and 2.6 GeV/c 
for pions. The expected number of photoelec- 
trons, -/Vpc, for a single counter is given by iVpo — 
No ■ L ■ sin^Oc where L is the distance tra- 
versed by the particle in the medium and No = 
370 cm-ieV-^ tco\{E)tAct{E)dE . The edot and EcoI 
terms are defined as the light detection and collec- 
tion efficiencies, respectively, and are functions of the 
energy E of the Cherenkov photon (in eV). Our 
design results in A'o ~ 200 cni^^ corresponding to 
A^pe ~ 0.6/cm m. 

The details of the luminosity measurement are de- 
scribed in Sec. IIIII 



F. Trigger systems 

The CDF trigger system (s^, [EHl was redesigned 
for Run II because of the changes in accelerator op- 
erating conditions. The upgraded trigger system re- 
duces the raw event rate in the detector (the nominal 
2.5 MHz beam crossing rate) to 75 Hz, the typical 
rate at which events can be recorded. 

The corresponding event rejection factor of roughly 
3 X lO'' is obtained using a three-level system where 
each level is designed to provide sufficient rejection 
to allow for processing with minimal deadtime at the 
subsequent level. The first level of the trigger system 
(Level 1) utilizes custom hardware to select events 
based on information in the calorimeters, tracking 
chambers, and muon detectors. Three parallel, syn- 
chronous hardware processing streams are used to 
create the trigger primitive data required to make 
the Level 1 decision. All detector data are fed into 
6 /iS pipelines to allow for processing time required 
at Level 1. The global Level 1 decision must be made 
and returned to the front-end detector hardware be- 
fore the corresponding collision data reach the end 
of the pipeline. Trigger decisions are made at the 
2.5 MHz crossing rate, providing dead-time free op- 
eration. 

One set of Level 1 hardware is used to find 
calorimeter objects (electrons and jets) and calcu- 
late the missing transverse energy and total trans- 
verse energy seen by the calorimeter in each event. 
At Level 1, electron and jet candidates are defined as 
single-tower energy deposits above some threshold in 
the electromagnetic or electromagnetic plus hadronic 
sections of trigger towers, respectively. Calorime- 
ter energy quantities are calculated by summing the 



transverse components of all single tower deposits as- 
suming a collision vertex of z = 0. A second set of 
hardware is utilized to select muon candidates from 
observed hits in the muon detector wire chamber and 
scintillator systems. A loose px threshold is applied 
based on differences in signal arrival times on pairs 
of projective wires in the CMU and CMX chambers. 
CMP primitives obtained from a simple pattern find- 
ing algorithm using observed hits on the four drift 
cell layers are matched to high px CMU candidates, 
and CSX hits within a certain time window consis- 
tent with collision-produced particles are matched to 
CMX candidates. 

An important element of the Run II CDF trigger 
upgrade is the third set of hardware which identifies 
COT track candidates within the tight Level 1 timing 
constraints. The eXtremely Fast Tracker (XFT) [5i| 
hardware examines hits on each axial superlayer of 
the COT and combines them into track segments. 
The found segments on the different layers are then 
linked to form tracks. The triggers used to collect 
the data samples utilized in these measurements are 
based on XFT tracks with reconstructed segments on 
all four COT axial superlayers. As discussed in more 
detail in Sec. IIVI this requirement has a small effect 
on the geometrical acceptance for lepton track candi- 
dates in our samples. The hit requirement for XFT 
track segments was changed from hits on 10/12 lay- 
ers to hits on 11/12 layers during the data collection 
period for the samples used in these measurements. 
This change led to a few percent drop in the trig- 
ger efficiency for high pT tracks but provided a sub- 
stantial increase in overall Level 1 event rejection. 
The XFT hardware reports tracks in 1.25° bins in 
(f). If more than one track is reconstructed within a 
given (f) bin, the track with the highest pT is used. 
The XFT feeds its lists of found tracks to another 
piece of hardware known as the track extrapolation 
unit (XTRP). The XTRP determines the number of 
tracks above certain pT thresholds and makes this 
information available for the global Level 1 decision. 
It also extrapolates each track based on its recon- 
structed pt into the calorimeter and muon detectors 
to determine into which slices of each system the 
track points based on the potential effects of mul- 
tiple scattering. This information is passed to the 
calorimeter and muon parts of the Level 1 trigger 
hardware in two sets of 2.5° (j) bins corresponding to 
groups of tracks above two programmable pt thresh- 
olds. Using this information, tracks are then matched 
to electron and muon primitives identified in those 
pieces of the Level 1 hardware to produce the final 
lists of electron and muon objects. 

The final Level 1 trigger decision is made based on 
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the number of physics objects (electrons, muons, jets, 
and tracks) found by the hardware and the calculated 
global calorimeter energy quantities. The maximum 
Level 1 event accept rate is roughly 20 kHz corre- 
sponding to an available Level 2 processing time of 50 
lis per event. Events accepted at Level 1 are stored in 
one of four buffers in the front-end readout hardware. 
Multiple event buffers allow for additional Level 1 
triggers to be accepted during the Level 2 processing 
of a previously accepted event. The Level 2 trigger 
system utilizes a combination of dedicated hardware 
and modified commercial processors to select events. 
There are two main pieces of dedicated Level 2 hard- 
ware. The first is the cluster finder hardware which 
merges the observed energy deposits in neighboring 
calorimeter towers to form clusters, and the second 
is the silicon vertex tracking hardware (SVT) [H^l 
which uses silicon detector hit information to search 
for tracks with displaced vertices. These systems are 
asynchronous in that processing time is dependent 
on the amount of input data associated with a given 
event. The output of these systems is passed to the 
global Level 2 processor along with the input data 
utilized in the Level 1 decision and additional hit in- 
formation from the CES to aid in low Et electron 
selection. The data are fed into the Level 2 processor 
board and simple selection algorithms, optimized for 
speed, are run to determine which events are passed 
to Level 3. The processor board has been designed to 
simultaneously read in one event while processing an- 
other which streamlines operation and helps to keep 
data acquisition deadtime at a minimum [55|. 



Events selected by the Level 2 trigger hardware 
are read out of the front-end detector buffers into 
the Level 3 processor farm. The current maximum 
Level 2 accept rate for events into Level 3 is roughly 
300 Hz. Level 3 processors run a speed-optimized 
version of the offline reconstruction code and impose 
loose sets of selection cuts on the reconstructed ob- 
jects to select the final 75 Hz of events which are 
recorded for further processing. The Level 3 proces- 
sor farm is made up of roughly 300 commercial dual 
processor computers running Linux to allow for one 
second of processing time for each event. The soft- 
ware algorithms run at Level 3 take advantage of the 
full detector information and improved resolution un- 
available at the lower trigger levels. The Level 3 algo- 
rithms are based on full three-dimensional track re- 
construction (including silicon hit information) which 
allows for tighter track matching with electromag- 
netic calorimeter clusters and reconstructed stubs in 
the muon detector for improved lepton identification. 



III. DATA SAMPLES AND LUMINOSITY 

The W ^ Iv and Z ^ li candidate event samples 
used to make the measurements reported here are se- 
lected from datasets collected using high Et lepton 
trigger requirements. Additional data samples used 
in the evaluation of efficiencies and backgrounds are 
discussed in further detail in the corresponding subse- 
quent sections. Here, we present the trigger require- 
ments for events contained within the datasets from 
which our candidate samples are selected. We also 
briefly describe data processing, the event quality cri- 
teria applied to our data samples, and the measure- 
ment of the integrated luminosities corresponding to 
our datasets. 



A. Trigger requirements 

The datasets used to select our candidate events 
are composed of events collected with well-defined 
trigger requirements at each of the three levels within 
the CDF trigger architecture (see Sec. HIl). The spe- 
cific trigger requirements associated with the datasets 
used to make our measurements are summarized 
here. The measured efficiencies of these trigger re- 
quirements are presented in Sec. IVfl 



1. Central electron trigger 

The trigger requirements for the dataset used to 
select W ^ ev and Z ^ ee candidate events are 
described here. Both candidate samples are selected 
from central, high Et electron triggered events, cor- 
responding to the region | ?7det | < 1.0. 

At Level 1, energies in physical calorimeter tow- 
ers of 0.1x15° in ridet-(f> space are first summed into 
0.2x15° trigger towers. At least one trigger tower is 
required to have Et > 8 GeV and the ratio of the 
hadronic to electromagnetic energies in that tower, 
Ehad/Ecm, must bc less than 0.125 (for measured 
Et < 14 GeV). In addition, at least one COT track 
with pt > 8 GeV/c pointing in the direction of the 
tower must be found by the XFT hardware. 

A clustering algorithm is run at Level 2 to 
combine associated energy deposits in neighboring 
calorimeter towers. Adjacent "shoulder" towers with 
Et > 7.5 GeV are added to the central "seed" tower 
found at Level 1. The total Et of the cluster is re- 
quired to be above 16 GeV and the E^i^d/Ecm ratio of 
the cluster is required to be less than 0.125. The pres- 
ence of an XFT track with pT > 8 GeV/c matched 
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to the seed tower of the central cluster is also recon- 
firmed. Finally, in Level 3 an electromagnetic clus- 
ter with Et > 18 GeV and £^had/-Eom < 0.125 must 
be found by the offline reconstruction algorithm. A 
track pointing at the cluster with pT > 9 GeV/c 
must also be found by the full three-dimensional COT 
track reconstruction algorithm run in the Level 3 pro- 
cessors. 

At each level of the trigger, the rate of accepted 
events is significantly reduced. At typical luminosi- 
ties (~ 2.5 X 10'^^cm~^s~^), the accepted rate of 
events for the above trigger requirements are 25 Hz, 
3 Hz, and 1 Hz for Levels 1, 2, and 3, respectively. 



2. Central muon triggers 

The dataset used to select our W ^ fiiy and 
Z —f fj,fi candidate samples is made of events col- 
lected using two analogous sets of trigger require- 
ments. In the most central region of the detector 
(I Vdct I < 0.6), trigger requirements are designed to 
select high px muon candidates which deposit hits 
in both the CMU and CMP wire chambers. An in- 
dependent but similar set of requirements is used to 
collect high px candidates in the extended central re- 
gion (0.6 < I 77dct I < 1-0) which produce hits in CMX 
wire chambers. 

The specific trigger requirements for the central re- 
gion at Level 1 are matched hits in one or more CMU 
projective wire pairs with arrival times within 124 ns 
of each other, a pattern of CMP hits on three of four 
layers consistent in (p with the observed CMU hits, 
and a matching COT track found by the XFT with 
Pt > 4 GeV/c. For the early part of the run period 
corresponding to our datasets we make no additional 
requirements at Level 2, but for the later portion we 
require at least one COT track with px > 8 GeV/c in 
the list of Level 1 XFT tracks passed to the Level 2 
processor boards. Because no muon trigger informa- 
tion was available at Level 2 during this run period, 
the higher px track was not required to match the 
CMU or CMP hits associated with the Level 1 trigger. 
Finally for Level 3, a reconstructed three-dimensional 
COT track with px > 18 GeV/c matched to recon- 
structed stubs in both the CMU and CMP chambers 
is required based on the offline reconstruction algo- 
rithms for muons. 

The analogous trigger requirements for the ex- 
tended central region at Level 1 are matched hits in 
one or more CMX projective wire pairs with arrival 
times within 124 ns of each other and a matching 
COT track found by the XFT with px > 8 GeV/c. 
For the latter part of our data collection period, a 



matching hit in the CSX scintillator counters consis- 
tent in time with a beam-produced particle is also 
required to help reduce the trigger rate from non- 
collision backgrounds. No additional requirements 
are made at Level 2. In Level 3, a reconstructed 
three-dimensional COT track with px > 18 GeV/c 
matched to a reconstructed stub in the CMX cham- 
bers is required based on the offline reconstruction 
algorithms for muons. 

At typical luminosities (~ 2.5 x 10'^^cm~^s~^), 
the accepted rate of events for the central trigger re- 
quirements are 30 Hz, 4 Hz, and 0.15 Hz for Lev- 
els 1, 2, and 3, respectively. For the extended central 
muon trigger requirements, the corresponding rates 
are 2 Hz, 2 Hz, and 0.1 Hz. 



B. Luminosity Measurement 

The total integrated luminosity (L) is derived from 
the rate of the inelastic pp events measured with 
CLC, Rpp, the CLC acceptance, cclCi a-nd the in- 
elastic pp cross section at 1.96 TeV, din, according to 
the expression 
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(7) 



The CLC acceptance is measured from tuned simu- 
lation and compared against the value obtained from 
a second method that relies on both data and simu- 
lation through the formula 
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where iVcLC+Piug is the number of inelastic events 
tagged by the CLC and plug calorimeter. New is a 
subset of those which contain an east-west hit co- 
incidence and pass the online selection criteria, and 
-^inelastic IS the total number of inelastic collisions. 
The fraction A^CLC+Piug/Mneiastic is extracted from 
simulation while the ratio iVEw/A^CLC+Piug is mea- 
sured from data. The acceptance calculated using 
this procedure is ecLC — 60.2 ± 2.6 % which is in 
good agreement with the value obtained directly from 
simulation. 

The value (Tin — 60.7 ± 2.4 mb is obtained by ex- 
trapolating the combined result for the inelastic pp 
cross section at ^/s = 1.8 TeV based on CDF and 
E811 measurements (59.3 ± 2.3 mb) [Hi to 1.96 TeV. 
Using these numbers, and restricting ourselves to 
runs with a good detector status, the total luminosity 
of our datasets is estimated to be 72.0 ± 4.3 pb""'^. 
The 6 % quoted uncertainty is dominated by the un- 
certainty in the absolute normalization of the CLC 
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TABLE II: Systematic uncertainties in the luminosity 
calculation based on the CLC measurement and the com- 
bined value of the CDF and E811 inelastic cross sec- 
tion measurements at y/s — 1.80 TeV extrapolated to 
y/s — 1.96 TeV. The total uncertainty in the CLC mea- 
surement is dominated by the uncertainty in the CLC ac- 
ceptance. The detector instability and calibration uncer- 
tainties are components of the overall CLC measurement 
uncertainty and therefore not included in the calculation 
of the total uncertainty. 



Effect 


Uncertainty Estimate 


Inelastic Cross Section 


4.0% 


CLC Measurement 


4.4% 


Detector Instability 


< 2.0 % 


Detector Calibration 


< 1.5% 


Total Uncertainty 


~ 6.0 % 



acceptance for a single pp inelastic collision . The 
complete list of systematic uncertainties, including 
uncertainties from the inelastic cross section and lu- 
minosity detector, is given in Table [TTl 



IV. EVENT SELECTION 

We search for W bosons decaying into a highly 
energetic charged lepton {£ = e,/z) and a neutrino, 
which is identified via large ]^rp in the detector. The 
Z £i {£ = e, iJ.) events are selected based on the two 
energetic, isolated leptons originating from a single 
event vertex. The two leptons produced in the decays 
are oppositely charged, and the charge information 
for leptons is included as part of the selection criteria 
when available. The reconstructed dilepton invariant 
mass is also required to lie within a mass window 
consistent with the measured Z boson mass. 

The complete set of selection criteria used to iden- 
tify W ^ tv and Z ^ II events are described here. 
As the selection of W and Z bosons proceeds from 
lepton identification, we also describe in some detail 
the variables used to select good lepton candidates. 



A. Track Selection 

The track quality requirements are common to 
electron and muon selection. As the silicon track- 
ing information is not vital to our measurements, 
we remove all silicon hits from the tracks and re- 
fit them, including the position of the beamline in 
the transverse direction as an additional constraint 
in the fit. The beamline position is measured in- 



dependently for each run period contributing to our 
datasets using the reconstructed COT track data con- 
tained within events from that period. The removal 
of silicon hits from tracks makes our measurements 
insensitive to the time-dependent efficiencies of the 
individual pieces of the silicon detector and allows 
us to include data from run periods when the silicon 
detector was not operational. The resulting beam- 
constrained COT tracks are used in the subsequent 
analysis work presented here. All of the kinematic 
track parameters used in these analyses with the one 
exception of the r — track impact parameter vari- 
able, do, used in muon selection, are based on these 
beam-constrained COT tracks. 

The reconstructed tracks obtained using the 
method described above have small residual curva- 
ture biases primarily due to COT misalignments that 
are not currently corrected for in our offline track- 
ing algorithms. We correct our track pT measure- 
ment for misalignment effects based on the observed 
0-dependence of the electron candidate Ejp distri- 
bution (see Sec. lIVBp . The form of the correction 
is 

= — - 0.00037 - 0.00110 • -smU + 0.28), (9) 
Vt Pt 

where and px are the transverse momenta in 
GeV/ c of the corrected and uncorrected track, respec- 
tively, Q is the charge of the track, and 4> is given in 
radians. 

We apply additional selection criteria on our re- 
constructed tracks to ensure that only high-quality 
tracks are assigned to lepton candidates. Each track 
is required to pass a set of minimum hit criteria. The 
reconstructed tracks are required to have a minimum 
of seven out of twelve possible hits on at least three 
of four axial and stereo superlayers within the COT. 
The minimum hit criteria for reconstructed tracks is 
less restrictive than that used to select Level 1 trig- 
ger track candidates (see Sec. IIIF[) to ensure high 
selection efficiencies for both triggerable and non- 
triggerable track candidates in our events. In ad- 
dition, to restrict ourselves to a region of high track 
reconstruction efhciency, we require the z coordinate 
of the lepton track intersection with the beam axis in 
the r ~ z plane, zq, be within 60 cm of the center of 
the detector. 

To help reduce real muon backgrounds from cosmic 
rays and ir/K decays, we impose additional quality 
requirements on muon track candidates. For muon 
track candidates only, we incorporate silicon hit in- 
formation in track reconstruction when available to 
calculate a more precise value for the r — impact 
parameter of the track, do- Cosmic ray muons and 



16 



muons produced m-K/K decays are less likely to point 
back to the event vertex and therefore will typically 
have larger measured impact parameters. We apply 
different cuts on the do of muon track candidates de- 
pending on whether or not the tracks contain any 
silicon hit information; |(io| < 0.2 cm for tracks with 
no silicon hits and jdol < 0.02 cm for tracks with sili- 
con hits. We also make a requirement on the quality 
of the final COT beam-constrained track fit for muon 
candidates. The track fit for muon backgrounds not 
originating from the event vertex will typically be 
worse when the additional constraint of the beamline 
position is included. For muon track candidates we 
require that x^/"-df < 2.0 where rtdf is the number 
of degrees of freedom in the track fit (the number of 
hits on the fitted track minus the five free parameters 
of the fit). 

We additionally restrict muon track candidates in 
to ensure that the tracks lie in a fiducial region of high 
trigger and reconstruction efficiency well-modeled by 
our detector simulation. We require that each muon 
track passes through all eight COT superlayers by 
making a minimum requirement on the exit radius 
of the track at the endplates of the COT tracking 
chamber. The exit radius is defined as 
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FIG. 5: The COT exit radius for CMX muons m Z ^ ^i^i 
candidate events. The points are the data and the his- 
togram is simulation. The selected CMX muons from 
data events are required to satisfy the high pr muon trig- 
ger criteria, but no trigger requirement is made on the 
muons selected from simulation. The two histograms are 
normalized to have the same number of events over the 
region 150 cm < pcoT < 280 cm. The arrow indicates 
the location of the muon track selection cut made on the 
pcoT variable. 



PCOT = {zcOT - zo) -iane, (10) 

where zcot is the distance of the COT endplates 
from the center of the detector (155 cm for tracks 
with > and —155 cm for those with 77 < 0). Here, 
T] and 6 are the previously defined pseudorapidity and 
polar angle of the track with respect to the directions 
of the colliding beams. A comparison of the pcoT 
distribution for CMX muons from Z fifi candidate 
events in data and Monte Carlo (MC) simulation (see 
Sec. El is shown in Fig. [51 The distributions do not 
match in the region pcOT < 140 cm due to a loss of 
data events in this region originating from the XFT 
track trigger requirements, at least ten (or eleven) 
hits out of a possible twelve for each of the four axial 
COT superlayers, which is not accounted for in the 
simulation. Based on this comparison, we require 
PcoT > 140 cm for muon track candidates. 

Track selection requirements are summarized in 
Table IIIII Distributions of the track quality variables 
used in the selection of all lepton tracks are shown in 
Fig. [HI and those used solely in the selection of muon 
tracks are shown in Fig[7l The distributions are con- 
structed from second, unbiased lepton legs in Z ^ ££ 
candidate data events. Based on these distributions, 
we expect the measured inefficiencies of our track se- 
lection criteria (see Sec. IVip to be on the order of a 
few percent. 



TABLE III: Summary of track selection requirements. 

Variable Cut 
All Tracks: 

# Axial COT Superlayers > 3 with > 7 hits 

# Stereo COT Superlayers > 3 with > 7 hits 
1^0 1 < 60 cm 
Muon Tracks: 

|do| < 0.2 cm (no silicon hits) 

I do I < 0.02 cm (silicon hits) 

x'Mdf < 2.0 

PCOT > 140 cm 



B. Electron Selection 

Electron candidates are reconstructed in either the 
central barrel or forward plug calorimeters. The clus- 
tering algorithms and selection criteria used to iden- 
tify electrons in the two sections are different, as we 
do not make use of tracking information in the for- 
ward detector region (| ?7det I > 1) where standalone 
track reconstruction is less reliable due to the smaller 
number of available tracking layers. Here, we discuss 
the specific identification criteria for both central and 
plug electrons. 
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FIG. 6: Distributions of the zo and number of axial and stereo COT superlayers contributing seven or more hits. These 
track quahty variables are from unbiased, second lepton legs oi Z U candidate events in data. The arrows indicate 
the locations of selection cuts applied on these variables. 



1. Central Electron Identification 

Electron objects are formed from energy clusters 
in neighboring towers of the calorimeter. An electron 
cluster is made from an electromagnetic seed tower 
and at most one additional tower that is adjacent to 
the seed tower in ijdet and within the same </> wedge. 
The seed tower must have Et > 2 GeV and a recon- 
structed COT track which extrapolates to that tower. 
The hadronic energy in the corresponding towers is 
required to be less than 0.125 times the electromag- 
netic energy of the cluster. 

Electron candidates for these measurements must 
lie within the well-instrumented regions of the 
calorimeter. The cluster position within the 
calorimeter is determined by the location of the asso- 
ciated CES shower. The CES shower must lie within 
21 cm of the tower center in the r — (j> view for the 
shower to be fully contained within the active region. 
We also exclude electrons reconstructed in the region 
where the two halves of the central calorimeter meet 



{\z\ < 9 cm) and the outer half of the most forward 
CEM towers (|z| > 230 cm) where there is substan- 
tial electron shower leakage into the hadronic part of 
the calorimeter. Finally, we exclude events in which 
the electron is reconstructed near the uninstrumented 
region surrounding the cryogenic connections to the 
solenoidal magnet (0.77 < ?7dct < 1.0, 75° < (/>< 90°, 
and \z\ > 193 cm). 

The selection requirements listed in Table IIVI 
are applied to electron candidates in the well- 
instrumented regions of the central calorimeter. We 
cut on the ratio of the hadronic to electromag- 
netic energies, i^hadZ-E-cm, for the candidate clusters. 
Electron showers are typically contained within the 
electromagnetic calorimeter, while hadron showers 
spread across both the hadronic and electromagnetic 
sections of the calorimeter. We require i^had/^^cm < 
0.055 + 0.00045 -E where E is the total energy of 
the cluster in GeV. The linear term in our selection 
criteria accounts for the increased shower leakage of 
higher-energy electrons into the hadronic calorimeter 
sections. 
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FIG. 7: Distributions of the do (with and without attached silicon hits), x^/^df, and pcoT- These track quality 
variables are for muons from unbiased, second muon legs of Z — > /xp candidate events in data. The arrows indicate the 
locations of selection cuts applied on these variables. 



We also cut on the ratio of the electromagnetic 
cluster transverse energy to the COT track transverse 
momentum, E /p. This ratio is nominally expected to 
be unity, but in cases where the electron radiates a 
photon in the material of the inner tracking volume, 
the measured momentum of the COT track can be 
less than the measured energy of the corresponding 
cluster in the calorimeter. In cases where the electron 
is highly energetic, the photon and electron will be 
nearly coUinear and are likely to end up in the same 
calorimeter tower. The measured COT track momen- 
tum will, however, correspond to the momentum of 
the electron after emitting the photon and thus be 
smaller than the original electron momentum. We 
require E/p < 2.0 which is efficient for the major- 
ity of electrons which emit a bremsstrahlung photon. 
Since this cut becomes unreliable for very large values 
of track pT, we do not apply it to electron clusters 
with Et > 100 GeV. 

The lateral shower profile variable, Lghr [13]) is 
used to compare the distribution of adjacent CEM 
tower energies in the cluster as a function of seed 



tower energy to shapes derived from electron test- 
beam data. We also perform a comparison of the 
CES lateral shower profile in the r—z view to the pro- 
file extracted from the electron test-beam data. For 
central electrons, we require Lghr < 0.2 and Xstrips < 
10.0. 

Since central electron candidates include a COT 
track, we can further reduce electron misidentifica- 
tion by cutting on track-shower matching variables. 
We define Q ■ Ax as the distance in the r — (j) plane be- 
tween the extrapolated beam-constrained COT track 
and the CES cluster multiplied by the charge of the 
track to account for asymmetric tails originating from 
bremsstrahlung radiation. The variable Az is the cor- 
responding distance in the r — z plane. We require 
— 3.0 cm < Q ■ Ax < 1.5 cm and |Az| < 3.0 cm. 

Distributions of central electron identification vari- 
ables are shown in Figs. [5] and [51 The plotted electron 
candidates are the unbiased, second electron legs in 
Z — > ee events in which both electrons are recon- 
structed within the central calorimeter and the first 
electron is found to satisfy the full set of identifica- 
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TABLE IV: Calorimeter variables and electron identifica- 
tion requirements. 



Variable 


Cut 


Central 


1 met 1 < 1.0 


Ehad/Ecm 


< 0.055 + 0.00045 --EfGeV] 


E/p (for Et < 100 GeV) 


< 2.0 


ishr 


< 0.2 




> -3.0 cm, < 1.5 cm 




< 3.0 cm 


2 

Xstrips 


< 10.0 


Plug 


1.2 < 1 r/dct 1 < 2.8 


-EhadZ-Eem 


< 0.05 


2 

Xpem 


< 10.0 



tion criteria. Based on these distributions, we expect 
a high efficiency for our central electron selection cri- 
teria (see Sec. IVip. 



2. Plug Electron Identification 

Electron candidate clusters in the plug calorimeter 
are made from a seed tower and neighboring towers 
within two towers in 77dct and from the seed tower. 
As for central electrons, the hadronic energy of the 
cluster is required to be less than 0.125 times the elec- 
tromagnetic energy. We also require plug electrons to 
be reconstructed in a well-instrumented region of the 
detector, defined as 1.2 < | T^dot I < 2.8. 

The additional selection criteria applied to plug 
electron candidates are summarized in Table IIVI 
Fewer variables are available for selecting plug elec- 
trons due to the lack of matching track information 
for candidates in the forward region of the detector. 
As in the case of central electrons, we cut on the 
ratio of hadronic to electromagnetic energies in the 
cluster, -Ehad/-Ecm, which is required to be less than 
0.05. We also compare the distribution of tower en- 
ergies in a 3x3 array around the seed tower to dis- 
tributions from electron test-beam data, forming the 
variable Xpem which we require to be less than 10.0. 
Distributions of the plug electron selection variables 
are shown in Fig. [TOl The plotted electron candi- 
dates are the unbiased, second plug electron legs in 
Z ^ ee events in which the first electron is recon- 
structed within the central calorimeter and found to 
satisfy a set of more restrictive cuts on the previously 
described central electron identification variables. 



C. Muon Selection 

Muon candidates used in these measurements must 
have reconstructed stubs in both the CMU and CMP 
chambers (CMUP muons) or a reconstructed stub 
in the CMX chambers. CMX chambers were offline 
for the first 16.5 pb^^ of integrated luminosity cor- 
responding to our datasets, and the reduced muon 
detector coverage during this period is taken into 
account in our measured acceptances for events in 
the muon candidate samples (see Sec. IV Cp . The 
muon candidate tracks are required to extrapolate 
to regions of the muon chambers with high single 
wire hit efficiencies to ensure that chamber-edge ef- 
fects do not contribute to inefficiencies in muon stub- 
reconstruction and stub-track matching (see Sec. lVip . 
We measure the location of an extrapolated muon 
track candidate with respect to the drift direction 
(local x) and wire axis (local z) of a given cham- 
ber. The extrapolation assumes that no multiple 
scattering takes place, and in some cases muons that 
leave hits in the muon detectors extrapolate to lo- 
cations outside of the chambers. In the CMP and 
CMX chambers, we require that the extrapolation is 
within the chamber volume in local x, and at least 
3 cm away from the edges of the chamber volume 
in local z. Studies of unbiased muons in Z — s- /i/i 
events show that these regions of chambers are max- 
imally efficient for hit-finding. No such requirement 
is needed for the CMU chambers. Some sections of 
the upgraded muon detectors were not yet fully com- 
missioned for the period of data-taking corresponding 
to our datasets, and we exclude all muon candidates 
with stubs in those sections. 

The selection criteria applied to muon candidates 
are summarized in Table|Vl We require that the mea- 
sured energy depositions in the electromagnetic and 
hadronic sections of the calorimeters along the muon 
candidate trajectory, Ecm and i?had, are consistent 
with those expected from a minimum-ionizing par- 
ticle. The positions of the reconstructed chamber 
stubs are required to be near the locations of the 
extrapolated tracks. The track-stub matching vari- 
able I AX I is the distance in the r — (f) plane between 
the extrapolated COT track and the CMU, CMP, or 
CMX stub. Fig. [Tl] shows the /S.X distributions for 
unbiased, CMU, CMP and CMX second muons in 
Z fifi events. 

Energetic cosmic ray muons traverse the detector 
at a significant rate, depositing hits in both muon 
chambers and the COT, and can in a small frac- 
tion of cases satisfy the requirements of the high px 
muon trigger paths and the offiine selection criteria. 
We remove cosmic ray events from our sample using 
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FIG. 8: Distributions of Lshr, E/p, -Ehad/Sem, and Ex°/Et (see Sec. IIVD|| central calorimeter electron selection 
variables from unbiased, second electron legs of Z ^ ee candidate events in data. The arrows indicate the locations of 
selection cuts applied on these variables. No arrow is shown on the iJiiad/iJcm distribution since the cut on this variable 
is dependent on the electron energy. 



the previously discussed track quality cuts for muon 
candidates and a cosmic ray tagging algorithm (see 
Sec. IVII E[) based on COT hit timing information. 



D. W ^ Selection 

W events are selected by first requiring a 

high-pT charged lepton in the central detectors, as de- 
scribed above. Electrons must have electromagnetic- 
cluster Et greater than 25 GeV and COT track px 
greater than 10 GeV/c. Muons must have COT track 
Pt greater than 20 GeV/c. The leptons from de- 
cays of W and Z bosons are often isolated from 



hadronic jets, in contrast to leptons originating from 
decays of heavy-flavor hadrons. We therefore re- 
quire that the calorimeter energy in a cone of radius 
AR = \J Arf + < 0.4 around the lepton ex- 
cluding the energy associated with the lepton, E^fP , 
be less than 10 % the energy of the lepton (Et for 
electrons and pt for muons). Fig. [5] shows the isola- 
tion distribution for unbiased central electrons from 
Z ^ ee decays. 

We also require evidence for a neutrino in W can- 
didate events in the form of an imbalance of the mea- 
sured momentum of the event since neutrinos do not 
interact with our detector. The initial state of the 
colliding partons has — 0, but unknown pz due 
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FIG. 9: Distributions of the zq, Xstripsi Q ' ^a;, and Az central calorimeter electron selection variables from unbiased, 
second electron legs of Z —> ee candidate events in data. The arrows indicate the locations of selection cuts applied on 
these variables. 



to the unknown value of initial parton momentum. 
Therefore, we identify the missing transverse energy, 
$Jj, , in the event with the px of the neutrino. In muon 
events, we correct the l^rp measured in the calorime- 
ter to account for the energy carried away by the 
muon, a minimum-ionizing particle. The muon mo- 
mentum is used in place of the calorimeter energy 
deposits observed along the path of the muon. For 
W ^ fxv candidate events we require If^j. > 20 GeV 
and tighten the requirement for W ^ ev events, 
> 25 GeV, to further reduce backgrounds from 
hadron jets. 

A background to W ^ tv is the Z ^ U channel, 
when one of the leptons falls into an uninstrumented 
region of the detector, creating false If^rp . This is a 



bigger problem in the muon channel since the cover- 
age of the muon detectors is in general less uniform 
than that of the calorimeter. Therefore, in the muon 
channel we reject events with additional minimum- 
ionizing tracks in the event with pT > 10 GeV/c, 
< 3 GeV (-Bern < 3 + 0.0140 • (p-100) GeV 
if p > 100 GeV/c) and E^^ < 6 GeV (^Bhad < & + 
0.0420 • (p-100) GeV if p > 100 GeV/c). Studies 
of simulated W iiv and Z — s- /i/i event samples 
(see Sec. |V| show that this additional rejection crite- 
ria removes 54.7% of the Z ^ fifi background while 
retaining 99.6% of the W ^ signal. Further dis- 
cussion of backgrounds to the W ^ iv channels is 
found in Section IVlII 
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FIG. 10: Distributions of the Xpemi (see Sec. lIVDil . and Shad/Scm plug calorimeter electron selection variables 
from unbiased, second electron legs of Z — > ee candidate events in data. The arrows indicate the locations of selection 
cuts applied on these variables. 



TABLE V: Calorimeter and muon chamber variables used 
in muon identification. The fiducial distance variables are 
defined as the extrapolated position of the muon track 
candidate with respect to the edges of a given muon cham- 
ber. The fiducial distance is negative if this position lies 
within the chamber and positive otherwise. CMUP muon 
candidates are those with reconstructed stubs in both the 
CMU and CMP detectors. CMX muon candidates have 
reconstructed stubs in the CMX detector. 



Variable 






Cut 


Minimum Ionizing Cuts: 






(GeV) 


£cm (p < 100 GeV/c) 






< 2 


Eem (p > 100 GeV/c) 


< 2 + (p- 


-100) 


■ 0.0115 


Shad (p < 100 GeV/c) 






< 6 


Shad (p > 100 GeV/c) 


< 6 + (p- 


-100) 


• 0.0280 


Muon Stub Cuts: 






(cm) 


lAXcMul (CMUP) 






< 3.0 


IAXcmpI (CMUP) 






< 5.0 


lAXcMxl (CMX) 






< 6.0 


CMP a;-fiducial distance (CMUP) 






< 0.0 


CMP z-fiducial distance (CMUP) 






< -3.0 


CMX a;-fiducial distance (CMX) 






< 0.0 


CMX z-fiducial distance (CMX) 






< -3.0 



E. 



££ Selection 



ses. We identify a second lepton in these events using 
less restrictive ( "loose" ) selection criteria to increase 
our Z ££ detection efficiency. The invariant mass 
of the two leptons is required to be between 66 and 
116 GeV/c2. 

After selecting the first electron, Z ee events are 
identified by the presence of another isolated electron 
in the central calorimeter with Et > 25 GeV passing 
less restrictive selection criteria or an isolated elec- 
tron in the plug calorimeter with Et > 20 GeV. The 
selection criteria for each type of electron are sum- 
marized in Table IVTl In the calculation of Et for the 
plug electron, the z-vertex is taken from the COT 
track of the central electron in the event. In the case 
of two central electrons, we also require they be op- 
positely charged, with both electron tracks passing 
the track quality criteria listed in Table IIIII 

After selecting the first muon, Z — > /i/i events 
are identified by the presence of another oppositely 
charged, isolated muon track with pT > 20 GeV orig- 
inating from a common vertex. The muon-stub cri- 
teria are dropped for the second leg to gain signal 
acceptance with very little increase in background; 
this second muon is merely a minimum-ionizing track. 
Table IVIII shows the complete set of selection crite- 
ria used to identify Z fifi events. Again, we re- 
quire both tracks to pass the track requirements of 
Table [ml 



We select events which contain an electron or muon 
that passes the same identification requirements as 
the lepton inW ^ £v candidate events. As described 
in Sec. IH systematic uncertainties are reduced by us- 
ing a common lepton selection in the W and Z analy- 



F. Event Selection Summary 

Using the selection criteria described here, we find 
a total of 37,584 W ^ ev candidate events. In the 
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FIG. 11: Distributions of the CMU, CMP, and CMX AX muon selection variables from unbiased, second muon legs 
of Z fifj. candidate events in data. The arrows indicate the locations of selection cuts applied on these variables. 



TABLE VI: Z ~> ee selection criteria. 

Variable "Tight" Central e "Loose" Central e Plug e 

Et > 25 GeV > 25 GeV > 20 GeV 

PT > 10 GeV/c > 10 GeV/c 

E!f° < 0.1 -£5}"=*" < 0.1 -^5,'"='" < 4 GeV 

-Ehad/Scm < 0.055 + 0.00045 • E < 0.055 + 0.00045 ■ E < 0.05 

E/p < 2.0 (or PT > 50 GeV/c) 

ishr < 0.2 

Q ■ Ax > —3.0 cm, < 1.5 cm 
lA^i < 3.0 cm 

Xstrips < 10.0 

Xpem < 10-0 



muon channel, we find 21,983 W boson candidates 
with CMUP muons and 9,739 with CMX muons for 
a grand total of 31,722 W ^ jiv candidates. In the 
Z boson decay channel, we find 1,730 events with 
two reconstructed electrons in the central calorime- 
ter and 2,512 events in which the second electron is 
reconstructed in the plug calorimeter giving a total 



of 4,242 Z ^ ee candidates. From our high pT muon 
dataset, we find 1,371 CMUP + track and 677 CMX 
+ track Z ^ip, candidates. There is an overlap of 
263 events between these two samples in which one 
candidate track is matched to stubs in the CMU and 
CMP muon chambers and the other is matched to a 
stub in the CMX chamber. Taking this overlap into 
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TABLE VII: Z fifj, selection criteria. 



Variable 






Cut 


Fiducial and Kinematic: 








\vl \ 


< 1.0 (CMUP+CMX) 


1 




< 1.2 (Track) 


Pt 






> 20 GeV/c 


m{2) 

Pt 






> 20 GeV/c 


Both Muon Legs: 








Ecm (p < 100 GeV/c) 






< 2 GeV 


Ecm (p > 100 GeV/c) 


< 2 + (p- 


-100) 


■ 0.0115 GeV 


£had (p < 100 GeV/c) 






< 6 GeV 


Shad (p > 100 GeV/c) 


< 6 + (p- 


-100) 


■ 0.0280 GeV 


r-iiso 






< 0.1 -pT 


First Muon Leg: 








1 AXcMul 




< 3.0 


cm (CMUP) 


1 AXcmpI 




< 5.0 


cm (CMUP) 


1 AXcMxl 




< 6.0 cm (CMX) 


CMP a;-fiducial distance 




< 0.0 


cm (CMUP) 


CMP 2-fiducial distance 


< 


: -3.0 


cm (CMUP) 


CMX x-fiducial distance 




< 0.0 cm (CMX) 


CMX z-fiducial distance 




< -3.0 cm (CMX) 



account, we obtain a total of 1,785 Z candidate 
events. 



V. SIGNAL ACCEPTANCE 
A. Introduction 

The acceptance terms in Eqs. [5] and |6] are defined 
as the fraction oi W —* iv or Z ^ i£ events produced 
in pp collisions at y/s — 1.96 TeV that satisfy the geo- 
metrical and kinematic requirements of our samples. 
Lepton reconstruction in our detector is limited by 
the finite fiducial coverage of the tracking, calorime- 
ter, and muon systems. Several kinematic require- 
ments are also made on candidate events to help re- 
duce backgrounds from non-W/Z processes. The re- 
constructed leptons in these events are required to 
pass minimum calorimeter cluster Et and/or track 
Pt criteria. In addition, a minimum requirement on 
the total measured missing Et is made on events in 
our W ^ candidate samples, and the invariant 
mass of Z ^ U candidate events is restricted to a 
finite range around the measured Z boson mass. 

The fraction of signal events that satisfy the ge- 
ometrical and kinematic criteria outlined above for 
each of our samples is determined using simulation. 
One geometrical cut on candidate events for which 



we measure the acceptance directly from data is the 
requirement that the primary event vertex for each 
event lies within 60.0 cm of the detector origin along 
the z-axis (parallel to the direction of the beams). 
Our simulation does include a realistic model of the 
beam interaction region, but we obtain a more ac- 
curate estimation of the selection efficiency for the 
event vertex requirement from studies of minimum 
bias events in the data as described in Sec. |Vll Since 
the geometrical and kinematic acceptance for candi- 
date events with a primary vertex outside our allowed 
region is significantly smaller, we remove the subset 
of simulated events with vertices outside this region 
from our acceptance calculations to avoid double- 
counting correlated inefficiencies. 

There is one additional complication involved in 
determining the kinematic and geometrical accep- 
tances for our Z ^ U candidate samples. Be- 
cause we make our ^* /Z i£ production cross sec- 
tion measurements in a specific invariant mass range, 
66 GeV/c^ to 116 GeV/c^, we need to account for 
events outside this mass range that are reconstructed 
in the detector to sit within this range due to the 
effects of detector resolution and final state radia- 
tions. To include events of this type in our accep- 
tance calculations, we use simulated ^* / Z ££ event 
samples generated over a wider invariant mass range 
{Mu > 30 GeV/c^). In order for an event to con- 
tribute to the denominator of our acceptance calcu- 
lations, we require that the invariant mass of the lep- 
ton pair at the generator level prior to application 
of final state radiative effects lies within the range of 
our measurement (66 GeV/c^ < M« < 116 GeV/c^). 
The generator-level invariant mass requirement is not 
made on events contributing to the numerator of our 
acceptance calculations, however, so that ^* /Z ^ ££ 
events generated outside the invariant mass range of 
our measurement which have reconstructed masses 
within this range are properly accounted for in the 
acceptance calculations. 



B. Event and Detector Simulation 

The simulated events used to estimate the accep- 
tance of our samples were generated with pythia 
6.203 [Hi. The default set of Parton Distribution 
Functions (PDFs) used in the generation of these 
samples is CTEQ5L [g^]- PYTHIA generates the hard, 
leading-order (LO) QCD interaction, q + q ^ 1* /Z 
(or q + q' ^ W), simulates initial-state QGD radia- 
tion via its parton-shower algorithms, and generates 
the decay, "f*/Z ££ (or W £v). No restrictions 
were placed at the generator level on the pt of the 
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final-state leptons or on their pseudorapidity. Both 
initial- and final-state radiation were turned on in the 
event simulation. In order to model the data accu- 
rately, the beam energy was set to 980 GeV, and the 
beam parameters were adjusted to provide the best 
possible match with data. The profile of the beam 
interaction region in z was matched to data by set- 
ting the mean of the vertex distribution to 3.0 cm in 
the direction along the beams and the correspond- 
ing Gaussian spread to 25.0 cm. The offset of the 
beam from the nominal center of the detector in the 
J — (/) plane is also taken into account. In the simula- 
tion, the position of the beams at z = is ofi'set by 
-0.064 cm in x and -1-0.310 cm in y to provide a rough 
match with the measured offsets in data. The loca- 
tion of a given vertex within the r — (f> plane is also 
observed to depend on its location along the z-axis 
due to the non-zero angle between the beams and the 
central axis of the detector. Slopes of -0.00021 and 
0.00031 are assigned in the simulation to the direc- 
tion of the beams relative to the y — z and x — z 
detector planes to model this effect. 

The intermediate vector boson px distribution in 
the simulation is tuned to match the CDF Run I 
measurement of the da/dpT spectrum for electron 
pairs in the invariant mass range between 66 GeV/ 
and 116 GeV/c^ [6H. The tuning is done using 
PYTHIA's nonperturbative "Kt smearing" parame- 
ters, PARP(91) and parp(93), and shower evolu- 
tion parameters, parp(62) and parp(64). The 
PARp(91) parameter affects the location of the peak 
in the da/dpr distribution in the vicinity of 3 GcV/c, 
and the parp(62) and parp(64) parameters affect 
the shape of the distribution in the region between 
7 GeV/c and 25 GeV/c. A comparison between the 
"tuned" 'y* / Z px distribution from simulation and 
the measured Run I spectrum is shown in Fig. 1121 
We assume that the optimized pythia tuning pa- 
rameters obtained from data collected at the Run I 
center of mass energy {y/s = 1.80 TeV) remain valid 
at the increased Run II center of mass energy (-^5 = 
1.96 TeV). The underlying event model in pythia is 
also tuned based on observed particle distributions in 
minimum bias events [63 |. 

The shape of the boson rapidity distribution is 
strongly dependent on the choice of PDFs. The shape 
of the da/dy distribution for j* /Z — > ee pairs in the 
mass region, 66 GeV/c^ < Mee < 116 GeV/c^, was 
measured by CDF in Run l[2^. The good agreement 
observed between the measured shape of da/dy with 
that obtained from simulation using CTEQ5L PDFs 
motivates the selection of this PDF set for our event 
generation. A comparison between the shape of the 
Run I measured da/dy distribution and the shape 
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FIG. 12: Tuned pythia 6.21 da/dpr in pb per GeV/c 
(on average) of 'y'/Z ee pairs in the mass region 
66 GeV/c^ < Mee < 116GeV/c2 (histogram) versus the 
measurement made by CDF in Run I (points). 

of the same distribution from pythia 6.21 simulated 
event samples generated with CTEQ5L is shown in 
Fig. [131 

A detector simulation based on GEANt3 [H, [g^l is 
used to model the behavior of the CDF detector. The 
GFLASH [6^ package is used to decrease the simula- 
tion time of particle showers within the calorimeter. 

C. Signal Acceptances from pythia 

Additional tuning is performed after detector sim- 
ulation to improve modeling of the data further. A 
detailed description of the techniques used to deter- 
mine the post-simulation tunings described here and 
the associated acceptance uncertainties is provided in 
sections EH] and El 

The tuned, simulated event samples are used to 
determine the acceptances of each W and Z event 
sample. As discussed in Sec. IV Al events with a pri- 
mary event vertex outside our allowed region (Izvtxl < 
60 cm) are removed from both the numerator and de- 
nominator of our acceptances. The W ^ tv accep- 
tance calculations are outlined in Tables IVIIII and HXl 
for the electron and muon candidate samples. The ge- 
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FIG. 13: Tuned pythia 6.21 da/dy in pb per 0.1 of 
7*/Z — > ee pairs in the mass region 66 GeV/c^ < 
Mee < 116 GeV/c^ (histogram) versus the measurement 
made by CDF in Run I (points). 



ometric and kinematic requirements listed in each ta- 
ble define the acceptances for the corresponding sam- 
ples. The number of simulated events which satisfy 
each of the successive, cumulative criteria is shown in 
the tables along with the resulting net acceptances 
based on the total number of events with primary 
vertices inside our allowed region. The W ^ 
events which contain reconstructed muons with stubs 
in the CMX region of the muon detector are assigned 
a weight of 55.5/72.0 in the numerator of the ac- 
ceptance calculation to account for the fact that the 
CMX detector was offline during the first 16.5 pb~^ 
of integrated luminosity that define our samples. The 
largest uncertainties attached to the individual lumi- 
nosity measurements (see Sec. IIII Bp cancel in our 
weighting ratio for CMX events and the residual un- 
certainty on this ratio has a negligible effect on the 
overall acceptance uncertainty. 

The Z ^ U acceptance calculations are outlined 
in Tables |X] and IXII for the corresponding electron 
and muon candidate samples. As previously stated 
in Sec. IV Ai the acceptances that we define for these 
samples are for 7*/^ — *■ it in the invariant mass range 
66 GeV/c2 < Mu < 116 GeV/c^. The simulated event 
samples used to estimate the Z boson acceptances are 



generated with a looser invariant mass requirement, 
Mu > 30 GeV/c^. Generated events with an invari- 
ant mass outside our allowed range do not contribute 
to the denominator of these acceptance calculations 
but can contribute to the numerator if the final recon- 
structed invariant mass turns out to lie within our al- 
lowed region due to radiative and/or detector resolu- 
tion effects. In order for an event to contribute to the 
denominator of the Z boson acceptance calculations, 
we require that the invariant mass of the dilepton pair 
at the generator level before application of any final 
state radiative effects lies within the correct invariant 
mass range, 66 GeY/cP < M«(Gen) < 116 GeV/c^. 
As in the case of the W ^ acceptance calculation, 
events in the numerator of the Z —^ fifi acceptance 
calculation must be weighted to account for the fact 
that the CMX portion of the muon detector was of- 
fline for the first subset of integrated luminosity that 
defines our samples. In order to account for this ef- 
fect, a weight of (55.5/72.0) is applied to events con- 
tributing to the numerator of the Z —^ fi/j. acceptance 
calculation which contain a CMX muon candidate 
satisfying the three muon geometric and kinematic 
requirements listed in Table IXII and no CMUP muon 
candidates satisfying these same three requirements. 



D. Improved Acceptance Calculations 

The tuned PYTHIA simulated event samples are de- 
signed to provide the best possible model for our 
W and Z boson candidate samples. However, the 
actual boson production cross section calculation 
made by pythia is done only at leading order (LO); 
see Fig. [T] The complex topologies of higher-order 
contributions are modeled using a backward shower 
evolution algorithm which includes initial-state ra- 
diative effects and a separate, post-generation al- 
gorithm for including final-state radiation. A bet- 
ter description of boson production can be obtained 
from recently developed NNLO theoretical calcula- 
tions of the double-differential production cross sec- 
tions, (Pa/dydM , as a function of boson rapidity (y) 
and mass (M), for both and Drell-Yan produc- 
tion [l3|. The calculations are based on the MRST 
2001 NNLL PDF set ^ and electroweak param- 
eters taken from The single differential cross 
sections, da/dy, are obtained by integrating over the 
mass range, 66 GeV/c^ < Mu < 116 GeV/c^ for Drell- 
Yan production and 40 GeV/c^ < Mi^, < 240 GeV/c^ 
for production. 

We use these NNLO theoretical calculations of 
da/dy for Drell-Yan and production to ob- 

tain improved acceptance estimates for our candidate 
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TABLE VIII: W ^ eu selection acceptance from pythia Monte Carlo simulation. Statistical uncertainties are shown. 



Selection Criteria 



Number of Events Net Acceptance 



Total Events 1933957 

l^vtxl < 60 cm 1870156 

Central EM Cluster 927231 0.4958 ± 0.0004 

Calorimeter Fiducial Cuts 731049 0.3909 ± 0.0004 

Electron Track pr > 10 GeV/c 647691 0.3463 ± 0.0003 

EM Cluster Et > 25 GeV 488532 0.2612 ± 0.0003 

Event > 25 GeV 447836 0.2395 ± 0.0003 



TABLE IX: W ~* selection acceptance from pythia Monte Carlo simulation. Statistical uncertainties are shown. 



Selection Criteria 



Number of Events Net Acceptance 



Total Events 2017347 

|2vtx| < 60 cm 1951450 

CMUP or CMX Muon 545221 0.2794 ± 0.0003 

Muon Chamber Fiducial Cuts 523566 0.2683 ± 0.0003 

Muon Track pr > 20 GeV/c 435373 0.2231 ± 0.0003 

Muon Track Fiducial Cuts 411390 0.2108 ± 0.0003 

Event > 20 GeV 383787 0.1967 ± 0.0003 



samples. First, the tuned PYTHIA event simulation is 
used to create acceptance functions for each candi- 
date sample as a function of boson rapidity, A{y). 
These functions provide the acceptance in each bo- 
son rapidity bin based on our modehng of the CDF 
detector contained in the event simulation. Fig. [14] 
shows the A{y) acceptance functions for each of our 
four candidate samples. The Z — > ee sample has 
larger acceptance at higher boson rapidity due to the 
plug calorimeter modules which provide additional 
coverage in the forward part of the detector for the 
second electron in these events. Based on these dis- 
tributions, the acceptance of each sample. A, is then 
calculated as 



_^J^.A{y)-dy 



(11) 



The acceptance values obtained with this approach 
are shown in Table IXIII and compared with values 
obtained directly from the pythia simulated event 
samples. The results all agree within 0.4 % indicat- 
ing that the shapes of the NNLO da/dy distributions 
are very similar to those computed with the pythia 
simulation. The acceptance values obtained using the 
NNLO theoretical differential cross section calcula- 
tions are used for our measurements. 



E. Uncertainties in NNLO Calculation 

Uncertainties in the NNLO calculations of the dif- 
ferential boson production cross sections lead to un- 
certainty on our calculated acceptance values. The 
theoretical calculations require a large number of in- 
put parameters taken from world average experimen- 
tal results that have their own associated uncertain- 
ties. The renormalization scale used in the calcula- 
tions is another source of uncertainty. The default 
renormalization scales used in the calculations are 
Mz for Drell-Yan production and Mw for pro- 
duction. To study the effect of this scale on our cen- 
tral acceptance values, we recalculate the da/dy pro- 
duction cross sections using renormalization scales 
twice and one-half of the default values. For both 
cases, we find the net change in our calculated accep- 
tances to be less than 0.1 % which has a negligible 
effect on our overall acceptance uncertainty. 

We perform a computational consistency check on 
the NLO component of the NNLO da/dy calcula- 
tion IITI w ith a different MS NLO computation of 
dcr/dy]ll|,[il,[il,M.]. We find that the resulting ac- 
ceptance values differ by no more than 0.1 %. Based 
on this agreement between the two calculations, we 
assign no additional uncertainty to our acceptance 
values based on the calculation itself. However, our 
default calculation is still susceptible to uncertainties 
from higher order effects beyond NNLO. To place a 



28 



TABLE X: Z ^ ee selection acceptance from pythia Monte Carlo simulation. Statistical uncertainties are shown. 



Selection Criteria Number of Events Net Acceptance 



Total Events 


507500 








Zvtx < 60 cm 


490756 








66 GeV/c^ < Mee(Gen) < 116 GeV/c^ 


376523 








Central EM Cluster 


363994 


0.9667 


± 


0.0003 


Calorimeter Fiducial Cuts 


299530 


0.7955 


± 


0.0007 


Electron Track pr > 10 GeV/c 


252881 


0.6716 


± 


0.0008 


EM Cluster Et > 25 GeV 


186318 


0.4948 


± 


0.0008 


Second EM cluster (Central or Plug) 


176417 


0.4685 


± 


0.0008 


Second Cluster Calorimeter Fiducial Cuts 


146150 


0.3882 


± 


0.0008 


Second Electron Track pr > 10 GeV/c (Central) 


138830 


0.3687 


± 


0.0008 


Second EM Cluster Et > 25 GeV (Central), 20 GeV (Plug) 


125074 


0.3322 


± 


0.0008 


Second EM Cluster Shad/^^em < 0.125 (Plug) 


124881 


0.3317 


± 


0.0008 


66 GeV/c^ < Me,(Rec) < 116 GeV/c^ 


120575 


0.3202 


± 


0.0008 


Opposite Charge (Central-Central) 


119925 


0.3185 


± 


0.0008 



TABLE XI: Z n/j, selection acceptance from pythia Monte Carlo simulation. Statistical uncertainties are shown. 



Selection Criteria Number of Events Net Acceptance 



Total Events 


507500 




Zvtx < 60 cm 


490755 




66 GeV/c^ < M^^(Gen) < 116 GeV/c^ 


375981 




CMUP or CMX Muon 


217041 


0.5773 ± 0.0008 


Muon Chamber Fiducial Cuts 


209693 


0.5577 ± 0.0008 


Muon Track Fiducial Cuts 


199940 


0.5318 ± 0.0008 


Muon Track pr > 20 GeV/c 


157244 


0.4182 ± 0.0008 


Second Track with pr > 10 GeV/c 


91048 


0.2422 ± 0.0007 


Second Track Fiducial Cuts 


62663 


0.1667 ± 0.0006 


Second Track pr > 20 GeV/c 


56459 


0.1502 ± 0.0006 


66 GeV/c^ < Mp^(Rec) < 116 GeV/c' 


52160 


0.1387 ± 0.0006 



TABLE XIL Central acceptance values for our candidate 
samples based on da/dy distributions obtained from both 
NNLO and pythia simulation. 

Acceptance NNLO Calc. pythia Difference (Vo) 



Aw-^t,,. 0.1970 0.1967 +0.15 

Aw^^„ 0.2397 0.2395 +0.08 

0.1392 0.1387 +0.36 

Az^ee 0.3182 0.3185 -0.09 

Az^^,^,/Aw-,^.. 0.7066 0.7054 +0.17 

Az^ee/Aw^,^ 1.3272 1.3299 -0.20 



conservative limit on the magnitude of higher-order 
uncertainties, we compare acceptance values based on 
the NLO and NNLO versions of our default da/dy 



production cross section calculations and assign an 
uncertainty based on the differences. The results are 
shown in Table IXIIII The largest difference is seen 
in the acceptance for the Z — > fifi candidate sample, 
which has the narrowest acceptance window in boson 
rapidity. 

F. Uncertainties from PDF Model 

The largest uncertainties on our acceptance values 
arise from uncertainties on the momentum distribu- 
tions of quarks and gluons inside the proton modeled 
in the PDF sets used as inputs to our theoretical 
calculations. The choice of PDF input has a signif- 
icant effect on the shape of the da/dy distributions, 
and consequently a significant effect on the calculated 
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FIG. 14: Acceptance as a function of boson rapidity, A{y), 
for our four candidate samples: W ^ ev (squares), W 
fiv (points), Z ^ ee (stars), and Z —> fifj, (triangles). 



TABLE XIII: Comparison of acceptances for our can- 
didate samples based on da/dy distributions from the 
NNLO and NLO versions of our default theoretical cal- 
culation. The difference is taken as an uncertainty on 
higher-order contributions. 



Acceptance 


NNLO 


NLO 


Difference (%) 




0.1970 


0.1975 


0.25 


Aw^ev 


0.2397 


0.2404 


0.29 


Az^fifi 


0.1392 


0.1402 


0.72 


Az^ee 


0.3182 


0.3184 


0.06 


Az^fj.ii/Aw^fii' 


0.7066 


0.7101 


0.50 


Az-i^el Aw^ev 


1.3272 


1.3246 


0.20 



acceptances for our candidate samples. As noted 
earlier, our theoretical calculations use the best-fit 
MRST 2001 NNLL PDF set [6^. The input PDF 
sets are created by fitting relevant experimental re- 
sults to constrain the parameters which describe the 
quark/gluon momentum distributions in the proton. 
Currently, the NNLL PDF set provided by the MRST 
group is the only one available to us. NLL PDF sets 
are available from both groups (MRSTOIE [H [13] 
and CTEQ6.1 [H]), however. To investigate differ- 
ences between the CTEQ and MRST PDF sets, we 
calculate da/dy at NLO using each group's NLL PDF 



TABLE XIV: Comparison of acceptances for our candi- 
date samples based on da/dy distributions from NLO the- 
oretical calculations using NLL MRST and CTEQ PDF 
sets. 

Acceptance MRST CTEQ Difference (%) 
Aw-,^„. 0.1976 0.1960 0.82 

Aw-,^,. 0.2405 0.2385 0.84 

Az^^^, 0.1401 0.1376 1.82 

Az-.^e 0.3183 0.3164 0.60 

Az-,^,^/Aw-,^^ 0.7088 0.7021 0.95 
Az-,ee/Aw-.e^ 1.3235 1.3264 0.22 



set and check for differences in the acceptance values 
for our candidate samples based on each calculation. 
The results are shown in Table IXIVI The differences 
are significant, especially for the Z — > /i/i candidate 
sample. 

Another recent development from the CTEQ and 
MRST groups is the release of "error" PDF sets 
at NLL which map out the space of potential PDF 
parameter values based on the uncertainties of the 
experimental results used to constrain them. The 
CTEQ PDF parameterization is based on 20 param- 
eters. Pi, which are tuned to their most likely values 
based on a minimization of the of a global fit to the 
experimental data. The equivalent MRST parame- 
terization uses only 15 parameters. As the covariance 
matrix of the Pi is not diagonal at the minimum, it is 
difficult to propagate fit errors on the Pi directly into 
uncertainties on experimentally measured quantities 
such as acceptances. However, both groups construct 
different sets of eigenvectors, Qi, which do diagonalize 
the covariance matrix of the fit in the vicinity of the 
minimum. The Qi are linearly independent by design, 
which allows experimental uncertainties based on de- 
viations in each parameter to be added in quadrature. 
The MRST and CTEQ groups transform individual 
± 1 cr variations of each Qi back into the Pi parame- 
ter space and generate sets of "up" and "down" error 
PDFs. This procedure outputs two PDF sets per pa- 
rameter for a total of 40 CTEQ (30 MRST) error 
PDF sets. The ± 1 ct variations of each eigenvec- 
tor for the MRSTOIE and CTEQ6.1 error PDFs are 
different. These variations are based on the following 
values for the global fit from its minimum: Ax^ = 
50 for MRSTOIE and Ax^ 100 for CTEQ6.1. 

To determine the uncertainty on the acceptance 
values for our candidate samples based on the CTEQ 
and MRST error PDF sets, we perform the NLO 
da/dy production cross section calculations for each 
error PDF set and check how much the acceptance 
values based on each calculation deviate from the val- 
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ues obtained using the best-fit NLL PDF set. The un- 
certainty associated with each Qi is determined from 
the changes in acceptance between the best-fit PDF 
set and both the "up" and "down" error PDF sets as- 
sociated with the given parameter, Av4| and A^^. In 
most cases the two acceptance differences he in op- 
posite directions and can be treated independently, 
but in a small number of cases both differences lie in 
the same direction and a different procedure needs to 
be followed. Table IXVl defines both the positive and 
negative uncertainties assigned to the acceptance un- 
certainty for each Qi based on the relative signs of 
AA\ and AA\ . 

The positive and negative uncertainties associated 
with each of the individual Qi (20 CTEQ and 15 
MRST) are summed in quadrature to determine the 
overall PDF model uncertainty on our acceptance val- 
ues. The results of these calculations using both the 
CT EQ a nd MRST error PDF sets are shown in Ta- 
ble IXVII We note that the MRST uncertainties are 
a factor of 2-3 lower than the CTEQ uncertainties 
which is most likely related to different choices for 
the Ax^ values used by the two groups to choose the 
± 1 cr points associated with each of the Qi. We 
choose to use the larger CTEQ uncertainties based 
on the fact that the magnitude of those uncertain- 
ties is more consistent with the differences observed 
between the acceptance values for our samples calcu- 
lated with the best-fit NLL CTEQ and MRST PDF 
sets (see Table IXTV]) . Based on the technique out- 
lined above, we also determine the PDF model un- 
certainties associated with three additional quantities 
useful in the calculation of r(iy) and .g^/(?e detailed 
in Sec. IVIIII These values are given in Table IXVlII 



G. Uncertainties from Boson pr Model 

As discussed in Sec. IV Bl the boson px distribu- 
tions in our PYTHIA simulated event samples are 
tuned based on the CDF Run I measurement of the 
da/dpT spectrum of electron pairs in the mass region 
between 66 GeV/c^ and 116 GeV/c^ (see Sec. 
The simulated j*/Z pr distribution at = 1.8 TeV 
after tuning is shown in Fig. [12] along with the mea- 
sured distribution from Run L The values for the four 
parameters we use in PYTHIA for this tuning are cho- 
sen using a comparison of the Z boson px spec- 
trum measured in Run I and the pythia generated 
spectra obtained while varying the values of our tun- 
ing parameters. 

The acceptance uncertainties related to our boson 
Pt model come primarily from the Run I measure- 
ment uncertainties. We quantify the effect of these 



uncertainties on our measured acceptances using the 
uncertainties returned from the fits used to obtain 
the four pythia tuning parameters. We choose to use 
conservative ± 3 cr fit errors since the fit values for 
each of the tuning parameters, parp(64) in particu- 
lar, are somewhat inconsistent with expectation. We 
study the effects of changes in the boson pT distribu- 
tions on our measured acceptances by re-weighting 
events in the default simulated event samples based 
on differences between the default boson px distri- 
bution and those obtained from ± 3 cr changes in 
our individual tuning parameters. Table |XVllT] sum- 
marizes the best fit values and ± 3 ct variations ob- 
tained for each tuning parameter and the correspond- 
ing acceptance uncertainties for each candidate sam- 
ple. Changes to the parp(93) tuning parameter were 
found to have a negligible effect on the boson pr spec- 
trum and the measured acceptances. Uncertainties 
associated with the other three tuning parameters 
are taken in quadrature to determine an overall ac- 
ceptance uncertainty associated with the boson pT 
model. 



H. Uncertainties from Recoil Energy Model 

An accurate model of the event recoil energy in the 
simulation is important for estimating the acceptance 
of the event criteria applied to our W ^ Iv can- 
didate events. Simulated recoil energy distributions 
are dependent on the models for hadronic showering, 
the boson recoil energy, and the underlying event en- 
ergy. In addition, the simulation used in these mea- 
surements does not model other mechanisms that also 
contribute to the residual recoil energy in data events 
such as multiple interactions and accelerator back- 
grounds. To account for the effects of these differ- 
ences, the simulated recoil energy distributions are 
tuned to match those in the data. 

As discussed in Sec. IV Ci the event recoil energy is 
defined as the total energy observed in the calorime- 
ter after removing the energy deposits associated 
with the high pT leptons in our W ^ Iv and Z ^ U 
candidates. To tune the simulated recoil energy dis- 
tributions, we separate the observed recoil Et in 
each event into components that are parallel and per- 
pendicular to the transverse direction of the highest 
Pt Icpton in the event. The two components, [/||'"^' 

and Vf^"^, are each assigned energy shift (C) and 
scale {K) corrections in the form 

([/|pi)' = (i^ii X J7|P') +q|, (12) 
iJJ'r^)' = {K^-kUT^)+C^. (13) 
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TABLE XV: Contributions to the positive and negative acceptance uncertainties based on acceptance differences 
between tfie "up" and "down" error PDF sets associated with a given Qi and the best-fit PDF set. 



Direction of Acceptance Shifts + Uncertainty — Uncertainty 

AA\ > and AA\ > y^(AA'^^ + Ay4'^^)/2 

AA\ > and AA\ < AA] AA\ 

AA\ < and AA\ > A^'^ AA\ 



AA'i < and AA\ < V (^^\ + )/2 



TABLE XVI: PDF model acceptance uncertainties based on the CTEQ and MRST error PDF sets. 





CTEQ 


CTEQ 


MRST 


MRST 


Acceptance 


+ Uncertainty 


— Uncertainty + Uncertainty 


— Uncertainty 




(%) 


(%) 


(%) 


(%) 


Aw^fiv 


1.13 


1.47 


0.46 


0.57 


Aw^ev 


1.16 


1.50 


0.48 


0.58 


Az^fj,/! 


1.72 


2.26 


0.67 


0.87 


Az^ee 


0.69 


0.84 


0.27 


0.33 


Az^fiij,/ Aw^iiv 


0.67 


0.86 


0.26 


0.31 


Az^ee/Aw^ei^ 


0.74 


0.56 


0.29 


0.23 



The scale corrections are used to account for prob- 
lems in the calorimeter response model and the effects 
of multiple interactions, the underlying event model, 
and accelerator backgrounds which are in principle 
independent of the lepton direction. The shift cor- 
rections are designed to account for simulation defi- 
ciencies that have a lepton-direction dependence such 
as the W boson recoil model and the model for lep- 
ton energy deposition in the calorimeter. Based on 
the nature of these effects, we expect that the scaling 
corrections in both directions, i^|| and K±, should be 
equivalent and that the shift correction in the per- 
pendicular direction, C±, should be zero. We check 
these assumptions, however, by keeping each param- 
eter independent in the fitting procedure used to de- 
termine the best values for tuning the recoil energy 
in simulated events. 

To determine the best values for these scaling and 
shifting constants, we perform fits between the 
data recoil energy distributions and corrected distri- 
butions from the simulation based on a range of scal- 
ing and shifting constants. An iterative process is 
used in which we first determine the best possible 
shifting constants and then fit for scaling constants 
based on those values. We repeat this process until 
the fits for both the scaling and shifting constants 
stabilize at set values. No effects are expected which 
can give rise to shifts in the energy perpendicular to 
the lepton momentum and the fitted shifts of these 



distributions are consistent with zero. We set C± to 
zero. We also find that the fitted scale factors for 
both recoil energy components agree well with each 
other in both the electron and muon candidate sam- 
ples. Based on this agreement, we also make a com- 
bined fit to both components for a single correction 
scale factor. We use this single scaling factor to cor- 
rect both recoil energy components. A comparison of 
the C/p'^' and U'^'^^ distributions for W ^ ev candi- 
date events in tuned simulation and data are shown 
in Figs. [15] and [TH 

The uncertainties on our measured acceptances re- 
lated to the recoil energy model in the simulation are 
estimated using the ± 3 a values of the scale and 
shift correction factors returned from our fit proce- 
dure. As in the case of boson pT model uncertain- 
ties, we choose to use the ± 3 cr values rather than 
the ± 1 (T values as we are using these parameters 
to cover a wide range of effects that are potentially 
incorrectly modeled in our simulated event samples. 
Since the tuning parameters are not directly related 
to the underlying mechanisms that affect the recoil 
energy distributions, we choose to be conservative in 
how we estimate the associated acceptance uncertain- 
ties via this procedure. We recalculate the acceptance 
of our candidate samples with each of the individual 
tuning parameters changed to its ± 3 ct values and as- 
sign an uncertainty based on the differences between 
these results and our default acceptance values. The 
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TABLE XVII: Additional PDF model acceptance uncertainties based on the CTEQ and MRST error PDF sets. 





CTEQ 


CTEQ 


MRST 


MRST 


Acceptance 


+ Uncertainty 


— Uncertainty + Uncertainty - 


~ Uncertainty 




(%) 


(%) 


(%) 


(%) 


aw ■ Aw^ij.v/cyz 


^Z^MM 1-03 


1.06 


0.52 


0.42 


aw ■ Aw^ev/az 


Az_ee 0.70 


1.06 


0.42 


0.62 


Aw^ev / Aw^fiv 


0.03 


0.04 


0.01 


0.01 



TABLE XVIII: Fit results for pythia boson pt tuning parameters and corresponding uncertainties on the measured 
acceptances of our candidate samples. 



± 3 a 

Parameter Best Fit Variation AAw^^v AAw^ei^ AAz-,^ti AAz-,ee 
(%) (%) (%) (%) 



PARP(62) 


1.26 


0.30 


0.01 


0.00 


0.01 


0.01 


PARP(64) 


0.2 


0.03 


0.03 


0.04 


0.08 


0.06 


PARP(91) 


2.0 


0.3 


0.02 


0.02 


0.00 


0.02 


PARP(93) 


14 


3 










Combined 
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FIG. 15: Comparison of Uf^'^'^ recoil energy distributions 
for W —* ev candidate events in tuned simulation and 
data. 



FIG. 16: Comparison of recoil energy distributions 

for W —* ev candidate events in tuned simulation and 
data. 



changes in acceptance found from modifying the over- 
all scale correction K and the shift corrections for 
both directions, C\\ and C±, are added in quadrature 
to estimate the total uncertainties on our measured 
acceptances due to the recoil energy model in the sim- 
ulation. To be conservative we choose to include an 
uncertainty based on fit results for Cj^ even though 
this parameter is set to zero for tuning the simulated 
recoil energy distributions. 

Table lXIXl summarizes the best fit values and ± 3 cr 



variations with respect to the best fit values obtained 
for each of the scaling and shifting parameters used 
to tune the recoil energy model in simulation and the 
corresponding acceptance uncertainties for the W — > 
Iv candidate samples. 
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TABLE XIX: Summary of simulation recoil energy tuning parameter values and uncertainties obtained from our fit 
procedure and the corresponding uncertainties on our measured acceptance values. 



Tuning W ^ ef W 


^ eu W ^ fiv W 






Parameter Fit Value ± 3 ct 


variation Fit Value ± 3 cr 


variation AAw^eu 










(%) 



/ClI 1.06 0.02 1.06 0.03 

K± 1.04 0.02 1.05 0.02 

K 1.05 0.02 1.05 0.02 0.17 0.20 

C|| -0.4 0.1 -0.1 0.1 0.18 0.29 

C± 0.0 0.1 0.0 0.1 0.00 0.00 

Combined 0.25 0.35 



I. Uncertainties from Energy and Momentum 
Scale/Resolution 

The modeling of COT track pr scale and resolution 
in the simulation affects our acceptance estimates for 
the minimum track px requirements made on muon 
and electron candidates in our samples. Similarly, the 
model of cluster Et scale and resolution for the elec- 
tromagnetic sections of the calorimeter can change 
the acceptance estimates for the minimum cluster Et 
requirements on electrons. Lepton energy and mo- 
mentum measurements can also alter the event ]frp 
calculation, and incorrect modeling of these quanti- 
ties can therefore also affect our acceptance estimates 
for the minimum ]frp criteria applied to our W £v 
samples. 

We check the scale and resolution of the track px 
and cluster Et measurements in the simulation using 
the invariant mass distributions of /Z fifj, and 
"f* /Z —^ ee candidate events. A direct comparison of 
these distributions in data and simulation is possible 
due to the small level of background contamination 
in these samples. We first define scale factors for 
COT track pT (Kpj.) and cluster Et {Ket) in the 
simulation via the expressions 



p't 



E'rp 



PT X PT, 

Ket X Et- 



(14) 
(15) 



The best values for these scale factors are determined 
by making a series of fits between the ^* /Z U 
invariant mass distributions in data and tuned simu- 
lation based on a range of values for the scale factors. 
The best fit for the track pT scale factor is Kp^ = 
0.997. Since the mean of the ^* /Z invariant 
mass peak in the simulation is centered on the mea- 
sured Z boson mass, the best fit value for i^Tp^ is 
indicative of the fact that the current pt scale for 
reconstructed tracks in the data is low. This result 
is consistent with track pT scaling factors for simula- 
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FIG. 17: 7*/^ ^ /i/i invariant mass distribution in data 
and tuned simulation normalized to the data. The ar- 
rows indicate the invariant mass range of our ^* /Z cross 
section measurement. 



tion obtained from similar fits to the J/^P ^ /i/i and 
T — !■ fifi invariant mass peaks indicating that the re- 
sulting scale factor is not pt dependent. The best 
fits for the cluster Et scale factors in the central and 
plug calorimeter modules are -R'e.j, (central) = 1.000 
and Ket (plug) = 1.025 indicative of a model that un- 
derestimates energy deposition in the plug modules 
but is accurate for the central module. 

Comparisons of the 'y* /Z U invariant mass dis- 
tributions in the data and the simulation are used 
to tune the track pT and cluster Et resolution in 
the simulation. We smear these values in simulated 
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using the ± 3 tr values to estimate acceptance uncer- 
tainties is conservatively based on the idea that these 
tuning parameters are not directly related to the un- 
derlying mechanisms that set the scale and resolu- 
tion of track px and cluster Et in the detector. The 
acceptance uncertainties are estimated by observing 
the changes in measured acceptance for each candi- 
date sample that occur when each individual tuning 
parameter is changed between its default and ± 3 ct 
values. A summary of the fitted values and uncertain- 
ties of the scale and resolution tuning parameters for 
track pt and cluster Et is given in Table IXXI along 
with the estimated uncertainties on the measured ac- 
ceptances of our candidate samples associated with 
each parameter. 



40 50 60 70 80 90 100110120130 
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FIG. 18: y* /Z — > ee invariant mass distribution in data 
and tuned simulation normalized to the data. The ar- 
rows indicate the invariant mass range of our y* /Z cross 
section measurement. 



events by generating a random number from a Gaus- 
sian distribution with mean equal to one and width 
equal to a chosen a for each lepton candidate in our 
samples. The resolution smearing is obtained by mul- 
tiplying the track pT and/or cluster Et by the dif- 
ferent random numbers obtained from our distribu- 
tion. Setting a equal to zero adds no smearing since 
each generated random number equals one by defi- 
nition. The best values for a are obtained from 
fits between the y* /Z ££ invariant mass distribu- 
tions in data and tuned simulation corresponding to 
a range of values for a. The best fits for track pT 
and central calorimeter Et resolution are found to be 
for the case of a equal to zero indicating that these 
resolutions are well-modeled in the simulation. The 
best fit for plug calorimeter Et resolution is for a 
value of a above zero indicating that the simulation 
model for Et resolution in the plug modules needs to 
be degraded to match the data better. Fig. [T7] and 
Fig. [T51 show comparisons between the j*/Z ^fj, 
and y* /Z — !■ ee invariant mass distributions in data 
and tuned simulation. 

The effects of uncertainties in the simulation model 
for the scale and resolution of track pT and cluster Et 
on our measured acceptances are estimated based on 
the ± 3 (T values of the corresponding tuning param- 
eters obtained from our fit procedure. Our choice of 



Uncertainties from Detector Material Model 

The acceptances of the kinematic selection crite- 
ria applied to electron candidates are dependent on 
the amount of material in the detector tracking vol- 
ume since electrons can lose a significant fraction 
of their energy prior to entering the calorimeter via 
bremsstrahlung radiation originating from interac- 
tions with detector material. 

The electron E/p distribution, because of its sen- 
sitivity to radiation, is used to compare the material 
description in the detector simulation in the central 
region with that of the real detector as observed in 
data. One measure of the amount of material that 
electrons pass through in the tracking region is the ra- 
tio of the number of events in the peak of the E/p dis- 
tribution (0.9 < E/p < 1.1) to the number of events 
in the tail of the distribution (1.5 < E/p < 2.0). 
We study the uncertainty in the amount of mate- 
rial in the simulation by varying the thickness of a 
cylindrical layer of material in the detector simula- 
tion geometry description in the region between the 
silicon and COT tracking volumes. We choose to use 
copper as the material for this cylindrical layer as 
it best describes the silicon tracker copper readout 
cables and is also supported by independent stud- 
ies of muon energy loss in the calorimeter. Based 
on electron candidates produced in decays of both 
W and Z bosons, we determine that the matching 
of the E/p distribution between data and simulation 
has an uncertainty corresponding to ± 1.5 % of a ra- 
diation length (Xq) of copper. This variation in the 
thickness of the cylindrical layer is used to model the 
acceptance uncertainties originating from the model 
of the detector material in the simulation. 

This result is cross-checked by counting the frac- 
tion of electrons inW ^ ev candidate events which 
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TABLE XX: Summary of simulation track pr scale and resolution tuning parameters and corresponding uncertainties 
on our measured acceptance values. 



Tuning Z — > /i/x Z — * fifi Z ^ ee Z ^ ee 



Parameter 


Fit Value ± 3 cr variation Fit Value ± 3 


a variation AAw^eu 




AAz->ee 










(%) 


(%) 


(%) 


(%) 




0.997 0.003 




0.03 


0.21 


0.04 


0.05 




1.000 0.003 




0.00 


0.00 


0.00 


0.00 


Kej. (Central) 


1.000 


0.003 


0.34 


0.00 


0.23 


0.00 


(Central) 


1.000 


0.015 


0.03 


0.00 


0.05 


0.00 


Ke^ (Plug) 


1.025 


0.006 


0.00 


0.00 


0.11 


0.00 


aET (Plug) 


1.027 


0.011 


0.00 


0.00 


0.05 


0.00 
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FIG. 19: Comparison oi E/p distribution for electron can- 
didates in Z ee events in data and simulation. The 
± 1 CT samples are simulated with ± 1.5 % of a radiation 
length of copper in the tracking volume. 



form "tridents" (see Sec. IVII|) . The probability of 
finding a trident, created when an electron radiates a 
photon which immediately converts into an electron- 
positron pair, is strongly dependent on the amount 
of material traversed by the electron inside the track- 
ing volume. We also compare the resolution of the 
Z — > ee invariant mass peak in data and simulation 
which is sensitive to the rate of radiative interactions 
within the tracking volume. The results of these stud- 
ies are consistent with the E/p results. 

Fig. [19] shows the E/p distributions for electron 
candidates in our Z — > ee data and simulated event 
samples. The ± 1 cr material samples are simulated 
using ± 1.5 % of a radiation length of copper. We 
observe good agreement between data and our default 
simulation in the region below E/p — 2.5. In the high 
E/p tail above this value, the comparison is biased 
by dijet background events in the data. 



The tracks associated with electron candidates in 
the calorimeter plug modules have a low reconstruc- 
tion efficiency due to the limited number of tracking 
layers in the forward region. Therefore, the plug pre- 
radiator detector is used to study the detector ma- 
terial in the simulation for plug electron candidates. 
The amount of energy deposited in the plug prera- 
diator depends on the shower evolution of the elec- 
tron in front of the calorimeter which is itself depen- 
dent on the amount of material the electron passes 
through before entering the calorimeter. On aver- 
age, electrons passing through more material inside 
the tracking volume will have more evolved showers 
at the inner edges of the calorimeter and therefore 
deposit more energy in the plug preradiator. 

To study the detector simulation material descrip- 
tion in the forward part of the tracking volume, we 
compare the ratio of energies observed in the plug 
preradiator and remaining plug calorimeter sections 
for forward electron candidates in data and simula- 
tion. As in the central region, we study the mate- 
rial in the simulation by varying the thickness of an 
iron disk in the volume between the tracking chamber 
endplate and the inner edge of the plug calorimeter. 
These studies indicate that our model for detector 
material in the forward region has an uncertainty cor- 
responding to ± 16.5% of a radiation length {Xq) in 
the iron disk. Fig. [20] shows the ratio of energies ob- 
served in the plug preradiator (PPR) and the plug 
electromagnetic calorimeter (PEM) for electron can- 
didates in data and simulation as a function of both 
the combined energy (PPR -I- PEM) and pseudora- 
pidity of the candidates. 

Acceptance uncertainties coming from the simu- 
lation material model are determined by generating 
simulated event samples with the thicknesses of the 
extra material layers set one at a time in the simula- 
tion to the lower and upper limits of their uncertainty 
ranges. The changes in measured acceptance for the 
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FIG. 20: Comparison of observed ratio of energies in plug preradiator and plug electromagnetic calorimeter for electron 
candidates as a function of the combined energy (left) and pseudorapidity (right) of candidates. Data distributions are 
denoted by the filled circles. The open triangles and associated shaded band show the distribution and uncertainty 
range obtained from simulation when ± 16.5 % of a radiation length thick iron disk is used in the detector material 
description. 



TABLE XXI: Summary of acceptance uncertainties due 
to detector tracking volume material model in simulation. 



Material Model A^w-ei. 


AAz->ee 


Central 0.73 % 


0.94 % 


Plug 


0.21% 



W ^ ev and ^* /Z ee samples relative to the de- 
fault simulation for the modified detector material 
models are summarized in Table IXXII 



K. Acceptance Uncertainty Summary 



The acceptance uncertainties on our event samples 
are summarized in Table IXXIII 



VI. EFFICIENCY 
A. Introduction 

The acceptance values estimated from our simu- 
lated samples are corrected for additional inefficien- 
cies from event selection criteria that are either not 
modeled in the simulation or are better measured 
directly from data. We determine a combined effi- 
ciency, etot, for each candidate sample based on mea- 
sured efficiencies for the individual selection criteria. 
We account for correlations between different selec- 
tion criteria by having a specific order in which in- 
dividual efficiency measurements are made. The effi- 
ciency measurement for a given selection criterion is 
made using a subset of candidates that passes the full 
set of selection criteria ordered prior to the one being 
measured. In addition, since the efficiency is applied 
as a correction to the acceptance, candidates used 
to measure efficiencies are also required to meet the 
geometrical and kinematic requirements used to de- 
fine these acceptances. The ordering and definitions 
of the individual selection criteria efficiencies are pre- 
sented in this introductory section. The following two 
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TABLE XXII: Summary of estimated uncertainties on the measured acceptances for our four candidate samples. 



Uncertainty Category 


lS.Aw—>ev 










/TV \ 
(%) 


(%) 


(%) 


AT/ \ 

(%) 


NNLO da/dy Calculation 


0.29 


0.25 


0.06 


0.72 


PDF Model (positive) 


1.16 


1.13 


0.69 


1.72 


PDF Model (negative) 


1.50 


1.47 


0.84 


2.26 


Boson pt Model 


0.04 


0.04 


0.06 


0.08 


Recoil Energy Model 


0.25 


0.35 


0.00 


0.00 


Track pr Scale/Resolution 


0.03 


0.21 


0.04 


0.05 


Cluster Et Scale/Resolution 


0.34 


0.00 


0.26 


0.00 


Detector Material Model 


0.73 


0.00 


0.96 


0.00 


Simulated Event Statistics 


0.13 


0.14 


0.24 


0.41 


Total (positive) 


1.46 


1.22 


1.23 


1.94 


Total (negative) 


1.75 


1.57 


1.26 


2.44 



sections describe how these individual efficiencies are 
combined to obtain the total event efficiencies for our 

W ^ £v and Z ^ U candidate samples. The re- 
maining sections describe how each of the individual 
efficiency terms is measured. 

The first efficiency term is Cvtxj the fraction of pp 
collisions that occur within ± 60 cm of the center of 
the detector along the 2;-axis. We impose this require- 
ment as a fiducial cut to ensure that pp interactions 
are well-contained within the geometrical acceptance 
of the detector. The ^-coordinate of the event ver- 
tex for a given event is taken from the closest in- 
tersection point of the reconstructed high pT lepton 
track(s) with the z-axis. Since event selection crite- 
ria can bias our samples against events originating in 
the outer interaction region, the efficiency of our ver- 
tex position requirement, evtx, is measured directly 
from the observed vertex distribution in minimum- 
bias events. 

We define etrk as the efficiency for reconstructing 
the track of the high pT lepton in the COT and e^^c 
as the efficiency for matching the found track to ei- 
ther a reconstructed electromagnetic cluster in the 
calorimeter (electrons) or a reconstructed stub in the 
muon chambers (muons). The erec term incorporates 
both the reconstruction efficiency for the cluster or 
stub and the matching efficiency for connecting the 
reconstructed cluster or stub with its associated COT 
track. 

For reconstructed leptons (tracks matched to clus- 
ters or stubs) , eid is the efficiency of the lepton iden- 
tification criteria used to increase the purity of our 
lepton samples. To increase the number of events 
in our Z U. candidate samples, we use a looser 
set of identification criteria on the second lepton leg 
in these events. The loose lepton selection criteria 



are a subset of the set of cuts applied to the single 
lepton in W ^ Iv events and the first lepton leg in 
Z ^ ££ events. The combined efficiency for the loose 
subset of cuts is referred to as eud, and we define 
etid as the efficiency for the set of remaining identi- 
fication cuts not included in the loose subset. The 
efficiency of our lepton isolation requirement, which 
helps to reduce non-W/Z backgrounds in our sam- 
ples, is defined independently as eiso- It is important 
to avoid double-counting correlated efficiency losses 
when measuring the efficiencies for our two sets of 
identification cuts and the isolation requirement. We 
eliminate this problem by defining a specific order- 
ing of these terms (eud, fiso, ftid) and measuring each 
efficiency term using the subset of lepton candidates 
that meets the requirements associated with all of 
the efficiency terms ordered prior to that being mea- 
sured. A natural consequence of using this procedure 
is that the total lepton identification efficiency, eia, is 
necessarily equal to the product of end and etid. 

As discussed previously, the high px electron and 
muon data samples used to make the production cross 
section measurements are collected with lepton-only 
triggers. We define etrg as the efficiency for an iso- 
lated, high quality reconstructed lepton to have sat- 
isfied all of the requirements of the corresponding 
lepton-only trigger path. CDF has a three-level trig- 
ger system, and the value of etrg is determined from 
the product of the efficiencies measured for each of 
the levels. The measured efficiency for a specific level 
of the trigger is based on the subset of reconstructed 
track candidates that satisfy the trigger requirements 
of the levels beneath it. This additional requirement 
is made to avoid double-counting correlated losses in 
efficiency observed in the different trigger levels. 

Finally, there are two efficiencies that are applied 
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only in measurements made in the muon decay chan- 
nels. We define Ecos as the efficiency for signal events 
not to be tagged as cosmic ray candidates via our tag- 
ging algorithm. The cosmic ray tagging algorithm is 
not based on the properties of a single muon, but 
rather on the full set of tracking data available from 
the COT in each event. As a result, ecos is determined 
as an overall event efficiency rather than an addi- 
tional lepton efficiency. Due to topological differences 
between W ^ and Z — > /i/i events, the fraction 
of signal events tagged by the algorithm as cosmic 
rays is different for the two candidate samples. We 
refer to the efficiency term for the W ^ fiv sample 
as eJJ^g and that for the Z — > /i/x sample as e^^. One 
additional event selection made only in the case of 
our W ^LV candidate sample is the Z-rejection cri- 
teria. Due to the non-uniform coverage of the muon 
chambers, we find cases in which only one of the two 
high muon tracks originating from a Z-boson de- 
cay has a matching stub in the muon detector. The 
additional selection criteria made to eliminate these 
events from our W ^ candidate sample has a 
corresponding efficiency defined as ez_iej- 

B. W ^ (.u Efficiency Calculation 

The efficiency of detecting a. W ^ Iv decay that 
satisfies the kinematic and geometrical criteria of our 
samples is obtained from the formula shown in Eg. [TBI 

£tot — £vtx X etrk X ercc X Ejcj 

XCiso X Etig X Ez— rcj X e^.^^ (16) 

As described in detail above, the ordering of the 
cuts, as shown by their left to right order in the for- 
mula, is important. Each efficiency term is an effi- 
ciency for the subset oiW^ii' events that satisfies 
the kinematic and geometric criteria of our samples 
as well as the requirements associated with each of 
the efficiency terms to the left of the term under con- 
sideration. For example, the trigger efficiency term 
in the formula, etrg, is an efficiency for reconstructed 
leptons that satisfy the geometrical, kinematic, iden- 
tification, and isolation criteria used to select the high 
Pt lepton in our W ^ tv candidate events. As noted 
previously, the ez-rc] and terms in the formula 
apply to the W candidate sample only. Ta- 

ble [XXlTT] summarizes the measurements of the indi- 
vidual efficiency terms (described in detail below) and 
the resulting combined efficiencies for our W Iv 
candidate samples. The electron efficiencies shown 
in Table IXXIIII are for central calorimeter electrons 
only since our W ^ ev cross section measurement is 



also restricted to candidates in this part of the detec- 
tor. 



C. Z U Efficiency Calculation 

For both electrons and muons, we define a loose 
set of lepton selection criteria for the second leg of 
Z ^ II events to increase the size of our candidate 
samples. The efficiency calculation for these samples 
is complicated by the fact that in many events both 
leptons from the Z boson decay can satisfy the tight 
lepton selection criteria which are required for only 
one of the two legs. 

In the electron channel, we allow for two different 
types of loose lepton legs. The second leg can be ei- 
ther a central calorimeter electron candidate passing 
a looser set of selection criteria or an electron recon- 
structed in the forward part of the calorimeter (plug 
modules). For Z — > /i/i candidates, a loose track leg 
is not required to have a matching reconstructed stub 
in the muon detectors. For this sample, the second 
muon leg is simply required to be a high pT-, isolated 
track satisfying the subset of muon identification cuts 
corresponding to the track itself. The breakdown of 
lepton identification cut efficiencies between the loose 
and tight criteria is shown in Table IXXIVI for both 
muons and central electrons. There is no reconstruc- 
tion inefficiency associated with loose muon legs since 
track candidates are not required to have a matching 
muon detector stub. 

Efficiencies for loose plug electrons are given in Ta- 
ble I XXVI There is no track reconstruction compo- 
nent in the plug electron selection efficiency since 
a matched track is not required for candidates in 
the plug region of the calorimeter. Also, since no 
matching between tracks and clusters is done in this 
region, the plug lepton reconstruction efficiency is 
100 %. There are no dead calorimeter towers in the 
data-taking period used in these measurements. We 
also find that kinematic distributions for tight central 
electron legs in our central-plug Z — > ee event sam- 
ple are somewhat different from those in the central- 
central sample. These kinematic differences have a 
small effect on the electron identification efficiencies 
for the central legs in central-plug Z — > ee events. In 
order to correct for this effect, we measure a central 
leg scale factor, S"^/"®, which is the ratio of central leg 
efficiencies in central-plug Z — > ee events to those in 
central-central events. The value of this scale factor 
given in Table IXXVI is determined from simulation 
and is applied as an extra term in the overall selec- 
tion efficiency for plug electrons. 

To determine a total event selection efficiency for 
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TABLE XXIII: Summary of the individual efficiency terms for W £iy. 

Selection Criteria Label W ev W ^ fiv 

Fiducial Vertex tvtx 0.950 ± 0.004 0.950 ± 0.004 

Track Reconstruction etrk 1.000 ± 0.004 1. 000 ± 0.004 

Lepton Reconstruction e^oc 0.998 ± 0.004 0.954 ± 0.007 

Lepton ID Qd 0.840 ± 0.007 0.893 ± 0.008 

Lepton Isolation eiso 0.973 ± 0.003 0.982 ± 0.004 

Trigger etrg 0.966 ± 0.001 0.925 ± 0.011 

Z-Rejection Cut e^-^oj - 0.996 ± 0.002 

Cosmic Ray Tagging e^s : 0.9999 ± 0.0001 

Total etot 0.749 ± 0.009 0.732 ± 0.013 



TABLE XXIV: Breakdown of loose and tight lepton iden- 
tification efficiencies. 

Selection Criteria Label Central Electron Muon 
Loose Lepton ID end 0.960 ± 0.004 0.933 ± 0.006 
Tight Lepton ID etid 0.876 ± 0.007 0.957 ± 0.005 
All Lepton ID eid 0.840 ± 0.007 0.893 ± 0.008 



TABLE XXV: Plug electron efliciencies. 

Selection Criteria Label Plug Electron 

Lepton Reconstruction efic^ 1.000 

Lepton ID ef^""^ 0.876 ± 0.015 

Lepton Isolation 6^'"*^ 0.993 ± 0.003 

Central Leg Scale Factor 5^/"*^ 1.014 ± 0.002 

Total eP'"s 0.883 ± 0.015 



Z ^ ee events, we first calculate efficiencies for the 
central-central and central-plug samples which are 
independent of one another by definition. The to- 
tal efficiency is a weighted sum of the efficiencies for 
the two samples. The weighting factors are deter- 
mined from the relative numbers of central-central 
and central-plug events in our simulated sample. The 
fraction of central-plug events, /cp, is determined to 
be 0.655 ± 0.001. Eg. [TTl shows the efficiency calcu- 
lation for central-central Z — s- ee events: 

^tot — ^vtx X Ctrk X ^rcc ^ ^lid ^ ^iso 

X [etid X (2 - etid)] 

x[et,gX (2-et,g)]. (17) 

The squared terms in the formula apply to efficiency 
terms that are applied twice (we require two recon- 
structed central electrons passing loose identification 
and isolation criteria). In order for this treatment to 
be correct, the efficiencies of the two electron legs 



in the Z —* ee candidates are required to be un- 
correlated. Using our sample of simulated Z — > ee 
events, we look for correlations between the efficien- 
cies for the two electron legs and find them to be 
negligible. The tight identification and trigger crite- 
ria can be satisfied by either of the two electrons. 
The combined efficiency for one of two objects to 
satisfy a particular requirement can be written as 
-h 2 X e X (1 - e) = e x (2 - e). The efficiency 
calculation for central-plug Z ^ ee events is given in 
Eq.[18l 

cp 

^tot — ^vtx ^ ^trk ^ ^rcc ^ ^lid ^ ^tid ^ Qso ^ ^trg 

xeP<l^s X eP'-s X Stt'' x e^^ (18) 

In these events only the central electron leg can sat- 
isfy the tight identification and trigger criteria so 
these efficiencies are only applied to the one central 
leg. Similarly, the plug efficiencies are applied only 
to the plug electron leg. Based on Eqs. [T7] and [TH] 
the event efficiency for our combined Z — > ee sample 
takes the form: 

^tot — ^vtx X etrk X ej-oc X eiitj X etid X eiso X etrg 
X [(1 - /cp) ^ ^tik ^rcc ^ ^lid ^iso 

X (2 -etid) X (2 -etrg) 

+/ep X eP',"s X efj"^ X St"^ x ^^\. (19) 

The calculation of the total selection criteria ef- 
ficiency for Z — > /i/i candidate events is similar to 
that for events in the electron channel but involves 
some additional compHcations. As discussed above 
we increase our acceptance for Z — > [L[i events by 
releasing the muon detector stub requirements for 
one of the two candidate track legs. The second 
muon leg in our candidate events can be any COT 
track passing the track quality, isolation, and mini- 
mum ionizing calorimeter energy deposition criteria 
used in this analysis for selecting muon track candi- 
dates. Since the track selection criteria are applied 
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to both muon legs in our candidate events, the cor- 
responding terms in the overall efficiency formula are 
squared. Only one of the two muon track candidates 
is required to have a matching stub in the muon de- 
tectors that satisfies our stub selection criteria. For 
roughly 40 % of our candidate events, both of the 
muon track legs point to active regions of the muon 
detectors. In these cases, either of the two legs can 
have a matching stub in the muon detectors and sat- 
isfy the tight leg criteria. In other cases, one of the 
two legs will not point to an active detector region, 
and the stub-matching criteria must be satisfied by 
the one leg that is pointed at the muon detectors. In 
order to determine the total efficiency for Z ^ fjLfi 
candidate events, we first determine the total selec- 
tion efficiencies for both of these event classes. The 
event selection efficiency for the combined sample is 
then extracted as the weighted sum of the efficiencies 
for the two different event types. 

The efficiency calculation for the subset of Z fifi 
events in which only one of the two muon tracks 
points to an active region of the muon detector is 
shown in Eq. [501 The efficiencies corresponding to 
selection criteria applied to both muon legs (track 
reconstruction, loose identification, and isolation) en- 
ter into the formula as squared terms. The track leg 
pointing at the inactive regions of the muon detector 
can not have an associated reconstructed stub so the 
other track leg in the event must have a matching 
stub for the event to satisfy the Z ^ nn selection 
criteria. This leg must also satisfy the tight muon 
identification and event trigger requirements since an 
associated reconstructed muon detector stub is a nec- 
essary pre-condition for a muon leg to meet these cri- 
teria. Since the muon stub reconstruction, tight iden- 
tification, and trigger selection criteria can only be 
satisfied by one of the two muon legs in these events, 
the corresponding efficiency terms enter into Eq. [501 
linearly. As previously mentioned, the efficiency for 
Z — !■ fi^ events not to be tagged as cosmics, e^g, is in- 
dependent of the measured value for W ^ events. 
As described subsequently, we measure this efficiency 
in Z ^ events to be ef^s = 0.9994 ± 0.0006. 

XErcc X etid X etrg (20) 

This situation is more complicated for the class of 
Z — > /i/i events where both muon legs point to active 
regions of the muon detector. For these events both 
legs can individually satisfy the stub reconstruction, 
tight identification, and trigger criteria of the sam- 
ple. In order to simplify the efficiency calculation, 
we require that at least one of the two muon legs in 



each candidate event satisfies the requirements asso- 
ciated with all three of the above criteria. With this 
additional restriction, the overall event selection ef- 
ficiency in the subset of Z — > /i/x candidates where 
both muon legs point at active regions of the muon 
detector can be written as shown in Eq. [^T] The 
combined efficiency for a muon leg to satisfy the stub 
reconstruction, tight identification, and trigger crite- 
ria (e* — free X etid X etrg) enters into Eq. [21] in the 
form e* x (2 — e*) which, as described above, is the 
resulting efficiency for a set of criteria required for 
one of two identical objects within an event. 

Hot — ^vtx X e(.Qg X eti.]j X eij(j x 6;^^ 

X [erec X etid X etrg 

X (2 - erec X etid X etrg)] (21) 

In order to combine Eos. [^ and [^ into a formula 
for the total event efficiency of our combined sample, 
we need to introduce an additional parameter, /dd, 
which is defined as the fraction of Z — > /i^ events 
within our geometric and kinematic acceptance in 
which both muon legs are found to point at active 
regions of the muon detector. This quantity is de- 
termined from the simulated event sample. For our 
candidate sample, we obtain /dd = 0.3889 ± 0.0021, 
which is a luminosity weighted average of the values 
for the different run periods in which the CMX was 
either ofHine or online. Using this additional factor, 
we determine a formula for the total event efficiency 
of our candidate sample by adding the expressions in 
Eqs. [^ and [^ weighted by factors of 1 — /dd and /dd 
respectively. Finally, we obtain the expression shown 
in Eq. [22| for the total selection efficiency for events 
in our Z jj,^ candidate sample. 

Z^p^ _ ^ ^ y, 2 X ,2 „ ,2 

^tot — ^vtx X e^Qg X etrk X eji^j x ej^o 

Xercc X etid X etrg X [1 + /dd 

x(l - eroc X etid X etrg)] (22) 

Based on the expressions in Eqs. [THj and [211 we 
can substitute our measured values for the individual 
efficiency terms and determine the combined event 
selection efficiencies for our Z ^ i£ candidate sam- 
ples. The resulting values are shown in Table IXXVII 

D. Vertex Finding Efficiency 

Our requirement that the z-position of the primary 
event vertex be within 60 cm of the center of the CDF 
detector (|Zvtx| < 60 cm) limits the event acceptance 
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TABLE XXVI: Results of Z 
ciency calculations. 



It combined event cfii- 



Candidate Sample 








Z ee 


0.713 


± 


0.012 


Z 


0.713 


± 


0.015 



to a fraction of the full luminous region for pp colli- 
sions. However, the luminosity estimate used in our 
cross section measurements is based on the full lu- 
minous range of the beam interaction region. We 
use minimum-bias data to measure the longitudinal 
profile of the pp luminous region, and this profile is 
subsequently used to estimate the fraction of interac- 
tions within our fiducial range in z. 

Fig. [21] shows the distribution of measured posi- 
tions along the z-axis (parallel to beams) for recon- 
structed primary vertices in minimum-bias events. 
The minimum-bias events are taken from the same 
set of runs from which our candidate samples are 
constructed. In addition, the minimum-bias data is 
weighted to ensure that it has the same run-by-run 
integrated luminosity as the cross section event sam- 
ples. We fit the distribution in Fig. [21] using the fol- 
lowing form of the Tevatron beam luminosity func- 
tion: 
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The five free parameters of the fit are TVq, <Jz, ^oii 
zo2, and (3* . The Zvtx distribution has some biases 
at large values of z due to increased contamination 
from Tiow-pp interactions such as those originating 
from beam-gas collisions and due to the decrease of 
COT tracking acceptance far away from the center 
of the detector. We avoid these biases by only fit- 
ting the measured Zvtx distribution to our function 
for dC{z)/dz in the region where \z\ < 60 cm. Within 
this finite range in z, the fraction of events not from 
pp collisions is negligible and the COT tracking ac- 
ceptance is high and uniform. 

The acceptance of our requirement on the z- 
position of the primary event vertex (|Zvtx| < 60 cm) 
is calculated as 



evtxd^l < 60) 



/-60 [d^z)/dz] dz 
I-^ [dC{z)/dz]dz 



(24) 



We perform the fit to the data and evaluate the ac- 
ceptance for both the full sample and several sub- 
samples of our minimum-bias data set. We observe 
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FIG. 21: The measured Zvtx distribution. The units on 
the horizontal axis are cm and there are a total of 100 
bins from —100 cm to -f 100 cm. The curve is the fit to 
the luminosity function (Eg. I23|) for \z\ < 60 cm, and the 
resulting fit with 55 degrees of freedom has a of 119. 



slight differences in the various sub-samples indicat- 
ing small changes over time in the z-profile of pp 
collisions in the interaction region of our detector. 
The maximum shift seen in the measured acceptance 
among the various sub-samples is 0.6%, and we assign 
half of this value as a systematic uncertainty on the 
efl[iciency measurement. The statistical uncertainty 
on the measurement is assigned based on fit errors 
obtained from the Zvtx fit for the full minimum-bias 
sample. Using the techniques described above, we 
measure the signal acceptance of our cut on the z- 
position of the primary event vertex to be 

e^t^ = 0.950 ± 0.002 {stat.) ± 0.003 (syst.) . 



E. Tracking Efficiency 

We define tracking efficiency as the fraction of high 
Pt leptons contained within our geometrical accep- 
tance for which our offline tracking algorithm is able 
to reconstruct the lepton track from hits observed in 
the COT. We measure this quantity using a sample of 
clean, unbiased W ^ ev candidate events based on a 
tight set of calorimeter-only selection criteria. The 
events for this sample were collected using a trig- 
ger path based on calorimeter Iflrp requirements to 
ensure that the sample is unbiased with respect to 
XFT tracking requirements in the hardware portion 
of the trigger and track reconstruction in the soft- 
ware portion. Events are required to have an elec- 
tromagnetic calorimeter cluster with Et > 20 GeV 
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and overall event If^j, > 25 GeV. Since we can not 
use a track matching requirement to help reduce non- 
electron backgrounds, we apply a tighter than normal 
set of electron identification criteria on the electro- 
magnetic cluster itself. We also remove candidate 
events containing additional reconstructed jets in the 
calorimeter with Et > 3 GeV and require that the 
Pt of the reconstructed W boson is above 10 GeV/c. 
These cuts are designed to remove background events 
in our sample originating from QCD dijet processes. 

Our tracking efhciency measurement is obtained 
from the fraction of events in this candidate sam- 
ple which have a COT track pointing to the electro- 
magnetic cluster. Matching, reconstructed tracks in 
the COT are required to point within 5 cm of the 
calorimeter electromagnetic cluster seed tower. In 
order to be absolutely sure that we are not including 
track-less background events in our efficiency calcula- 
tion, we also require that our candidate events have 
a reconstructed track based entirely on hits in the 
silicon tracking detector (independent of the central 
outer tracking chamber) pointing at the electromag- 
netic cluster. A total of 1368 candidate events in 
our 72.0 pb~^ sample have a matching silicon track. 
Of these, 1363 events also contain a matching re- 
constructed track based solely on hits in the central 
outer tracker yielding a COT tracking efficiency of 
etrkiData) = 0.9963 to.oo40- The uncertainty on the 
measurement is primarily systematic and is based on 
studies of both silicon-only track fake rates and cor- 
related failures in COT and silicon based tracking 
algorithms. 

We compare the tracking efficiency measured in the 
data with an equivalent measurement based on our 
W ^ ev simulated event sample. Using the same 
technique, we obtain a simulation tracking efficiency 
of etrk(MC) = 0.9966 t^^^^i, consistent with our 
measured value from data. A study of failing sim- 
ulated events reveals that the small tracking ineffi- 
ciency we measure is mainly due to bremsstrahlung 
radiation where the silicon-only track points in the 
direction of the hard photon and the COT track fol- 
lows the path of the soft electron (pointing away 
from the high Et electromagnetic cluster). Since 
the loss of events due to this process is already ac- 
counted for in our acceptance calculation, we avoid 
double-counting by taking the ratio of the tracking 
efficiency measured in data to that measured in sim- 
ulation as our net tracking efficiency. Based on this 
approach, our final value for the COT tracking effi- 
ciency is etrk — 1.000 ± 0.004 where the uncertainty is 
the combined statistical and systematic uncertainty 
of our measurement technique. 



F. Reconstruction Efficiency 

The lepton reconstruction efficiency is defined as 
the fraction of real leptons that is within our geo- 
metrical acceptance and has matching, reconstructed 
COT tracks which are subsequently reconstructed as 
leptons by our offline algorithms. In the case of elec- 
trons, this efficiency corresponds to the combined 
probability for forming the electromagnetic cluster 
and matching it to the associated COT track. For 
muons, it is the probability for reconstructing a stub 
in the muon detectors and matching the stub to the 
corresponding COT track. 

The reconstruction efficiency is measured using the 
unbiased, second legs oi Z ^ U decays. The leptons 
from Z boson decays have a similar momentum spec- 
trum to those originating from W boson decays and 
are embedded in a similar event environment. Events 
are required to have at least one fully reconstructed 
lepton leg that satisfies the complete set of lepton 
identification criteria used in the selection of our can- 
didate samples. The same lepton leg must also satisfy 
the requirements of the corresponding high pj- lepton 
trigger path to ensure that the second leg is unbi- 
ased with respect to the trigger. A lepton leg satisfy- 
ing these requirements is then paired with a second 
opposite-sign, high pT track in the event. If the in- 
variant mass of a lepton-track pair lies within the Z 
boson mass window, 80 GeV/c^ < M« < 100 GeV/c^, 
the second track leg is utilized as a candidate for test- 
ing the lepton reconstruction efficiency. In the case 
of Z —> /i/i candidate events only, the second track 
leg is also required to have associated calorimeter en- 
ergy deposition consistent with a minimum ionizing 
particle which reduces backgrounds from fake muons 
without biasing the measurement. In the subset of 
Z U candidate events where each track leg is a 
reconstructed lepton passing the full set of identi- 
fication and trigger criteria, both legs are unbiased 
lepton candidates and included in the efficiency mea- 
surement. 

Each candidate track leg is extrapolated to deter- 
mine if it points at an active area of the calorimeter 
or muon detectors as appropriate. If the track does 
point at an active detector region, it is expected to 
be reconstructed as a lepton. The fraction of this 
subset of candidate tracks which are in fact recon- 
structed as leptons provides our measurement of the 
reconstruction efficiency. Fig. [22] shows the invari- 
ant mass distributions for muon-track pairs in cases 
where the second track is and is not reconstructed 
as a muon. The small peak seen near the Z boson 
mass in the latter case indicates that we do observe a 
non-negligible muon reconstruction inefficiency in the 
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FIG. 22: Invariant mass of muon-track pairs for the muon 
reconstruction-efficiency measurement. We show the dis- 
tribution for pairs in which the track is a reconstructed 
as a muon track (open histogram) and for pairs in which 
the track is not reconstructed as a muon track (sohd 
histogram). Only the region between 80 GeV/c^ and 
100 GeV/c^ is used for the efficiency calculation. 



data. However, some of the measured inefficiency is 
attributable to the effects of multiple scattering. A 
particle associated vifith a track that points at an ac- 
tive detector region can in some cases pass outside 
of this region due to the cumulative effects of inter- 
actions with material in the detector. This effect is 
modeled using the simulated event samples. All real 
reconstruction inefficiencies observable in the simula- 
tion are accounted for in the acceptance calculation 
and must not be double-counted in the lepton re- 
construction efficiency measurement. Therefore, we 
determine our net lepton reconstruction efficiency by 
dividing the value measured in data by the value ob- 
tained from an equivalent measurement using sim- 
ulation. The lepton reconstruction efficiency mea- 
surements for electrons and muons are summarized 
in Table IXXVffI Plug electron candidates are not 
required to have a matching reconstructed track and 
therefore by our definition have a fixed reconstruc- 
tion efficiency of 100 %. We make additional checks 
to confirm that the leptons in our test samples are 
a good match for the leptons in our candidate sam- 
ples and based on this agreement take the statistical 



uncertainty of our measurements as the total uncer- 
tainty on the reconstruction efficiencies. 



G. Lepton Identification and Isolation Cut 
Efficiencies 



The efficiencies of our lepton identification and 
isolation cuts are also determined directly from the 
data using Z ^ ££ events. We use slightly different 
techniques for measuring these efficiencies for elec- 
trons and muons. The motivation for using separate 
methods is the non-negligible fraction of background 
events in the Z ee candidate sample in which at 
least one of the reconstructed electrons is either a fake 
or the direct semileptonic decay product of a hadron. 
In order to accurately measure the selection efficien- 
cies for electrons originating from W and Z boson 
decays, it is important to correct for the contribution 
of these background events to our efficiency calcula- 
tion. Since these types of backgrounds are negligible 
in our Z — > fi^ candidate sample, we are able to use 
a more aggressive approach which maximizes the sta- 
tistical size of the muon candidates used to determine 
these efficiencies. 

As previously mentioned, the identification and 
isolation efficiencies for leptons are determined in a 
specific order to avoid double-counting correlated in- 
efficiencies between different groups of selection cri- 
teria. The order we employ in making these measure- 
ments is the following: efficiencies from loose identifi- 
cation cuts, isolation cut efficiencies, and efficiencies 
from tight identification cuts. This ordering is cho- 
sen to simplify the extraction of combined selection 
efficiencies for Z ££ events from our individual, 
measured efficiency terms. To protect this ordering, 
we require that lepton candidates used to measure 
each group of selection efficiencies satisfy the selec- 
tion criteria associated with all groups defined to be 
earlier within our assigned order. 

To minimize backgrounds in the Z ee event sam- 
ple used to make the efficiency measurements, we re- 
quire that at least one of the two reconstructed elec- 
trons passes the full set of identification and isola- 
tion criteria used in the W ^ ev analysis. The sec- 
ond electron leg in each event, referred to here as the 
probe electron, is simply required to satisfy the geo- 
metric and kinematic cuts that define the acceptance 
of our candidate samples. In addition, the invari- 
ant mass of the electron pair is required to be within 
a tight window centered on the measured Z boson 
mass (75 GeV/c^ < M^e < 105 GeV/c^), which fur- 
ther reduces non-Z backgrounds in the sample. By 
definition, the electron passing the complete set of 
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TABLE XXVII: Summary of lepton reconstruction efficiency measurements. Because plug electron candidates are not 
required to have a matching reconstructed track, the corresponding reconstruction efficiency is one by definition. 

Lepton Data Efficiency Simulation Efhciency Net Efficiency 

Central Electons 0.990 ± 0.004 0.992 ± 0.001 0.998 ± 0.004 

Plug Electrons 1.000 1.000 1.000 

Muons 0.935 ± 0.007 0.980 ± 0.001 0.954 ± 0.007 



identification and isolation criteria is a central elec- 
tron. Central-central Z — > ee events satisfying the 
criteria listed above are used to measure central elec- 
tron efficiencies, and central-plug events are used to 
measure plug electron efficiencies. 

We define the number of central-central Z ^ ee 
candidates passing our criteria as N^c- As mentioned 
above, each event has at least one electron passing the 
full set of identification and isolation criteria. Elec- 
trons of this type are referred to as tight. In some 
number of events in our candidate sample, A'tt, both 
electrons are found to satisfy the tight criteria. In the 
remaining events, the probe electron necessarily fails 
at least one part of our selection criteria. However, 
some number of these remaining events will satisfy a 
particular subset of the identification and isolation re- 
quirements corresponding to an efficiency term that 
we want to measure. The total number of events 
where the probe leg is found to satisfy a given sub- 
set of cuts is referred to as N^. In this case, the 
corresponding efficiency for the subset of cuts being 
studied is determined from the expression given in 
Eq. [25l The variable i in this expression refers to 
the three sets of selection cut efficiencies to be mea- 
sured (1 = loose identification, 2 ~ isolation, and 3 
= tight identification). In the second two cases, we 
limit our sample of probe electrons to those that sat- 
isfy the criteria associated with the lower numbered 
efficiency terms to avoid the double-counting problem 
discussed above. The net result is that for the second 
two cases iVtc = iVt(i_i) and Nu is re-defined as the 
number of events where the probe leg is found to pass 
the requirements associated with the efficiency term 
being measured and those numbered below it. This 
new definition implies that for the final case iVt,; is 
simply equal to A^tt- 



One additional complication is that we must sub- 
tract the contribution of background to each of the 
input event totals in Eq. [21] to accurately measure 
the efficiencies for electrons produced in W and Z 
boson decays. For central-central Z — > ee events, the 



background in each event subset is determined from 
the number of equivalent same-sign events observed 
in the data sample. A correction for tridents (real 
Z — > ee events where the charge of one electron is 
measured incorrectly due to the radiation of a hard 
bremsstrahlung photon) in the same-sign event totals 
is made based on the relative numbers of opposite- 
sign and same-sign events in our Z ^ ee simulated 
event sample. The event counts and background cor- 
rections for each of the input parameters used in the 
efficiency calculations are given in Table IXXVIIII 

The fraction of background events in the central- 
plug Z — > ee candidate sample used to measure plug 
electron efficiencies is much larger than that in the 
central-central sample. In order to eliminate some of 
this additional background, we make an even tighter 
set of requirements on the isolation and electron qual- 
ity variables associated with the central electron to 
pick the candidate events used to measure these effi- 
ciencies. As the probe leg in these candidates is the 
only plug electron of interest in the event, efficien- 
cies are measured simply as the fraction of probe legs 
that satisfy the associated set of selection criteria. In 
the analyses reported here, plug electrons are utilized 
only as loose second legs for selecting Z ^ ee can- 
didate events. There is therefore no corresponding 
tight identification cut efficiency to measure for plug 
electrons. However, the ordering of the loose identifi- 
cation and isolation cuts for plug electrons is identical 
to that used for electrons in the central region. We 
account for this ordering by requiring that the probe 
electrons used to measure the efficiency of the iso- 
lation cut satisfy the full set of loose plug electron 
identification cuts. We correct the number of probe 
legs in both the numerator and denominator of our ef- 
ficiency calculations for the residual backgrounds re- 
maining in our candidate sample. These backgrounds 
are estimated using electron fake rate calculations 
outlined in Sec. IVIII Based on this method, we ob- 
tain independent estimates for the background con- 
tributions from both QCD dijet and W ^ ev plus 
jet processes and sum them to obtain our final back- 
ground estimates. The inputs to our plug electron 
efficiency measurements and the resulting efficiency 
values are summarized in Table IXXIXl 
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TABLE XXVIII: Z — > ee event counts used as inputs to the calculation of electron identification and isolation efficien- 
cies. 



Efficiency Measurement 


Symbol 


i 


iVtc 


Nu 


iVtt 


iVtc Ntt 


Efficiency 














Background Background Background 




Loose Identification Cuts 




1 


1901 


1751 


1296 


28.3 6.1 0.6 


0.960 ± 0.004 


Isolation Cut 


^iso 


2 


1751 


1663 


1296 


6.1 -0.4 0.6 


0.973 ± 0.003 


Tight Identification Cuts 


etid 


3 


1663 


1296 


1296 


-0.4 0.6 0.6 


0.876 ± 0.007 



TABLE XXIX: Input parameters to plug electron identification and isolation efficiency measurements using central-plug 
Z ~* ee candidates. 



Efficiency 
Measurement 


Symbol 


Number of Number passing 
Probe Electrons Selection Criteria 


Probe Electron Passing Electron 
Background Background 


Efficiency 


Identification Cuts 
Isolation Cut 


,plug 

plug 
iso 


2517 2126 
2126 2111 


108.4 15.0 
15.0 14.1 


0.876 ± 0.015 
0.993 ± 0.003 



The calculation of muon identification and isola- 
tion efficiencies is simplified by the lack of signif- 
icant backgrounds in our Z — s- /i/i candidate sam- 
ples. To obtain the largest possible sample of probe 
muons for measuring these efficiencies, we make only 
a minimal set of requirements on the first muon 
leg in these events. In order to avoid selection bi- 
ases, we simply require that at least one muon leg 
in each event satisfies both the trigger requirements 
and loose cuts used to select events into our high 
muon sample from which the candidate events 
are chosen. The second muon leg in each of these 
events is then utilized as an unbiased probe leg for 
measuring our selection efficiencies. In the subset of 
candidate events where both muon legs satisfy the 
trigger and loose selection requirements of our sam- 
ple, both muons are unbiased and included in our 
sample of probe muons. To ensure that we are se- 
lecting probe muons from a clean (low background) 
sample of Z — s- /i/i candidate events, we do require 
that the invariant mass of each muon pair lies within 
a tight window around the measured Z boson mass 
(80 GeV/c^ < M^,^, < 100 GeV/c^) and remove any 
events identified by our tagging algorithm as cosmic 
ray candidates. After applying these criteria, we find 
that only 3 of over 1,500 probe muons come from 
same-sign candidate events confirming the negligi- 
ble background fraction in the event sample used for 
these measurements. 

As in the case of electrons, the full set of muon 
identification cuts is divided into loose and tight sub- 
sets to simplify the calculation of the combined event 
selection efficiency for Z — > fifj, candidate events. The 
second muon track leg in these events is not required 
to have a matching stub in the muon detector. There- 



fore, the identification cuts for muons which we refer 
to as loose are those that are applied to the track 
itself. The remaining tight selection cuts are those 
applied only to muon track legs with matching muon 
detector stubs. In some sense, the reconstruction of 
a matching stub in the muon detector is therefore 
also a tight selection criteria although we choose to 
treat the efficiency for this requirement separately. 
We use the same ordering of selection criteria (loose 
identification, isolation, and tight identification) as 
that used for electrons to avoid the double-counting 
of correlated muon inefficiencies. Muon probe legs 
used to measure the efficiency for each set of selec- 
tion criteria are required to satisfy all selection cuts 
corresponding to previously ordered efficiency terms. 
Table IXXXI summarizes the inputs to the muon effi- 
ciency calculations and the resulting efficiency values. 



H. Trigger Efficiency 

As described in Sec. IIIIl the data samples used to 
select our candidate events are collected via high pT 
lepton-only trigger paths. The three-level trigger sys- 
tem utilized by the upgraded CDF data acquisition 
system reduces the 2.5 MHz beam-interaction rate 
into a final event collection rate on the order of 75 Hz. 
The first two levels utilize dedicated hardware to se- 
lect events for readout from the detector, and the 
third level is a processor farm that runs a fast ver- 
sion of the full event reconstruction to pick out the 
final set of events to be written to tape. Level 1 lep- 
ton triggers are constructed from high pT COT tracks 
identified in the fast tracking hardware matched with 
single tower electromagnetic energy deposits in the 
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TABLE XXX: Input parameters to muon identification and isolation efficiency measurements using Z fifj, candidates. 



Efficiency Measurement 


Symbol 


Number of Number passing 
Probe Muons Selection Criteria 


Efficiency 


Loose Identification Cuts 


eiid 


1574 1469 


0.933 ± 0.006 


Isolation Cut 


f^iso 


1469 1443 


0.982 ± 0.003 


Tight Identification Cuts 


etid 


1443 1381 


0.957 ± 0.005 



calorimeter (electrons) or groups of hits in the outer 
wire chambers (muons). Level 2 hardware is used 
to perform a more sophisticated calorimeter energy 
clustering algorithm on electron candidates to obtain 
improved Et resolution. The improved Et variable is 
utilized at Level 2 to make tighter kinematic cuts on 
the electron candidates. No additional requirements 
are made on muon candidates at Level 2. Events se- 
lected at Level 2 are read out of the detector and 
passed to the Level 3 processor farm. A fast ver- 
sion of the offline lepton reconstruction algorithms 
are run on each event, and the identified leptons are 
subjected to both kinematic and loose quality selec- 
tion cuts. 

The measurement of trigger efficiencies for elec- 
trons is simplified by the availability of secondary 
trigger paths that feed into our W ^ ev candidate 
sample. A trigger path based solely on calorimeter 
quantities is used to measure the efficiency of track- 
ing requirements at each of the three trigger lev- 
els. This path utilizes identical calorimeter cluster 
requirements to those in the default electron path 
but does not require matching tracks to be found at 
any level. Instead, events are selected based on the 
presence of large $rp in the calorimeter (15 GeV at 
Level 1/Level 2 and 25 GeV at Level 3) associated 
with the high energy neutrino in the W boson decays. 
For W ^ nv candidate events, the muon deposits 
only a small fraction of its energy into the calorime- 
ter and hence the residual If]rp in the calorimeter is 
too small to allow for an equivalent trigger path for 
muon candidates. To measure the efficiencies of the 
electron trigger path track requirements, we select 
events from the secondary trigger path that pass the 
complete set oiW ^ ev selection criteria. The frac- 
tion of events in this unbiased sample that satisfy the 
track requirements of our lepton-only trigger path at 
each of the three levels gives the corresponding effi- 
ciency for those requirements. The double-counting 
of correlated inefficiencies between the different trig- 
ger levels is avoided by requiring that events used 
to measure higher level trigger efficiencies pass all of 
the tracking requirements associated with levels be- 
low that being measured. 

Due to slight changes in the track trigger require- 



ments over time, the corresponding efficiencies are 
measured in three run ranges. A final efficiency is de- 
termined by taking the luminosity weighted average 
of the results obtained for each run range. The event 
samples used to make these measurements were stud- 
ied to look for possible trigger efficiency dependencies 
on other event variables such as electron isolation, 
number of additional jets in the events, total event 
energy, and electron charge. No dependencies were 
found for these variables, within the statistical uncer- 
tainties of our sample. We did observe a small trigger 
efficiency dependence as a function of the measured 
pseudorapidity of the electron track. We observe a 
small inefficiency for tracks near yydot ~ due to wire 
spacers in the tracking chamber and reduced overall 
charge collection due to the shorter track path length 
through the chamber. However, the effect of this de- 
pendence on our final efficiency results was found to 
be negligible within our measurement uncertainties. 
The final efficiency results for the electron trigger 
path tracking requirements at each trigger level are 
shown in Table IMXD 

In order to measure the total efficiency of our elec- 
tron trigger path, we additionally need to measure 
the efficiencies of the calorimeter cluster requirements 
at each level of the trigger. The requirement of an 
electromagnetic cluster with Et > 8 GeV at Level 1 is 
studied using reconstructed electromagnetic objects 
found in muon-triggered events. We determine the 
highest energy trigger tower associated with each ob- 
ject and check to see if the Level 1 trigger bit cor- 
responding to this tower is turned on in the data. 
We measure a turn-on efficiency of 99.5 % for trig- 
ger towers with a measured electromagnetic energy 
between 8 GeV and 14 GeV and 100% for those mea- 
sured above 14 GeV. The small inefficiency observed 
for towers with measured energies below 14 GeV is 
due to an additional Level 1 requirement placed on 
the ratio of hadronic and electromagnetic energies 
{Ehad/Ecm < 0.05) in towers with energies below 
this cut-off value. The effect of this inefficiency on 
the fully reconstructed electrons in our W ev can- 
didate events is determined by checking how often 
the associated trigger tower with the highest electro- 
magnetic Et has a measured energy below 14 GeV. 
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TABLE XXXL Efficiencies for tracking requirements in high Et electron trigger path. 

Trigger Level Track Requirement Measured Efficiency 

Level 1 Fast Tracker (pT > 8 GeV/c) 0.974 ± 0.002 

Level 2 Fast Tracker {pr > 8 GeV/c) 1.000 ± 0.000 

Level 3 Full Reconstruction {pr > 9 GeV/c) 0.992 ± 0.001 

Combined Level 1 ^ Level 3 0.966 ± 0.002 



We find that less than 1 % of the reconstructed elec- 
trons in our candidate sample {Et > 25 GeV) do not 
have at least one associated trigger tower with Et > 
14 GeV. Based on these numbers, we estimate the 
overall trigger efficiency for the Level 1 electromag- 
netic cluster requirement to be 100 % for the events 
in our candidate samples. 

Additional secondary trigger paths are used to 
measure the efficiencies of the Level 2 and Level 3 
cluster requirements in our default electron trigger 
path. The efficiency of the Level 2 cluster require- 
ment is obtained using events collected with two ad- 
ditional secondary trigger paths that have no Level 2 
selection requirements other than simple prescales. 
The Level 1 and Level 3 trigger requirements in these 
paths are equivalent to in one case those of the de- 
fault path and in another those of the path used to 
collect events for measuring the efficiencies of track 
requirements. The subset of these events that pass 
our full set oi W ^ ev selection criteria are also 
found to satisfy the Level 2 cluster trigger criteria. 
Based on these samples we conclude that the effi- 
ciency of the Level 2 electron cluster requirement is 
100 % for reconstructed electrons also satisfying our 
selection criteria for tight central electrons. Since the 
electron clustering algorithm run in the Level 3 pro- 
cessor farms is nearly identical to that used in of- 
fline reconstruction, we expect candidate events with 
high Et electrons to also satisfy the Level 3 clus- 
ter requirements of our trigger path. However, due 
to slight differences in the calorimeter energy correc- 
tions applied at Level 3 and offline, it is possible that 
we could observe trigger inefficiencies close to the Et 
threshold utilized for Level 3 clusters. To check for 
this inefficiency, we collect events on an additional 
secondary trigger path which is based on the Level 1 
and Level 2 requirements of our default electron trig- 
ger path but no requirements at Level 3 other than 
a simple prescale. We find that all of the events col- 
lected on this path which satisfy our event selection 
criteria also satisfy the Level 3 cluster criteria of our 
default trigger path. Based on this study, the effi- 
ciency of the Level 3 cluster requirement for events 
in our candidate samples is also 100 %. Since we 
do not measure inefficiencies for the cluster require- 



ments of our trigger path at any of the three levels, 
we conclude that the overall efficiency of our default 
trigger path for electrons is completely determined by 
the measured efficiencies of the track criteria given in 
Table KXXll 

As mentioned above we do not have the benefit 
of an equivalent set of secondary trigger paths for 
collecting W ^ candidate events to measure the 
efficiencies of our muon trigger path requirements. 
Instead, we use Z — > /i/i candidate events in which 
both muons satisfy the full set of isolation and iden- 
tification cuts used to define our samples. To avoid 
background events we require that the invariant mass 
of the dimuon pair lies in a tight window around the 
Z boson mass (76 GeV/c^ < M^^ < 106 GeV/c^) and 
that the event has not been identified as a cosmic ray 
by our tagging algorithm. In this sample we know 
that at least one of the two muons in the event sat- 
isfied the muon trigger path requirements and can 
make a measurement of the muon trigger efficiency 
based on the fraction of events in which both muons 
meet the criteria of our trigger path. If we define etrg 
as the single muon trigger efficiency we want to mea- 
sure, then (ctrg)^ is the fraction of events containing 
two triggered muons, and 2 (etrg) (1 — Ctrg) is the frac- 
tion of events with only one triggered muon. There is 
also a remaining fraction of events (1 — Ctrg)^ which 
contain no triggered muons, but these events do not 
make it into our Z ~f fifj, candidate sample. Based 
on these definitions, the number of candidate events 
in our sample in which both muons meet the trigger 
criteria, N2trg divided by the total number of events 
in the sample, iVtot, can be expressed with Eq. 1261 
From this expression we obtain the formula shown in 
Eq. [27] which gives the muon trigger efficiency as a 
function of this fraction F. 

p ^ -^2trg ^ ('^trg)^ .26) 

A^tot (efg)2 + 2(etrg)(l-efg) ^ ' 

To check whether an individual muon in our can- 
didate sample satisfies the requirements of our muon 
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trigger path, we first look at the hits on the recon- 
structed muon stub to determine the position of the 
muon with respect to the 144 Level 1 muon trigger 
towers (2.5 degrees each in 0) defined in the hard- 
ware. We then check to see if the trigger bits corre- 
sponding to each individual requirement of our trig- 
ger path are set for the matched trigger tower. The 
Level 1 requirements of our trigger path include both 
a high pt COT track identified in the fast track- 
ing hardware and a sufficient set of matching hits 
in the muon detector wire chamber(s) along the path 
of the reconstructed muon. Matching CSX scintilla- 
tor hits are additionally required in the region of the 
muon detector between 0.6 and 1.0 in T^dot (CMX re- 
gion). No significant additional trigger requirements 
are made at Level 2 for muon candidates. In order 
to measure the efficiency of the muon reconstruction 
algorithms at Level 3, we use the subset of events 
in the Z — > ^/u candidate sample described above in 
which both muons are found to satisfy the Level 1 
trigger criteria. This restriction is made to ensure 
that we do not double-count correlated inefficiencies 
between the different trigger levels. In addition, we 
require that one of the two muons is found in the 
region of the muon detector between 0.0 and 0.6 in 
r^det (CMUP region) while the other is found in the re- 
gion between 0.6 and 1.0 in T^dot (CMX region). Since 
different Level 3 muon reconstruction algorithms are 
run in these two regions, it is simple to check if both 
or only one of the muons in these events satisfy the 
Level 3 requirements of our muon trigger path. The 
input parameters to our muon trigger path efficiency 
calculations are shown in Table IXXXIII along with 
the final results of these calculations. 



I. Cosmic Tagger Efficiency 

The tagging algorithm used to remove cosmic ray 
events in our W — > fiv and Z — s- fifi candidate sam- 
ples is discussed in Sec. IVII El We measure the frac- 
tion of real events tagged as cosmic rays by this al- 
gorithm for both candidate samples using the corre- 
sponding electron decay mode samples. The tagging 
algorithm is based solely on the hit timing informa- 
tion associated with reconstructed tracks in the COT. 
Since the kinematics of W and Z boson decays into 
electrons and muons are nearly identical, we expect 
that the reconstructed electron tracks in — > et^ 
and Z ^ ee candidate events are a good model for 
the muon tracks in the corresponding decay channels. 
Unlike the muon channels, however, the electron de- 
cay mode candidate samples have a negligible cosmic 
background. Therefore, we obtain a measurement of 



the fraction of real W ^ fiiy and Z fifi signal 
events tagged as cosmic ray candidates directly from 
the observed fraction of events in the corresponding 
electron channels which our algorithm identifies as 
cosmic ray candidates. In order to make the tracks in 
the electron events match as closely as possible with 
those in the muon events, we first apply the muon 
track impact parameter cut described in Sec. IIVI to 
each of the electron candidate tracks in these samples. 
This additional requirement reduces the number of 
events in the W ^ ev candidate sample to 37,070. Of 
the remaining events, only five are tagged as cosmic 
ray candidates by our modified version of the cosmic 
tagging algorithm. The resulting efficiency for a W 
boson decay not to be tagged as a cosmic by our algo- 
rithm is = 0.9999 ± 0.0001. Applying the track 
impact parameter cut to the Z ^ ee sample reduces 
the total number of candidate events to 1,680. Of 
these events, only one is tagged as a cosmic by our 
modified tagging algorithm. The resulting efficiency 
for a Z boson decay not to be tagged as a cosmic by 
our algorithm is ef^, = 0.9994 ± 0.0006. 



J. Over-Efficiency of Z-Rejection Criteria 

The criteria for rejecting Z — > /i/i events in our 
W fiv candidate sample are defined in Sec. IIVI A 
small fraction of real signal events are also removed 
from our candidate sample via this selection criteria. 
We measure the efficiency for signal events to survive 
the Z-rejection cuts directly from simulation. The 
resulting value, 0.9961 ± 0.0001, is determined by the 
number of — *■ fiv candidate events in our simulated 
sample that exclusively fail the .Z-rejection criteria. 

The systematic uncertainty on this efficiency is 
based on a comparison of the shape of the invari- 
ant mass spectrum for the muon plus track candidate 
events rejected solely due to this criteria to the shape 
of the same spectrum obtained from j*/Z — > fj.fi 
simulated events. A comparison of the ratio of re- 
jected events inside and outside the Z-mass window 
(66 GeV/c2 < M^^ < 116 GeV/c^) to that found in 
the J* /Z simulation sample provides a good measure 
whether our rejected events are a relatively pure sam- 
ple of J* /Z decays. Based on this approach, we mea- 
sure an additional systematic of ± 0.17% to apply to 
the Z-rejection efficiency value obtained from simu- 
lation. The final resuh is ez-rcj = 0.9961 ± 0.0017. 
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TABLE XXXIL Efficiencies for high pr muon trigger path. 

Trigger Level Number of Z ~> jj.^ Number of Events Efficiency 

Candidate Events with 2 Muon Triggers 

Level 1 338 293 0.929 ± 0.011 

Level 3 138 137 0.996 ± 0.004 

Combined - - 0.925 ± 0.011 



VII. BACKGROUNDS 

Other physics processes can produce events that 
mimic the signature oiW ^ Iv and Z ^ events in 
our detector. Some processes have similar final-state 
event topologies to those of our signal samples and 
others can fake similar topologies if a non-lepton ob- 
ject within the event is misidentified as an electron or 
muon. In this section, the sources of backgrounds to 
W and Z events are discussed. We separate the back- 
ground sources into three main categories: events in 
which hadronic jets fake leptons; events from other 
electroweak processes; and events from non-collision 
cosmic ray backgrounds. The techniques used to esti- 
mate the contribution to our candidate samples from 
each background source are given in this section along 
with the final estimates. 



A. Hadron Jet Background m W li' 

Extracting the contribution of events to the W — > 
Iv candidate samples in which real or fake leptons 
from hadronic jets are reconstructed in the detec- 
tor is one of the more challenging components of 
our measurements. Real leptons are produced both 
in the semileptonic decay of hadrons and by photon 
conversions in the detector material. Some events 
also contain other particles in hadronic jets which 
are misidentified and reconstructed as leptons. Typ- 
ically, these types of events will not be accepted into 
our W candidate samples because we require large 
event $j, . In a small fraction of these events, how- 
ever, a significant energy mismeasurement does re- 
produce the $rp signature of our samples. Because 
of the large total cross section for hadronic jets in our 
detector, even this small fraction results in a substan- 
tial number of background events in our W candi- 
date samples. These events are particularly difficult 
to model in the simulation since the associated en- 
ergy mismeasurement makes them unrepresentative 
of typical hadronic events. In order to estimate the 
background contribution of these sources to our sam- 
ples, we release the selection criteria on lepton iso- 
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FIG. 23: Eif°/ET versus event $j. for W ^ eu candi- 
dates (no cuts on the lepton isolation fraction variable or 
the event ). The definitions of regions A, B, and C 
which are used in the calculation of the hadronic back- 
ground are provided in the text. 

lation and event $Jrp and use events with low lepton 
isolation and low $Jrp as a model of the background in 
the low lepton isolation and high l^j^ W signal region. 
The contributions in the low and high $!rp regions are 
normalized to the number of events in those regions 
with high lepton isolation based on the assumption 
that there is no correlation between lepton isolation 
and l^j^ in the hadronic background. 

Fig. [23] shows the lepton isolation fraction variable 
plotted against event $]rp for W ^ ev candidates (no 
cuts on lepton isolation fraction or event Tf^j, ). In the 
lepton isolation fraction versus If^j, parameter space, 
we define four regions as follows: 

• Region A: EfjEr < 0.1 and ^j. < 10 GeV 

• Region B: E}^°/Et > 0.3 and < 10 GeV 

• Region C: E}f°/ET > 0.3 and > 25 GeV 

(20 GeV for W ^ /ii/) 

• Region W: E'^°/Et < 0.1 and > 25 GeV 

(20 GeV for W ^iv) 
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Region W is the W ^ £v signal region and the 
others contain mostly hadronic background events. 
The background contribution to the W signal region, 

W 



Nu/^, is estimated using 



(28) 



where N^^^., 



TV^t, iV^t are the number of events in 



regions A, B and C, respectively, as defined above. 
This technique has been pre viously described in [13, 
[37} and more recently in [69| . 

A simple approach would be to assume that all of 
the events in regions A, B, and C are hadronic back- 
ground events. In that case, the observed number 
of data events in each region can be used directly in 
Eq. [28] to extract the hadronic background contribu- 
tion to the W signal region. We further improve our 
estimate, however, by accounting for the fact that 
these regions contain small fractions of signal events 
and events from other electroweak background pro- 
cesses such a.s Z M and W ^ tu in addition to 
hadronic background events. Fig. 1241 shows distribu- 
tions of lepton isolation fraction versus for sim- 
ulated events passing the full set of selection criteria 
(no cuts on lepton isolation fraction or ) for the 
W ^ ev signal, Z ^ ee background, and W ^ tv 
background samples. From these distributions, and 
the equivalent ones for W ^ fj^v candidates, we ob- 
tain modeled event fractions in regions A, B, and C 
relative to the signal region for the signal and other 
electroweak background processes. Based on these 
fractions and our estimates for the relative contribu- 
tions of W ^ ii^, Z ^ U, and W ^ m the sig- 
nal region (see Sec. IVII C|) . we correct the observed 
numbers of events in regions A, B, and C to remove 
the contributions from non-hadronic backgrounds. A 
more accurate estimate of the hadronic background 
in the W signal region is then obtained from Eq. [55] 
using these corrected inputs. Table [XXXTTT] summa- 
rizes both the corrected and uncorrected hadronic 
background estimates for the W signal region ob- 
tained from Eq. [^5] for the W ^ ev and W ^ 
decay channels. 

Since the lower limit on lepton isolation fraction 
and upper limit on event If^rp used to define regions 
A, B, and C are arbitrary choices, we check the ro- 
bustness of our technique for obtaining the hadronic 
background estimates by raising and lowering the 
cuts used to define these regions. We take observed 
changes in the estimated hadronic backgrounds as 
a systematic uncertainty on our measurement tech- 
nique. Fig. shows the dependence of the estimated 
hadronic background contribution to the signal re- 
gion as a function of the lepton isolation fraction and 
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FIG. 24: E^t°/Et versus event for the simulated 
W ^ eu signal, W —> Tf background, and Z — > ee back- 
ground samples. We correct the observed number of data 
events in regions A, B, C to account for events from these 
processes when estimating the hadronic background in 
the W ^ iv candidate samples. 



event Iflj^ values used to define the non-signal regions 
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TABLE XXXIII: Summary of hadronic background event contribution estimates to the W ^ eu and W ^ fiv candidate 
samples. The statistical and systematic uncertainties are indicated. 





Uncorrected 
W ^eu 


Corrected 
W ^eu 


Uncorrected 
W ^ tiu 


Corrected 
W ^ fiu 


Region A 


30023 


26655 


3926 


3575 


Region B 


5974 


5972 


5618 


5615 


Region C 


228 


131 


496 


345 


Region W 


37584 


37584 


31722 


31722 


Hadronic Background 


1146 


587 


346 


220 


Statistical Error 


78 


52 


17 


13 


Systematic Error 




294 




110 


Background Fraction 


3.0 ± 0.2% 


1.6 ± 0.8% 


1.1 ± 0.1% 


0.7 ± 0.^0 



both before and after correcting the number of ob- 
served events in regions A, B, C for W ^ ev signal 
and other background processes. We observe similar 
dependencies using the W ^ iiv candidate sample. 
The background estimate is mostly independent of 
the selection of the lower $rp border for regions A 
and B but does depend on the location of the up- 
per lepton isolation fraction border for regions B and 
C. Although we observe some evidence from simu- 
lated event samples that the observed fluctuations 
are a feature of the hadronic background, we choose 
to use a conservative systematic uncertainty that cov- 
ers the full range of the fluctuations seen in Fig. [25l 
We estimate the range of the observed fluctuations 
to be within 50 % of our central values correspond- 
ing to uncertainty estimates of ± 294 events in the 
W ^ ev candidate event sample and ±110 events in 
the W ^ jiv candidate sample (see Table |X"XXIII|) . 

We make an independent cross-check of the esti- 
mated hadronic background in W eiy events by 
applying a measured rate for jets faking electrons to 
a generic hadronic jet sample. The rate for jets faking 
electrons is measured from events with at least two 
jets with Et > 15 GeV, $rp < 15 GeV, and no more 
than one loose electron. These requirements ensure 
that the input sample has a negligible contribution 
from real W and Z events. From this sample, the jet 
fake rate is defined as the fraction of reconstructed 
jets with Et > 30 GeV that are also found to pass 
the standard set of tight electron cuts. 

We use the Et dependence of the jet fake rate 
in the background estimate. Because of differences 
in the clustering algorithms used for electrons and 
jets, the reconstructed energies of electrons originat- 
ing from hadronic jets are smaller than the recon- 
structed energies of the jets. Scale factors are applied 
to convert the measured jet energies into correspond- 
ing electron cluster energies, and as a consequence the 
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FIG. 25: Dependence of hadronic background estimate 
on the Ei}° /Et and event Ipj, cut values used to define 
the control regions for W — > ev. The results both be- 
fore and after corrections for signal and electroweak back- 
ground contributions to regions A, B, and C are applied 
are shown in triangles and circles, respectively. 



lowest Et bins are not included in the fitted constant 
for the jet fake rate. A significant uncertainty on the 
final background estimate is assigned, however, based 
on the results obtained using different models for fit- 
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FIG. 26: Ipj. distribution for inclusive jet sample 
weighted by measured jet fake rate. The arrow indicates 
the location of the selection cut on used to select 
W —> eu candidate events. 

ting the Et dependence of the jet fake rate. The 
measured fake rate is applied to jets in an inclusive 
jet data sample to determine how often these types 
of hadronic events with fake electrons satisfy the ad- 
ditional selection criteria of our W ev candidate 
sample. Jets in the inclusive sample are required to 
have i;!?^'^^ > 25 GeV where E^''^'"^ is the jet Et 
scaled down to the Et of the fake electron to match 
the electron selection criteria of our sample. The dis- 
tribution in Fig. [26] is the resulting ^t distribution 
for the inclusive jet sample weighted by the jet fake 
rate. The events above the candidate sample ^t cut 
of 25 GeV are integrated giving 800 ± 300 background 
events, consistent with the result obtained using our 
default technique. 



B. Hadron Jet Background in "/* / Z ^ ££ 

Our Z ^ £i candidate samples have smaller over- 
all contributions from background sources than the 
W £v samples. One common background source is 
events in which one or both Icptons are either real or 
fake leptons from hadronic jets. We expect that the 
two leptons in these types of events have no charge 
correlation so that the numbers of opposite-sign and 
same-sign lepton pairs from this source are roughly 
equal. Based on this assumption, we use the num- 
ber of same-sign lepton pair candidates to place an 
upper limit on the number of hadronic background 
events in our opposite-sign dilepton pair candidate 
samples. This approach is only viable for events with 
two central leptons where the lepton charge is taken 
from the reconstructed track. As discussed later in 
this section, the background contribution to Z ^ ee 



events with one central electron and one plug elec- 
tron is measured using a variation of the jet fake rate 
method described previously. 

Since the calorimeter energy associated with muon 
candidates is required to be consistent with a 
minimum-ionizing particle, the probability for a 
hadronic jet to fake a muon is significantly smaller 
than that for an electron. Despite the fact that we 
make no opposite-sign charge requirement on the two 
muon legs in our Z —^ fifj. candidate events, none of 
the 1,785 events in this sample are observed to con- 
tain a same-sign muon pair. Based on finding no 
such events, we estimate a background contribution 
of 0.0 g events from muons produced in hadronic 
jets. 

The number of same-sign events observed in the 
Z — > ee candidate sample needs to be corrected for a 
fraction of real Z ee events that are reconstructed 
as same-sign electron pairs. We observe a total of 
22 events with same-sign electron pairs correspond- 
ing to our sample 1730 Z ee candidate events 
with two central electrons. The invariant mass dis- 
tributions for both the opposite-sign and same-sign 
electron pairs in our candidate sample are shown in 
Fig. [18] and Fig. [27] Both distributions show a peak 
in the Z boson mass window indicating that at least 
some fraction of the same-sign electron pairs are pro- 
duced in Z decays. These events result from decays 
in which one of the electrons radiates a high Et pho- 
ton which subsequently converts in the detector ma- 
terial producing an electron-positron pair. We call 
this type of event a "trident" event. If the track as- 
sociated with the positron from the photon conver- 
sion is matched to the corresponding electron cluster, 
both electrons in the event will be assigned the same 
charge. We remove the contribution of real Z — > ee 
events from the number of observed same-sign elec- 
tron pairs by subtracting the observed number of 
opposite-sign events in the data scaled by the frac- 
tion of same-sign to opposite-sign candidates in our 
simulated samples. The remaining number of same- 
sign electron pair candidate events is then used to 
estimate the background contribution from electrons 
produced in hadronic jets to the opposite-sign candi- 
date sample. 

Using this technique, we estimate 20.4 same-sign 
events from Z decays in the invariant mass window 
between 66 GeV/c^ and 116 GeV/c^. Subtracting this 
estimate from the total number of observed same- 
sign events (22), we estimate the contribution from 
electrons originating from hadronic jets to be 1.6 j^^ g 
where the uncertainty is based solely on the statistics 
of our sample. 

The dominant source of systematic uncertainty on 
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FIG. 27: Reconstructed invariant mass of two central elec- 
trons in Z ^ ee candidate events in data (points) and 
simulation (histogram). This distribution is for events 
in which the electrons are reconstructed with the same 
charge. The distribution for events with two electrons of 
opposite sign is shown in Fig. 1181 The number of events 
in the simulated distributions are normalized so that the 
number of opposite-sign events in the simulated sample is 
equal to the number of opposite-sign events in the data. 
The arrows indicate the location of the invariant mass 
cuts used to select our candidate samples. 



the background contribution from events with elec- 
trons originating from hadronic jets comes from the 
simulation detector material model. The probability 
for an electron to radiate a bremsstrahlung photon 
prior to entering the calorimeter is strongly depen- 
dent on the amount of material in the tracking vol- 
ume. We study the effect of the material model us- 
ing the two previously described samples of simulated 
events generated with ± 1.5 % of a radiation length 
of copper added in a cylinder between the silicon and 
COT tracking detectors. We estimate the system- 
atic uncertainty based on differences in the number 
of same-sign events observed after subtracting the 
predicted number of real Z ~* ee events based on 
the default and modified simulations. The resulting 
systematic uncertainty on our estimate is 5.2 events 
which when added in quadrature with the statistical 
error results in a final background estimate of 1.6^} g. 

This technique outlined above can not be used to 
estimate the background contribution from electrons 
originating from hadronic jets in Z ee candidates 
with one central and one plug electron owing to the 



FIG. 28: Measured tight central and plug electron jet fake 
rates as a function of jet Et- 



undetermined charge of the plug candidate. Instead, 
we estimate the background contamination based on 
a variation of the previously described method us- 
ing measured jet fake rates. In order to measure the 
background contribution to the combined Z — > ee 
sample from hadronic events producing two fake elec- 
trons, we need to measure the jet fake rates for tight 
central, loose central, and plug electrons. We remove 
W and Z boson candidates from the inclusive jet sam- 
ple used to make the fake rate measurements by se- 
lecting events with no more than one loose electron 
and ]^rp < 15 GeV. Based on this inclusive sam- 
ple, the jet fake rates are defined as the fractions of 
central jets reconstructed as either tight or loose elec- 
trons and plug jets reconstructed as plug electrons. 
The measured jet fake rates for reconstructed tight 
central and plug electrons as a function of jet Et are 
shown in Fig. 

As previously mentioned, the reconstructed energy 
of the electrons produced by hadronic jets is smaller 
than the reconstructed energy of the jets themselves. 
To account for these differences, we fit the distri- 
butions of E^°/El^^ to a Gaussian for the jets re- 
constructed as tight central, loose central, and plug 
electrons. The means of the fits are used as scal- 
ing factors to convert raw jet energies into scaled 
electron energies, E^^^°'^. To obtain the background 
contribution of events with two electrons originating 
from hadronic jets to the Z —> ee sample, we apply 
the measured jet fake rates and energy scalings to a 
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generic multi-jet data sample. Events containing ei- 
ther two central jets with E^''^'"^ > 25 GeV or one 
central jet and a plug jet with E^'^^°'^ > 20 GeV are 
used to extract dijet invariant mass distributions to 
model the hadronic background for Z —^ ee. The 
weights assigned to each event in these distributions 
is set equal to the product of the jet fake rates for 
the two jets based on the parameterizations shown in 
Fig. [28l The final weighted dijet invariant mass dis- 
tributions for central-central and central-plug events 
are shown in Fig. 1291 The resulting distributions 
are integrated over the invariant mass window of our 
measurements (66 GeV/c^ < M^e < 116 GeV/c^) to 
obtain an estimate for the number of background 
events in the Z —^ ee candidate sample (after scaling 
upward by the trigger prescale used to collect events 
in the generic multi-jet sample). 
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FIG. 29: Di-jet invariant mass distributions for central- 
central and central-plug events in generic multi-jet data. 
Events are weighted by the product of the measured jet 
fake rates for each jet. The scaled energies of both jets 
must pass the electron Et requirements of our Z — > 
ee candidate sample (25 GeV for central electrons and 
20 GeV for plug electrons). 



As illustrated in Fig. [28l the jet fake rates mea- 
sured as a function of jet Et need to be assigned an 
additional uncertainty based on the assumed shape 
of the fit. We fit the jet fake rate distributions using 
several different functional forms and assign an ad- 
ditional systematic uncertainty of 30 % based on the 
resulting spread in background estimates. Based on 
this technique the measured background contribution 
of events with two electrons originating from hadronic 
jets to the Z ^ ee candidate sample is 2.4 ± 1.0 
central-central and 39 ± 17 central-plug events. The 
estimated number of central-central events is in good 



agreement with the result obtained using the ob- 
served number of same-sign events in our candidate 
sample. Using the central-central background esti- 
mate based on same-sign events and the central-plug 
estimate based on the jet fake rate method, we obtain 
a combined estimate for the background contribution 
of events with two electrons originating from hadronic 
jets of 41 ± 18 events. 



C. Electroweak Backgrounds mW^lv 

Z ^ U events mimic the signature oi W ^ £v 
events in cases where one of the two leptons passes 
through an uninstrumented region of the detector 
creating an imbalance in the observed event Et- 
The W —f signature can also be reproduced by 
W ^ TV events in which the r lepton subsequently 
decays into an electron or muon. Background con- 
tributions from both diboson and tt production pro- 
cesses are negligibly small. 

The contribution of these electroweak background 
sources to our W ^ iv candidate samples are ob- 
tained from simulation. The 7*/Z U and W tv 
simulated event samples are obtained from the equiv- 
alent PYTHIA event generation and detector simula- 
tion used to produce the signal samples (see Sec. IVT) . 
The complete set oi W ^ iv selection criteria as de- 
scribed in Sec.|lV]are applied to the simulated events 
in these samples to obtain the fraction of events from 
each process that satisfy the criteria of our candidate 
samples. Then, based on Standard Model predictions 
for the relative production rates of our signal process 
and the two background processes, we use the esti- 
mated acceptances from simulation to obtain the rel- 
ative contributions of each process to our candidate 
samples. 

The Standard Model predicts equivalent produc- 
tion cross sections for W ^ ev, W ^ and 
W TV, while the Z ^ U production cross sec- 
tions are related to the corresponding W ^ £v cross 
sections via the ratio R defined in Eq. [TJ In order 
to extract the relative contributions of /Z i£ 
events to our W ^ iv candidate samples, an input 
value for R is required. We choose to use the value 
R = 10.67 ± 0.15 for W and Z boson production at 
= 1.96 TeV obtained from the NNLO theoretical 
calculation [l^[3,[l5, 16]. However, to be conserva- 
tive we inflate the uncertainty on the predicted value 
for R based on the CDF Run I measured value of 
R = 10.90 ± 0.43 [H,!!!. The difference in the val- 
ues of i? at = 1-80 TeV and 1.96 TeV is expected 
to be negligible. Based on this prediction, the mea- 
sured value is in good agreement with the theoretical 
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TABLE XXXIV: Estimated W ^ £1^ backgrounds from 
other electroweak production processes. 



Source W ^ ev W ^ fiu 

Background Background 
Z -^U 426 ± 19 2229 ± 96 
W -*TU 749 ± 17 988 ± 24 



value. To account for the current level of experimen- 
tal uncertainty, we add an additional 3.9 % system- 
atic uncertainty to the NNLO prediction resulting in 
a value of i? = 10.67 ± 0.45. 

The relative contributions of ^ tv, Z ££, 
and — !■ rz^ in our W ^ £i> candidate samples are 
estimated based on the above value for R and the 
simulated acceptances for each process. The rela- 
tive acceptances are normalized to the total number 
of events in each candidate sample after subtract- 
ing contributions from non-electroweak backgrounds 
(events with reconstructed leptons originating from 
hadronic jets and cosmic rays). The final background 
estimates for electroweak backgrounds in the W ^ £v 
candidate samples are summarized in Table IXXXIVI 




60 80 100 120 140 160 180 200 

Mee (GeV/c') 

FIG. 30: Reconstructed invariant mass distribution for 
simulated ')* /Z — > ee (open histogram) and 7*/^ tt 
(solid histogram) events satisfying the Z — > ee candidate 
sample selection criteria. 



D. Electroweak Backgrounds in j* /Z ^ ££ 

Several electroweak processes also contribute back- 
ground events to our Z ^ ££ candidate samples. 
Z TT events mimic the Z ^ ££ event signature 
when both r leptons decay into or are reconstructed 
as an electron or muon pair with a reconstructed in- 
variant mass within the mass window of our Z ^ ££ 
measurements. As in the previous section, this back- 
ground is estimated using a simulated Z — > rr event 
sample obtained from the equivalent pythia event 
generation and detector simulation used to produce 
the Z ^ ££ signal samples. The full set oi Z ^ ££ 
selection criteria is applied to the simulated Z — > rr 
and Z ££ samples to determine the relative ac- 
ceptances. Based on the Standard Model prediction 
of equivalent production cross sections for Z — + ee, 
Z iifi, and Z — > rr, the number of Z rr back- 
ground events in each candidate sample is extracted 
using the relative acceptances from the total number 
of events after removing non-electroweak background 
contributions. 

A comparison of the reconstructed invariant mass 
distributions for simulated 7*/Z — > ee and 7*/Z — > 
TT events passing the Z ee selection criteria is 
shown in Fig. 1301 The majority of 7*/Z tt 



events are observed to have a reconstructed invariant 
mass below the mass window used in our measure- 
ments. As a result, the contribution of this back- 
ground source to our candidate samples is small, 
3.7 ± 0.4 events in the Z ^ ee sample and 1.5 ± 0.3 
events in the Z — > /x^ sample. An identical approach 
is used to estimate Z ^ ££ background contributions 
from both top quark and diboson production. The es- 
timated background contributions from each of these 
sources is found in all cases to be less than one event 
and therefore considered to be negligible. 

An additional source of background events to the 
Z — !■ ee candidate sample is W ev events with an 
associated hadronic jet that results in a second recon- 
structed electron within the event. We use our sim- 
ulated W ^ ev sample to estimate the background 
contribution from this source by applying previously 
determined jet fake rates for the hadronic jets in these 
events with scaled Et above the corresponding elec- 
tron thresholds. The relative acceptance of simulated 
W ^ eu events, weighted by the measured jet fake 
rates, and Z ^ ee signal events are used to extract 
the number of background events from this process 
based on the value for R presented in the previous 
section. Once again, the relative acceptances are ap- 
plied to the final candidate sample after subtract- 
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ing the estimated number of background events from 
non-electroweak sources. The estimated number of 
W ev background events in the Z — > ee sample is 
16.8 ± 2.8 events. 



E. Cosmic Ray Backgrounds in W —> fiv 

Energetic cosmic ray muons traverse the detector 
at a significant rate, depositing hits in both muon 
and COT chambers, and in some cases can mimic 
the signatures of our W — > fii^ and Z — > fifi can- 
didate events. A cosmic ray muon passing through 
the detector is typically reconstructed as a pair of in- 
coming and outgoing legs relative to the beam line of 
the detector. The reconstructed muon legs tend to 
be isolated and pass our muon selection criteria. In 
some cases, one of the two cosmic legs is not recon- 
structed due to fiducial and/or timing constraints. 
These events typically satisfy both the .^-rejection 
and criteria of our W fj,i> candidate sample 
due to the lack of an additional track and the result- 
ing transverse momentum imbalance. 

We remove cosmic ray events from our W ^ 
candidate sample using a tagging algorithm based on 
the timing information associated with hits in the 
COT. The algorithm uses a multi-parameter fit over 
the full set of hits left by the incoming and outgoing 
cosmic legs. The leg belonging to the reconstructed 
muon serves as the seed track for the fit. The other 
leg is referred to as the opposite-side track. The al- 
gorithm performs the following steps to determine if 
an event is consistent with the cosmic ray hypothesis: 

• Hits belonging to the seed track are refitted 
with the five helix parameters and a floating 
global time shift, Iq. 

• Based on the best fit values, an incoming or out- 
going hypothesis is assigned to the seed track. 

• The refitted seed track is used to search for the 
hits belonging to the second cosmic leg on the 
opposite side of the COT. 

• If enough hits are found on the other side of 
the COT, a similar fit procedure is performed 
to identify the opposite-side track. 

• If both legs are found, a simultaneous fit is per- 
formed to combine all hits from the seed and 
opposite-side legs into a single helix. 

The final decision of the cosmic tagger depends on 
the quality of the simultaneous fit to the hits on both 
legs. If one leg is recognized as incoming and fits well 



to an outgoing leg on the other side of the detector, 
the event is tagged as a cosmic ray. As described 
in greater detail below, we observe that our tagging 
algorithm identifies most of the cosmic background 
events in our candidate sample. We also find that 
the algorithm tags very few real events as cosmic rays 
(see Sec. ED). 

After removing tagged events from our W fj,v 
sample, we need to estimate the remaining back- 
ground from cosmic rays. This estimate is made by 
searching for hits in the muon chambers on the oppo- 
site side of the reconstructed muon track in our final 
candidate events. These hits are present for a large 
fraction of cosmic ray muons even in cases where the 
second leg is not identified by our algorithm. Since 
the muon chamber hits are not used in the tagging 
algorithm, their presence is unbiased with respect to 
its decision. The distribution for matched hits 
produced by cosmic ray muons with respect to the di- 
rection of the muon candidate track is sharply peaked 
in the region around 180°. These events sit on top of 
a flat event background in Acjj originating from ran- 
dom coincidences between the muon track and unre- 
lated matched hits in the muon chambers. The con- 
tribution of cosmic ray events to the candidate A(j) 
distribution is determined by counting the number of 
events with matched muon chamber hits in a 10° win- 
dow centered on A(/) = 180 ° and subtracting a fitted 
contribution from the flat background. Using this 
approach, we would estimate a cosmic background 
contribution of 54.7 ± 5.0 events in our 31,722 event 
W — > fiiy candidate sample. 

Some of the cosmic ray background events in our 
candidate sample, however, do not have opposite-side 
muon chamber hits due to gaps in the muon detec- 
tor coverage. In order to estimate the total cosmic 
ray background in our candidate sample from the ob- 
served number of events with matched opposite-side 
hits, we apply an acceptance correction based on the 
fraction of W candidate events in which the 

reconstructed muon track points at an active region 
of muon chambers when extrapolated to the opposite 
side of the detector. We extrapolate the 31,722 muon 
tracks in our W ^ candidate events to the oppo- 
site side of the detector and find that 58 ± 30 % point 
at active regions of the muon chambers. Our accep- 
tance correction assumes that the spacial distribution 
of muons originating from cosmic rays is similar to 
that of our W ^ candidate sample. We assign a 
large systematic uncertainty on the measured accep- 
tance to account for the non-uniform spacial distribu- 
tion (most enter from the top side of the detector) of 
cosmic rays and the reconstruction biases associated 
with their entry locations and angles of incidence on 
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the detector. 

To complete the cosmic background measurement 
for our W — > iiiy candidate sample, we also need 
to estimate the contribution of Z — > fifj, events to 
the observed excess of events in the window around 
A(p — 180 °. Z ^ ^fjL events that contain a sec- 
ond reconstructed track passing a loose set of min- 
imum ionizing cuts are rejected from our candidate 
sample via the .^-rejection selection criteria. How- 
ever, a small fraction of muon tracks from Z fj,^ 
events are embedded in jets and fail the loose min- 
imum ionizing cuts or in other cases are simply not 
reconstructed. Since the muons in Z — » /x/it de- 
cays are typically produced in roughly opposite di- 
rections to one another, the non-identified tracks in 
these events can also produce muon chamber hits 
on the opposite side of the one reconstructed muon 
in these events. This background is estimated from 
our simulated Z — > /i/i event sample. Based on this 
sample, we estimate the number of Z — > /i/i back- 
ground events in our W —^ fiv candidate sample with 
matched muon chamber hits in the 10 ° window cen- 
tered on A(j) = 180° to be 35.4 ± 9.1. The uncertainty 
assigned to this background is based on our use of dif- 
ferent techniques for looking at opposite side muon 
chamber hits in data and simulation. 

The final estimate of the cosmic ray background in 
the W ^ jjLv sample, N^'^, is obtained from 

^^bg — 75pp l^yj 

where is the number of ^ iiiy candidate 

events with matched hits in the tight window cen- 
tered on A(p = 180°, N^^^^ is the predicted number 
of Z fifx background events with matched hits in 
the same window, and A^^^ is muon chamber accep- 
tance for muon tracks in — > fiv candidate events 
extrapolated to the opposite side of the detector. Us- 
ing the input values obtained above, we estimate a to- 
tal cosmic background of 33.1 ± 22.9 events for our 
W ^ IJLv candidate sample. 



F. Cosmic Ray Backgrounds in Z —> fifi 

Cosmic rays also contribute to the Z —f fj,fj. can- 
didate sample. The majority of these events are re- 
moved using the cosmic ray tagging algorithm de- 
scribed in the previous section. The remaining cosmic 
ray background is estimated based on the distribu- 
tion of the three-dimensional opening angle between 
the muon tracks in candidate events. The two re- 
constructed muon legs in the cosmic ray background 




Opening Angle (rad) 

FIG. 31: Comparison of the three-dimensional opening 
angle distribution for muon tracks in ^ — > /i/t candidate 
events with the same distribution from simulated events. 
The simulated distribution is scaled to match the data in 
the region below 2.8 radians. 



events are typically back-to-back with opening angles 
at or near 180°. The residual background is esti- 
mated by fitting the opening angle distribution for 
data events in the region below 2.8 radians (assumed 
to be background free) to the same distribution for 
simulated Z ~* ^ii events. The output of the fit 
is a scale factor for the distribution from simulation 
which is also applied in the region above 2.8 radians. 
The number of scaled simulation events with an open- 
ing angle greater than 2.8 radians is compared to the 
number of data candidate events in the same region. 
The observed excess in data over simulation is taken 
as our estimate of the cosmic ray background. Using 
this technique, we estimate a total of 12 ± 12 cosmic 
ray background events in our Z /i/t candidate sam- 
ple where the quoted uncertainties are based on the 
statistics of our data sample. A comparison of the 
opening angle distribution between data and scaled 
simulation is shown in Fig. 1311 

G. Background Summary 

Based on the information presented in the preced- 
ing sections, the estimated background contributions 
to the W — > ev, W ^lv, Z ^ ee and Z /i/i 
candidate samples are summarized in Table IXXXVI 

VIII. RESULTS 

Using the measured event counts, kinematic and 
geometric acceptances, event selection efficiencies. 
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TABLE XXXV: Summary of background event estimates for the W ^ Iv and Z ^ U candidate samples. 



Background source W ev W ^ ni' Z — > ee Z ^ fijj, 



Multi-jet 


587 ± 299 


220 ± 112 


41 ± 18 




z ^ee 


426 ± 19 


2229 ± 96 






Z TT 


negl. 


negl. 


3.7 ± 0.4 


1.5 ± 0.3 


W -^TU 


749 ± 17 


988 ± 24 


negl. 


negl. 


W ^fM 






16.8 ± 2.8 


negl. 


Cosmic rays 


negl. 


33 ± 23 


negl. 


12 ± 12 



Total 1762 ± 300 3469 ± 151 62 ±18 13 ± 13 



background estimates, and integrated luminosities for 
our candidate samples, we extract the W and j* /Z 
boson production cross sections multiphed by the lep- 
tonic (e and /i) branching ratios. We also determine 
a value for the ratio oi W ^ iv to Z ^ M cross sec- 
tions, i?, taking advantage of correlated uncertainties 
in the two cross section measurements which cancel 
in the ratio. To test for lepton universality, we use 
the measured ratio oiW^lv cross sections in the 
muon and electron channels to extract a ratio of the 
W ^ coupling constants, (/^/^e- Then, based on 
the assumption of lepton universality, we increase the 
precision of our results by combining the production 
cross section and cross section ratio measurements 
obtained from the electron and muon candidate sam- 
ples. The resulting combined value of R is used to 
extract the total decay width of the W boson, r(T/F), 
and the W leptonic branching ratio, Br{W iiy), 
which are compared with Standard Model predic- 
tions. The measurement of r(M^) is also used to con- 
strain the CKM matrix element Vcs- 



TABLE XXXVL Summary of the input parameters to 
the W ^ Iv cross section calculations for the electron 
and muon candidate samples. 



W 



W 



Ai-bck 

Aw 
ew 

jcdt (pb-^; 



31722 
3469 ± 151 



37584 

1762 ± 300 
0.2397 to;oo42 0.1970 
0.749 ± 0.009 0.732 ± 0.013 

72.0 ± 4.3 72.0 ± 4.3 



0.0024 
0.0031 



and 



(Tw ■ Br{W fiu) = 2.722 ± 0.0l5{stat.) 

± 
± 



0.061^ 

0.163{lum.) nb 



(31) 



We compare our measurements to a recent NNLO to- 



W 



Cross Section 



The cross section cr{pp W) times the branching 
ratio Br{W £1/) is calculated using Eq. [5] given in 
Sec. m The measurements of the required input pa- 
rameters for the electron and muon candidate sam- 
ples are described in the previous sections and sum- 
marized in Table IXXXVII Based on these values, we 
obtain 



aw ■ Br{W -> eiy) = 2.771 ± O.OU{stat. 

± 
± 



0.056 V 

0.166(/wm.) nb 



tal cross section calculation for y/s = 1.96 TeV 
which utihzes the MRST 2002 NNLL PDF set 
\6T\ . The resulting predicted W ^ Iv cross section is 
2.687 ± 0.054 nb, which agrees well with our mea- 
sured values in both lepton channels. The uncer- 
tainty on the predicted cross section is mostly due 
to PDF model uncertainties derived from the MRST 
error PDF sets. We also perform an independent 
calculation of the uncertainty on the total W tv 
cross section originating from uncertainties in the 
PDF model using the method described in Sec. IVl 
Based on this method, we obtain a consistent 1.3 % 
uncertainty based on the MRST error PDF sets and 
a 3.9 % uncertainty based on the CTEQ6 error PDF 
sets. 

Distributions of electron E^, muon p^^ 



(30) event Iprp , and W transverse mass (Mt = 



59 



u4000 
> 

O3500 
^3000 

•4—' 

1 2500 
lU 

2000 
1500 
1000 
500 





. W^|iv Data 
□Sum 

W^liv MC 
nyVZ^IiLi MC 
□W^xv MC 
■QCD 



20 40 60 80 100 

Pj (GeV/c) 



> 

(§5000 

CM 

^4000 
> 

UJ 

3000 



2000 
1000 




20 



40 



. e V Data 

□Sum 
Vy^ e V MC 
y/Z^ e e MC 

□W^ X V MC 

■QCD 



60 



80 100 
ET(GeV) 



FIG. 32: Muon pt (left) and electron Et (right) distributions for W ^ lu candidate events in data (points). The 
solid lines are the sum of the predicted shapes originating from the signal and background processes weighted by their 
estimated contributions to our candidate samples. The separate contributions originating from the signal and each 
individual background process are also shown. 



^2[Et$t - {Ex$T,x + Ey$T y )]) for events in 
our W ^ iv candidate samples are shown in Figs. [32] 
- [331 The data distributions are compared against a 
sum of the predicted shapes of these distributions for 
the W ^ £iy signal and each contributing background 
process {Z ^ U, W ^ tv, and hadronic jets). The 
predicted shapes are obtained from our simulated 
event samples except in the case of the background 
arising from hadronic jets, which is modeled using 
events in the data containing non-isolated leptons 
that otherwise satisfy the W ^ £v selection criteria. 
In the sum, the predicted shape obtained for each 
process is weighted by the estimated number of 
events in our W ^ lu candidate samples originating 
from that process (see Table IXXXV[) . In the case 
of the IfJrp distribution, we remove the selection 
cut on the If^rp variable to include events with low 
]prp in the comparison and highlight the significant 
background contribution from hadronic jets in this 
region. 



B. 'y* /Z ^ U Cross Section 

Similarly, the cross section <7{pp J* /Z) times 
the branching ratio Br{'-f* /Z ££) is calculated us- 
ing Eq. [6| given in Sec. [H The measurements of the 



required input parameters for the electron and muon 
candidate samples are described in the previous sec- 
tions and summarized in Table IXXXVlH Based on 
these values, we obtain 

a^*/z ■ Br{-f*/Z -* ee) ^255.8 ± 3.9(siat.) 

± liisyst.) 
± 15.3{lum.) pb 
(32) 

and 

a^,/z ■ Br{j*/Z ~* fi^i) = 248.0 ± 5.9(stat.) 

± ?:2(s2/sf.) 

± 14.8(;Mm.) pb . 

(33) 

These measurements are the cross sections for dilep- 
tons produced in the mass range 66 GeV/c^ < 
Mu < 116 GeV/c^ where both 7* and Z bo- 
son exchange contribute. A correction factor of 
F = 1.004 ± 0.001 determined from a NNLO da/dy 
calculation, phozpr [H, [H, [H , using MRST 2002 
NNLL PDFs [63, is needed to convert these mea- 
sured cross sections into those for pure Z boson ex- 
change over the entire dilepton mass range; the mea- 
sured cross sections need to be multiplied by F. We 
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FIG. 33: Event distributions for W ^ iv candidate events in data (points). The selection requirement on event 
$j. has been removed to include candidate events with low . The solid lines are the sum of the predicted shapes 
originating from the signal and background processes weighted by their estimated contributions to our candidate 
samples. The separate contributions originating from the signal and each individual background process are also 
shown. The arrows indicate the location of the event selection criteria used to define our candidate samples. 



compare the corrected cross sections for pure Z bo- 
son exchange to the recent NNLO total cross section 
calculations for = 1.96 TeV The Z ^ U 

production cross section predicted by these calcula- 
tions is 251.3 ± 5.0 pb, which is in good agreement 
with the corrected, measured values obtained in both 
lepton channels. The uncertainty on the predicted Z 
boson production cross section is also primarily due 
to uncertainties in the PDF model derived from the 
MRST error PDF sets. Our independent estimates 
for these uncertainties using the method described in 
Sec. |V] are a consistent 1.2% uncertainty based on the 
MRST error PDF sets and a somewhat larger 3.7% 
uncertainty based on the CTEQ6 error PDF sets. 



TABLE XXXVII: Summary of the input parameters to 
the y* /Z ££ cross section calculations for the electron 
and muon candidate samples. 



■ ee 



"f'/Z fin 



Az 

/ £dt (pb- 



4242 1785 
62 ±18 13 ± 13 

n -^lao +0 0039 r, -.nQn +0.0027 
U.OiOZ -0.0041 U.±0»z _0.0033 

0.713 ± 0.012 0.713 ± 0.015 
72.0 ± 4.3 72.0 ± 4.3 



Ratio ofW^ivtoZ 



U 



Fig. [T7] and Fig. [TH] show the invariant mass dis- 
tributions for events in our Z ^ It candidate sam- 
ples. The data distributions are compared against 
predicted shapes from our simulated Z ^ U event 
samples. The predicted shapes are normalized to 
the total number of events in the candidate samples. 
In making these comparisons, we ignore background 
processes which account for less than 1 % of the events 
in these samples (see Table IXXXVp . 



Precision measurements of the ratio oiW^lv 
to Z ^ M production cross sections, R, are used to 
test the Standard Model. The Standard Model pa- 
rameters r(M^) and Br{W Iv) can be extracted 
from our measured values of this ratio and are sen- 
sitive to non-Standard Model processes that result 
in additional decay modes for the W boson. A new 
high-mass resonance which decays to either W ov Z 
bosons could also have a direct effect on the measured 
value for R. 
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FIG. 34: Transverse mass (Mt) distributions for W ^ Iv candidate events in data (points). The solid lines are the sum 
of the predicted shapes originating from the signal and background processes weighted by their estimated contributions 
to our candidate samples. The separate contributions originating from the signal and each individual background 
process are also shown. 



TABLE XXXVIII: Summary of the input parameters to 
the R calculations for the electron and muon candidate 
samples. 



jy-obs 
Aw 



37584 
1762 ± 300 
4242 
62 ± 18 



31722 
3469 ± 151 
1785 
13 ± 13 



1.3272 ± 0.0109 0.7066 ± 0.0068 
0.952 ± 0.011 0.974 ± 0.010 
1.004 ± 0.001 1.004 ± 0.001 



The ratio of cross sections can be expressed in 
terms of measured quantities: 



R 



1 

F 



Ajohs 



Ajbck 



^bck 



Aw 



(34) 



where F is the correction factor for converting the 
measured ^*/Z ^ £i cross section into the cross sec- 
tion for pure Z boson exchange and the other pa- 
rameters are as defined for the W and Z production 
cross section measurements. The integrated luminos- 
ity terms in the W and Z cross section calculations 



along with their associated uncertainties cancel com- 
pletely in the R calculation, allowing for a signifi- 
cantly more precise measurement of the ratio than is 
possible for the individual cross sections. In addition, 
we take advantage of many correlated uncertainties in 
the event selection efficiencies and kinematic and ge- 
ometric acceptances of our W and Z candidate sam- 
ples which cancel in the ratios Az/Aw and ez/ew- 
For example, uncertainties on the acceptances arising 
from the PDF model are significantly smaller for the 
ratio of the Z ^ ££ and W ^ acceptances than 
for either individual acceptance. The calculation of 
Az/Aw and tz/^w for our electron and muon can- 
didate samples and the treatment of the correlated 
uncertainties in these ratios are discussed in Sees. IVl 
and I VII The event counts and background estimates 
for the W ^ £iy and Z U candidate samples are 
the same as those used in the individual cross section 
calculations. Table IXXX Villi summarizes the input 
parameters used to calculate R using the electron and 
muon candidate samples. Substituting these values 
into Eq. [Ml we obtain 



Re = 10.79 ± Q.YI{stat) ± 0.16(st/sl 



(35) 



and 



i?^ = 10.93 ± 0.27(stot.) ± 0.18(s7/si.) . (36) 
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Based on the calculations of the production cross 
sections for W ^ and Z ^ U provided by 
m El, H E [H, the expected value for R at 
y/s = 1.96 TeV is 10.69. To obtain an accurate esti- 
mate for the uncertainty on this prediction, we need 
to account for correlated uncertainties in the individ- 
ual cross section predictions. The error originating 
from PDF model uncertainties has the largest con- 
tribution to the total uncertainty. We estimate the 
magnitude of this contribution using the previously 
defined method in Sec. |V] and obtain a 0.45 % un- 
certainty based on the MRST error PDF sets and a 
larger 0.56 % uncertainty based on the CTEQ6 er- 
ror PDF sets. We also need to account for the ef- 
fect of additional uncertainties in the values of the 
electroweak parameters and CKM matrix elements 
used in the cross section calculations. We estimate 
these uncertainties using the MS NNLO total cross 
section calculation, zwprod We have up- 

dated the calculation code to incorporate the CTEQ 
and MRST PDFs and variations of the electroweak 
parameters and CKM matrix elements. We obtain 
an uncertainty of 0.15 % for the az calculation and 
0.40% for the aw calculation. The larger uncertainty 
associated with the aw calculation is due primar- 
ily to experimental uncertainties on the CKM ma- 
trix values. To be conservative, we add the larger 
PDF model uncertainty (0.56%) in quadrature with 
the individual cross section calculation uncertainties 
(0.15% and 0.40%) to obtain a combined uncertainty 
on the prediction for R of 0.70%. The resulting pre- 
diction, 10.69 ± 0.08, agrees with the measured val- 
ues of R in both lepton channels. 



D. fj,-e Universality in W Decays 

Stringent tests of lepton universality at LEP pro- 
vide strong evidence for lepton universality in Z ^ i£ 
production. We make a similar test for lepton univer- 
sality inW^lv production by extracting the ratio 
oi W —> couplings, gfi/ Qei from the measured ra- 
tio of the W ^ jiv and W ^ ev cross sections. The 
W ^ iiy couplings are related to the measured ratio 
U of the cross sections, defined as 

gyy ■ BrjW ^xv) ^ T{W ^ H ^ 

~ aw ■ Br{W ev) T{W ev) gl ' 

(37) 

As in the case of the R measurements described in 
the previous section, many of the uncertainties asso- 
ciated with the individual cross section measurements 
cancel in the ratio. Table lXXXIXI summarizes the un- 
correlated uncertainties between the two cross section 



TABLE XXXIX: Uncertainties on the measured ratio of 
W HI' and W ev cross sections, U . 



Category 


Uncertainty 


Statistical Uncertainty 


0.0075 


Acceptance Ratio: 




Simulation Statistics 


0.0019 


Boson pt Model 


0.0001 


PDF Model 


+0.0003 
-0.0004 


Pt Scale and Resolution 


0.0018 


Et Scale and Resolution 


0.0034 


Material Model 


0.0072 


Recoil Energy Model 


0.0010 


Efficiency Ratio: 




Uncorrelated 


0.0199 


Backgrounds: 




Hadronic 


0.0043 


Electroweak 


0.0030 


Cosmic Ray 


0.0008 



measurements that contribute to the overall uncer- 
tainty on gpi/ge- The uncertainties due to the PDF 
model cancel almost completely in the ratio. The 
major remaining contributions to the systematic un- 
certainty come from the uncorrelated event selection 
efficiencies for the electron and muon candidate sam- 
ples. Since these efhciencies are measured directly 
from Z ^ U candidate events in the data, the asso- 
ciated uncertainties will decrease as additional data 
are analyzed. In this sense, the remaining uncertainty 
on gp,/ ge is primarily statistical in nature and can be 
reduced with larger data samples. Using the input 
parameters to our W iv cross section measure- 
ments, we obtain 

^ = 0.991 ±0.012 . (38) 

Using Eq. [37] and the current world average of exper- 
imental results for Br{W fiv) = 0.1057 ± 0.0022 
and Br{W ev) = 0.1072 ± 0.0016 20], the ex- 
pected value of g^j./ ge is 0.993 ± 0.013 which is in 
good agreement with our measured value. 

E. Combined Results from the Electron and 
Muon Channels 

Since our measurement of g^/ge supports the con- 
clusion of lepton universality in W ^ £v produc- 
tion, we proceed to combine our measurements of 
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TABLE XL: Uncertainty categories for the inclusive W 
cross section measurements. These values are absolute 
contributions to aw in pb. The uncertainties in the elec- 
tron and muon channels for each category are treated as 
either 100% correlated (1.0) or uncorrelated (0.0). 



Category 


Electron Muon Correlation 


Statistical Uncertainty 


14.3 


15.3 


0.0 


Acceptance: 








Simulation Statistics 


3.6 


3.9 


0.0 


Boson pt Model 


1.2 


1.0 


1.0 


PDF Model 


36.9 


35.4 


1.0 


Pt Scale and Resolution 


0.8 


5.6 


1.0 


Et Scale and Resolution 


9.5 


0.0 


0.0 


Material Model 


20.2 


0.0 


0.0 


Recoil Energy Model 


6.8 


9.4 


1.0 


Efficiencies: 








Vertex zo Cut 


11.7 


11.5 


1.0 


Track Reconstruction 


11.1 


10.9 


1.0 


Trigger 


2.9 


32.7 


0.0 


Lepton Reconstruction 


11.1 


19.7 


0.0 


Lepton Identification 


24.1 


23.8 


0.0 


Lepton Isolation 


9.4 


9.7 


0.0 


Z-rejection Cut 


0.0 


4.6 


0.0 


Cosmic Ray Algorithm 


0.0 


0.3 


0.0 


Backgrounds: 








Hadronic 


23.1 


10.8 


1.0 


Z -^U 


1.5 


9.2 


1.0 


W ->TV 


1.3 


2.3 


1.0 


Cosmic Ray 


0.0 


2.2 


0.0 



the W Iv and Z U production cross section 
measurements in the electron and muon channels to 
increase the overall precision of these results. We 
also combine our measurements of i?e and i?^ to de- 
termine a precision value for R which is used to test 
the Standard Model. 



1. Combination of the Cross Sections 

We use the Best Linear Unbiased Estimate 
(BLUE) [tOi Izil method to combine measurements 
in the electron and muon channels. For the W ^ Iv 
measurements, we identify twenty categories of un- 
certainties, several of which are correlated in the elec- 
tron and muon channels. Table IXLI lists these cate- 
gories and summarizes the raw contribution of each 
(in pb) to the W cross section measurements in the 



electron and muon channels. Based on the informa- 
tion in this table, we combine the measurements in 
the two lepton channels and obtain 

aw ■ Br{W £v) = 2.749 ± 0.010(siai.) 

± Om^isyst.) 
± 0.165(km.) nb, 

(39) 

which has a precision of 2.0 %, not including the un- 
certainty associated with the measured integrated lu- 
minosity of our samples. The uncertainty on luminos- 
ity is not included in the calculation of the combined 
value. 

The combination of the Z ^ It cross section mea- 
surements in the electron and muon channels is based 
on the same procedure. In this case, we identify sev- 
enteen categories of uncertainties, some of which are 
correlated between channels. Table IXLII provides a 
list of these categories and summarizes the raw con- 
tribution of each (in pb) to the Z cross section mea- 
surements in the electron and muon channels. The 
additional acceptance for forward electrons in the 
plug calorimeter modules reduces the statistical un- 
certainty associated with the Z cross section mea- 
surement in the electron channel, which thus has a 
larger weight in the final combination. The combined 
result is 

(7^,./^ • Br(77Z ^ «) = 254.9 ± 3.3(stoi.) 

± 4.6(s?/sf.) 
± \h.2{lum) pb, 
(40) 

which has a precision of 2.2 %, not including the un- 
certainty associated with the measured integrated lu- 
minosity of our samples. As discussed previously, the 
combined cross section given here is the cross sec- 
tion for dileptons in the mass range 66 GeV/c^ < 
M« < 116 GeV/c^ including contributions from both 
7* and Z boson exchange. In order to convert the 
measured cross section into a cross section for pure 
Z boson exchange over the entire mass range, one 
must multiply the measured value by the correction 
factor presented earlier, F = 1.004 ± 0.001. 

A comparison of the predictions from [13, [13, [3, 
[H [H for aw ■ Br{W iv) and az ■ Br{Z U) 
as a function of the pp center of mass energy, i?cM , 
with our measured values and other experimental re- 
sults [H, m, [ill are shown in Fig. [35l 
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FIG. 35: W ^ and Z U cross section measurements as a function of the pp ce nter-of-mass energy, Ecm- The 
solid lines correspond to the theoretical NNLO Standard Model calculations from [l^. [l^. [isl. [l^ . 



2. Combination of the R measurements 

The same BLUE method is also used to combine 
our measurements of Re and R^. For our cross sec- 
tion ratio measurements we identify fifteen categories 
of uncertainties, some of which are correlated be- 
tween our measurements in the electron and muon 
channels. Table IXLIII lists these categories and sum- 
marizes the raw contribution of each to the i?e and 
i?^ measurements. Since most of the uncertainties 
related to efhciency factors are uncorrelated in the 
electron and muon channels, the corresponding un- 
certainties are combined into a single net uncertainty 
for uncorrelated efficiencies. The exception is the 
uncertainty on COT track reconstruction efficiency 
which is 100 % correlated between the two channels. 
The combined result is 

R ^ 10.84 ± 0.15{stat.) ± O.U{syst.) (41) 

which is precise to 1.9%. 



F. Extraction of Standard Model Parameters 

As previously discussed, the precision value for R 
obtained from the combination of our measurements 



in the electron and muon channels can be used to 
measure various Standard Model parameters and in 
the process test the predictions of the model. The 
ratio of cross sections can be expressed as 

a{pp^W)T{W T{Z) 

a{pp -> Z) T{Z -> U) r{W) ■ ^ ^ 

Using the precision LEP measurements for T(Z — > 
££)/r{Z) at the Z pole mass and the NNLO calcula- 
tion ofa{pp w)/<j{pp -> z) by 

we extract the Standard Model parameter Br{W 
£v) = T[W iu)/T{W) from Eq.gH using our mea- 
sured value of R. Using the Standard Model pre- 
diction for V{W iv), we also make an indirect 
measurement of r(T/F) and based on this value place 
a constraint on the CKM matrix element Vcs- 



1. Extraction of Br (W ^ Iv) 

The required parameters to extract Br{W Iv) 
from our measured R value using Eq. are the pre- 
dicted ratio of W and Z production cross sections 
and the measured value of Br{Z ££) = r(Z 
££)/T{Z). The value of aw/c^z obtained from the 
NNLO calculations provided by [H, [H, [ll, [H, [3 is 
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TABLE XLI: Uncertainty categories for the inclusive Z 
cross section measurements. These values are absolute 
contributions to az in pb. The uncertainties in the elec- 
tron and muon channels for each category are treated as 
either 100% correlated (1.0) or uncorrelated (0.0). 



Category 



Electron Muon Correlation 



Statistical Uncertainty 


3.93 


5.87 


0.0 


Acceptance: 








Simulation Statistics 


0.61 


1.01 


0.0 


Boson 'p^ IVIodel 


0.16 


0.19 


1.0 


PDF Model 


1.96 


4.94 


1.0 


Pt Scale and Resolution 


0.10 


0.13 


1.0 


Et Scale and Resolution 


0.67 


0.00 


0.0 


Material Model 


2.45 


0.00 


0.0 


Recoil Energy Model 


0.00 


0.00 


0.0 


Efficiency: 








Vertex Cut 


1.08 


1.04 


1.0 


Track Reconstruction 


1.42 


1.98 


1.0 


Trigger 


0.17 


2.05 


0.0 


Lepton Reconstruction 


1.43 


1.24 


0.0 


Lepton Identification 


3.39 


3.48 


0.0 


Lepton Isolation 


1.21 


1.77 


0.0 


Cosmic Ray Algorithm 


0.00 


0.15 


0.0 


Backgrounds: 








Hadronic 


1.10 


0.08 


1.0 


Z ^ TT 


0.02 


0.04 


1.0 


W ~*eiy 


0.17 


0.00 


1.0 


Cosmic Ray 


0.00 


1.76 


0.0 



3.3696 with associated relative uncertainties of 0.0056 
coming from the PDF model and 0.0043 coming from 
electroweak and CKM matrix parameters used in 
the calculations (see Sec. I VIII C[) . The experimen- 
tal value of Br{Z U) = 0.033658 ± 0.000023 as 
measured at LEP is taken from [20|. 

When extracting Br{W iv) from R, it is impor- 
tant to consider correlated uncertainties in the ratio 
of predicted cross sections and the ratio of accep- 
tances, Az/Aw, used in the measurement of R. In 
a sense, we measure Br{W tv) by equating i?phys 
to i?meas, where 



and 



phys 



A^obs 



aw Br{W -> tv) 
Br{Z ie) 



ATbck 



Azez 



Awew 



jyohs _ ^bck 



(43) 



(44) 



TABLE XLII: Uncertainty categories for the R measure- 
ments. The uncertainties in the electron and muon chan- 
nels for each category are treated as either 100 % corre- 
lated (1.0) or uncorrelated (0.0). 



Category 


Electron 


Muon 


Correlation 


Statistical Uncertainty 


0.1748 


0.2659 


0.0 


Acceptance Ratio: 








Simulation Statistics 


0293 


0.0472 


0.0 


Boson Pt Model 


0.0020 


0.0044 


1.0 


PDF Model 


0.0701 


0.0836 


1.0 


Pt Scale and Resolution 


0.0012 


0.0167 


1.0 


Et Scale and Resolution 


0.0184 


0.0000 


0.0 


Material Model 


0.0322 


0.0000 


0.0 


Recoil Energy Model 


0.0267 


0.0377 


1.0 


Efficiency Ratio: 








Uncorrelated 


0.1204 


0.0999 


0.0 


Track Reconstruction 


0.0169 


0.0437 


1.0 


Backgrounds: 








Hadronic 


0.0437 


0.0399 


1.0 


Uncorrelated Electroweak 


0.0089 


0.0094 


0.0 


Correlated Electroweak 


0.0057 


0.0369 


1.0 


Cosmic Ray 


0.0000 


0.0689 


0.0 


Correction Factor, F 


0.0107 


0.0109 


1.0 



Then, 

Br{W iv) 



ATobs 



N^'' ez 



ew 



jyhck 



Azcrz 
Aw aw 



Br{Z ^ «) . 



(45) 



The ratio of the acceptance times the cross section for 
Z and W bosons on the right-hand side of Eq.|45]is af- 
fected by uncertainties in the PDF model. To account 
properly for correlations between the PDF uncertain- 
ties associated with each of these four quantities, we 
independently calculate a PDF model uncertainty for 
the quantity contained within the parentheses using 
the method described in Sec.|Vl The measured PDF 
model uncertainties on this quantity are found to be 
slightly larger than for those on Az/Aw alone (0.9% 
versus 0.6% in the electron channel and 1.0% versus 
0.8% in the muon channel). These correlated uncer- 
tainties are separately accounted for in our extraction 
of Br{W iv) from the measured value of R. We 
obtain 



Br{W tv) = 0.1082 ± 0.0022 



(46) 
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Standard Model 
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(includes Run I results) 
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Br(W^lv) 



Standard Model 
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(includes Run I results) 



(from R) 
CDF II (e+i^) 
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FIG. 36: Comparison of our measured value of Br(W —> 
iu) with previous hadron collider measurements [2a. l27l . 
|32| |. the current world average of experimental results [20| |, 
and the Standard Model expectation [20I ]. 



FIG. 37: Comparison of our measured value of FfWO with 
previous hadron collider measurements [2^ . [271 . [32 . [7p . 
|73 |. the current world average of experimental results |20(]. 
and the Standard Model expectation [20l |. 



where the uncertainty contributions are from R 
(± 0.00212), the predicted ratio of cross sections 
(± 0.00047), and the Z ^ U branching ratio 
(± 0.00007). The Standard Model value for this pa- 
rameter is 0.1082 ± 0.0002, and the world average 
of experimental results is 0.1068 ± 0.0012 [loj, both 
of which are in good agreement with our measured 
value. A summary of Br{W — > Iv) measurements is 
shown in Fig. 1361 



2. Extraction ofV{W) 

An indirect measurement of r(T4^) can be made 
from our measured value of Br{W — *■ £v) using the 
Standard Model value for the leptonic partial width, 
T(W iv). We use the fitted value for T{W ^ Iv) 



of 226.4 ± 0.4 MeV [20]. Based on this value, we 
obtain 

T{W) = 2092 ± 42 MeV (47) 

which can be compared to Standard Model predic- 
tion of 2092 ± 3 MeV ^ and the world average of 
experimental results, 2118 ± 42 MeV [l^j. A sum- 
mary of r(l/l^) experimental measurements is shown 
in Fig. 1371 Our indirect measurement is in good 
agreement with the fit [l^] and the theoretical pre- 
diction as well as other measurements in literature. 

An alternative approach for obtaining r(VF) is to 
first use the predicted values for both T{W tv) 
and r(Z U) to extract a ratio of the total widths, 
T{W)/T{Z), from the measured value of R. The pre- 
cisely measured value of V{Z) from the LEP experi- 
ments (2495.2 ±2.3 MeV ^fl]) is then used to extract 
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a value for r(M^). Using this approach we obtain 
T{W) 



nz) 



0.838 ±0.017 



(48) 



for the ratio of total widths, which can be com- 
pared to the Standard Model prediction of 0.8382 ± 
0.0011 [23|. Based on the measured value of T{Z) we 
obtain 



T{W) = 2091 ± 42 MeV 



(49) 



where the uncertainty on the measured value for 
T{Z) makes a negligible contribution to the total un- 
certainty. Since the measurement of T{Z) is inde- 
pendent of the measurement of the branching ratio 
Br{Z ££), both extracted values of r(VF) are in- 
dependent to some degree. 



3. Extraction of Vcs 

In the Standard Model the total W width is a sum 
over partial widths for leptons and quarks where the 
latter subset involves a sum over certain CKM matrix 
elements [2011: 



r 



w 



3r^ 



(i + ^^ 

V TT 



1.409(— )^ 

TT 



12.77(— )^ 

TT 



IK, 



[no top] 



(50) 



Only the first two rows of the CKM matrix contribute 
as decays to the top quark are kinematically forbid- 
den. Thus the relevant CKM matrix elements are 
Kd, Vus, Vcd, Vcs, Kb, and Kb- Of these, Ks con- 
tributes the largest uncertainty. We use the indirect 
measurement of r(TK) from our measured value of 
Br{W iv) as a constraint on Vcs based on world 
average measurements of all the other CKM matrix 
elements and find 



muon decay channels based on 72 pb~^ oipp collision 
data at y/s = 1.96 TeV. We calculate the ratio of the 
W and Z cross sections, R, in each lepton channel 
and combine them to obtain a value which is precise 
to 1.9%. The precision will improve when more data 
are analyzed. From this ratio we extract the leptonic 
W branching ratio, the W width, the ratio of the 
W and Z widths, and constrain the CKM matrix 
element Ks- A summary of extracted quantities is 
given in Table |2£LIIIl 
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IK. 



0.976 ±0.030 



(51) 



using Us = 0.120 and = 226.4 MeV i20|. Our 
measured value is more precise than the direct mea- 
surement at LEP, \Vcs\ = 0.97 ± 0.11 (zl, IH, but 
not as precise as the combined value from LEP and 
Run I at the Tevatron, |KJ = 0.996 ± 0.013 fzl. 



G. Summary 



We have performed measurements for the W and 
Z boson production cross sections in the electron and 
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TABLE XLIII: Standard Model parameters extracted from the measured ratio of W and Z production cross sections, 
R. 



Quantity Our Measurement World Average SM Value 

Br{W iv) 0.1082 ± 0.0022 0.1068 ± 0.0012 0.1082 ± 0.0002 
T{W) in MeV 2092 ± 42 2118 ± 42 2092 ± 3 

V{W)/T{Z) 0.838 ± 0.017 0.849 ± 0.017 0.838 ± 0.001 

V^s 0.976 ± 0.030 0.996 ± 0.013 N/A 

g^/ge 0.991 ± 0.012 0.993 ± 0.013 1 
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